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Abstract 

We  study  the  diffraction  of  a  two-dimensional  Gaussian  wavepacket  through  a  rectangular  aperture  in  a  finite 
potential  wall  (one  slit  e^^riment).  For  wavepacket  with  inadent  wavevector  Icq  samfying  the  diffiacdon 
condition,  Icq  =  2k/w,  w  being  the  slit  width,  the  near  field  (Fresnel-like)  di^atmon  pattern  behind  the  slit  can 
be  clearly  seen  for  snuQl  time  duration  (  <  Oilps).  At  later  time  steps,  the  difiracted  beam  is  fiagmented  into 
lobes  (perpendicular  to  the  Section  of  incidence  of  the  wavepacket)  as  a  result  of  the  multiple  refiections  of 
the  wavepacket  inside  the  slit  (assumed  to  be  of  finite  thickness).  At  later  time,  no  far-field  Fraunhofier 
diffraction  pattern  is  observed  in  our  numerical  simulations. 
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Introduction 

With  the  advent  of  sophisticated  growth  technit^ues  such  as  Molecular  Beam  ^ita»  and  Metal  Oiganic 
Chemical  Vapor  Deposition,  there  has  been  an  mcreased  theoretical  interest  in  various  quantum  medianical 
tuimeling  problems  mduding:  (1)  the  resonant  tunneling  of  electrons  throu^  double  bamer  heterostiuctures, 
a  problem  of  primary  importance  in  asserting  the  switching  time  of  resonant  hmneling  device  (RTD)  (Hwge 
and  Stovnena  1989;  Goldberg  and  coworkers,  1967),  (2)  the  electron  propagation  throu^  narrow  ballistic 
constriction  (defined  by  a  spht-gate)  in  the  two-dimensional  electron  gas  of  a  GaAs-ALGa^.j^As 
heterostructure.  The  latter  followed  the  recent  eiqierimental  discovery  (Van  Wees  and  ooworkers,  1988)  that 
the  conductance  of  such  constrictions  increases  in  a  sequence  of  steps  of  height  2e^/h  (at  sufficiently  low 
temperature).  More  recently,  the  possibility  of  using  narrow  split-rate  for  transistor  a^lications  h^  been 
suggested  by  Kriman  and  coworkers  (1988)  in  their  newly  proposed  QUADFET. 

Both  RTD’s  and  QUADFETs  have  potential  hi^-speed  device  applications  within  the  terahertz  regime 
^aura  and  Stovneng,  1989;  Bandyopadhyay  and  coworkers,  1989).  On  one  hand,  the  fast  switching  response  of 
KTD^s  has  been  w  .dely  investigated  using  pui  ely  ore  dimensional  analysis  of  wavepacket  propagatibn  throu^ 
double  barrier  heterostructures  (Hauge  and  Stovneng,  1989;  Goldbeig  and  coworkers,  1967).  In  practice 
however,  electrons  are  injected  from  3D  contacts  and  tunneling  through  the  quantum  well  of  the  kTD  is 
characterized  by  two-dimensional  dynamics.  On  the  other  hand,  in  view  of  their  potentiality  for  high-speed 
applications  (such  as  in  the  QUAD^T),  there  is  now  a  call  for  a  transient  analysis  of  quantum  transport 
through  narrow  ballistic  constrictions  to  supplement  the  two  dimensional  steady-state  analysis  completed 
recently  by  several  authors  (Szafer  and  Stone,  1989;  Kirezenow,  1989  and  references  therein). 

Hereafter,  we  limit  our  numerical  simulations  to  the  diffraction  of  a  gaussian  wavepacket  through  a  narrow 
aperture  in  a  potential  wall  of  finite  height  The  thickness  d  and  len^  w  of  the  slit  are  assumed  to  be  100/\  and 
200A  respectively,  which  corresponds  to  an  aspect  ratio  a  =  d/w  equal  to  0.5.  The  choice  of  such  a  narrow  slit 
was  imposed  to  make  the  problem  tractable  numerically  without  having  to  use  an  excessive  number  of  grid 
points.  However,  the  time  evolution  of  the  diffracted  beam  shows  interesting  features  which  could  stilTbe 
present  in  the  electron  difiiraction  through  a  more  realistic  mlit  rate  realizable  with  present-day  technology 
(Thornton  and  coworkets,  1986;  Zheng  and  coworkers,  198o).  The  use  of  electron  diffiraction  through  a 
split-gate  was  recently  proposed  by  Krunan  and  coworkers  (1988)  for  interesting  device  applications. 

The  Numerical  Approach 

Recently,  Ancilatto  and  coworkers  (19891  have  developed  a  method  to  solve  the  multi-dimensional 
Schroedinger  equation  based  on  the  Chewchev  polynomial  eimansion  of  the  time  evolution  operator.  Each 
term  in  this  eimansion  is  calculated  using  FFT  computations.  The  subsequent  effort  to  calculate  the 

Juantum-mecnanical  wavefunction  scales  rou^ly  as  MN^*  In  N,  where  M  is  the  number  of  terms  in  the 

!hebychev  expansion,  d  is  the  number  of  space  dimensions  and  N^  b  the  total  number  of  grid  points.  One 
drawback  of  this  approach  is  that  it  requires  the  use  of  uniformly  spaced  grid  points.  Anomer  drawback  of  the 
Chebychev  expansion  of  the  time-evolution  operator  is  that,  being  a  truncated  serIv.o  expansion,  it  d«(es  not 
conserve  the  unitarity  of  the  time-evoltition  operator  vnl?'s  the  Uuncated  series  converges.  As  a  consequence, 
energy  and  norms  of  liit  wavefunction  need  not  be  conserved  at  any  time  step. 
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In  Ais  paper,  we  used  an  alternative  approach  to  solve  the  time-dependent  Schroedinger  equation  based  upon  a 
finite-oifterence  solution  procedure  to  solve  a  set  of  coupled  equations  governing  the  real  and  imaginary 
components  of  the  wavefunction.  The  use  of  Crank-Nicolson  otfferenang  scheme  insures  conservation  of  the 
norm  of  the  wavefunction  at  all  times.  The  resulting  coupled  difference  equations  are  then  solved  usii^  a  block 
alternating  direction  implicit  (ADI)  technique  following  the  scalar  ADI  development  of  Eiouglas  and  Gunn 
(1964).  Recently,  a  similar  algorithm  has  been  used  by  Barker  (1989)  to  study  the  wavepacket  propagation 
through  Aharanov-Bohm  rings.  The  technique  can  be  used  with  non-uniform  grid  spaangs,  allows  for  an 
expliat  time-dependence  of  the  potential  energy  profile  and  can  readily  be  extended  to  include  thepresence  of 
a  spatially  varying  effective-mass,  of  an  extern^ magnetic  field  and  to  tnree  dimensional  configuraUons  (Cahay 
and  coworkers,  1989). 


Numerical  Examples 

As  an  illustrative  example,  we  consider  the  diffraction  of  a  two-dimensional  Gaussian  wavepacket  through  a 
narrow  constriction  such  as  the  one  shown  in  the  upper  left  frame  of  Fig.  1.  The  subsequent  fra.mes  show  the 
time-evolution  of  the  wavepacket  impinging  on  a  100  A  wide  potential  barrier  of  height  0.06eV  containing  a  200 
A  wide  slit  where  the  potential  is  assumed  to  be  zero.  The  slit  is  disposed  symmetric^ly  with  respect  to  the 
center  of  the  simulation  domain,  a  square  of  dimensions  3,000A  x  3,000A.  Our  numeric^  simulations  were 
performed  using  a  non-uniform  grid  spacing  with  69  and  79  grid  points  in  the  x  and  y  directions  respectively.  In 
an  actual  split-gate,  the  slit  can  actually  be  varied  from  a  few  hundred  to  a  few  thousand  angstroms  while 
sweeping  the  gate  voltage.  The  potential  in  the  2D  electron  gas  is  also  different  from  the  sn^  potential  wall 
assumeoin  the  present  work.  We  will  comment  on  this  point  in  our  conclusion  section.  The  initial  wavepacket 
is  assumed  to  be 


i^(x,y,t=o) 


exp  [  ikoX ]  exp  [ - { ( x-x^)  "  -  ( y 'Yq )  M/ 2o  ^ ) 


(1) 


Where  (xq.Yq)  are  the  coordinates  of  the  wavepacket  center  and  o  is  pqual  to  lOOA.  Finally,  the  electron 
wavevector  1^  is  0.0094  A*‘  (with  this  value  of  Icq,  a  free  electron  (m  =0.06701^)  travels  a  length  of  l.OOOA  in 
0.6  ps).  The  average  kinetic  energy  of  an  electron  in  the  state  (1)  is  about  10  m^  for  the  values  of  the 
parameters  listed  above  and  is  therefore  about  l/6th  of  the  barrier  height  A  Fermi  energy  of  about  10  meV 
corresponds  to  a  ^ical  2D  electron  gas  sheet  density  of  -  3  x  lO^*  cm’A  In  Fig.  1,  the  frames  arc  taken  at 
various  physical  time  steps  equal  to  a  multiple  of  0.05  ps  (which  is  equal  to  10  times  the  computational  time 
step).  We  nave  plotted  the  logarithmic  contour  plots  (logm)  of  the  probability  density  |  ^(t)]  * .  For  = 
O.OO^A*^  the  electron  De  Broglie  wavelength  (X  =  2x/0  is  about  660A  and  therefore  bi^er  than  the  slit 
width.  Consequendy,  even  though  part  of  the  wavepacket  (a  35%)  is  transmitted  on  the  right-hand  side  of  the 
wall  at  time  ta0.25pi  no  diffracuon  lobes  are  detectable  in  the  transmitted  waveform.  For  an  electron  with  an 
incident  energy  of  45  meV  (such  an  electron  is  far  in  the  tail  of  the  Fermi  distribution  in  the  previous  exanmle), 
Iq,  =  0.0265A'^,  which  is  close  to  the  diffraction  condition,  kp=2itAv.  In  that  case,  as  can  be  seen  in  Fig.  2,  tnere 
is  actually  a  diffrartion  of  the  electron  wavepacket  through  tne  slit  as  indicated  by  the  existence  of  three  distinct 
lobes  in  die  contour  plots  of  the  charge  density  profiles  in  the  immediate  vicinity  (on  the  right)  of  the  slit  at 
time  t  equal  0.15  ps.  However,  for  later  time  steps,  the  diffracdon  lobes  actually  smear  into  a  main  one  which 
(for  this  specific  case)  precludes  the  actual  observaUon  of  a  Fraunhoffer-like  diffracdon  pattern  far  from  the 
slit  In  fact,  the  diffracdon  condition,  Iq,  =  2xAv,  is  only  met  by  a  small  fracdon  of  the  various  Fourier 
components  of  the  wavepacket  incident  on  the  slit  As  a  result,  the  diffracdon  pattern  cannot  be  as  sharp  as  in 
the  case  of  an  incident  plane  wave,  which  is  the  idealistic  situadon  equivalent  to  the  one  used  iu  opdcs  to  study 
light  diffraction  througn  a  narrow  slit  More  numerical  simuladons  involving  modificadon  of  the  shape  of  the 
wavepacket  potential  walls  and  slit  dimensions  and  extension  of  the  simulation  domain  need  to  be  performed 
however  before  determining  if  an  appropriate  set  of  parameters  can  eventually  lead  to  a  Fraunhoffer-like 
diffraction  pattern  far  from  the  slit  The  overall  shape  of  the  wavepacket  behind  the  slit  changes  drasdcally 
while  varying  the  direcdon  of  incident  wavevector  1^.  This  is  illustrated  in  Fig,  3  where  the  wavepacket  with  a 
kinedc  energy  equal  to  10  meV  was  assumed  to  be  incident  at  a  45*  angle  on  the  slit  This  example  stresses  the 
importance  of  collimadng  the  electron  beam  in  order  to  observe  a  diffracdon  pattern  with  different  lobes 
beyond  the  slit  This  point  was  already  stressed  by  Kriman  and  coworkers  (1988)  in  their  steady  state  analysis  of 
the  QUADFET. 

The  width  of  the  potential  barrier  below  the  split  gate  is  of  primary  importance  in  determining  the  diffracdon 
pattern  behind  the  split-gate.  As  a  result  of  the  finite  width,  the  electron  wavepacket  suffer  muldple  reflecdons 
between  the  potci*Lial  v-lls  defined  by  the  slit  The  resuldng  spread  in  transit  dmes  introduces  aadidonal 
structure  in  tne  diffracted  electron  beam.  For  instance,  the  time  frame  t  =  035  ps  in  Fig.  2  above  clearly  shows 
the  presence  of  successive  maxima  in  the  transmitted  part  of  the  wavepacket  in  tne  direcdon  perpendicular  to 
the  potential  wall.  Finally,  additional  lobes  can  be  clearly  seen  in  the  reflected  part  of  the  wavepacket,  i.e..  on 
the  left  side  of  the  slit.  Such  "Back  diffraction'  would  probably  also  e^t  while  considering  realistic 
split-gate  and  it  could  be  eliminated  by  appropriately  collimating  the  incident  electron  beam. 


Conrlusion<; 

P'-eliminary  results  describing  the  time-evolution  and  diffraaion  of  a  wavepacket  through  a  narrow  slit  in  a 
finite  potential  wall  have  been  reported.  In  order  to  save  computadonal  time,  the  size  of  the  slit  was  chosen  to 
be  smaller  than  the  actual  split-gate  realized  with  present  day  technology  (Thornton  and  coworkers  (1986);  Van 
Wees  and  coworkers  (1986).  Several  intere'ting  features  in  the  time  evolution  of  the  wavepacket  have  however 
been  pointed  out:  ( 1 )  A  Fresnel-like  diffraction  pattern  is  seen  at  early  time  once  the  diffraction  requirement  is 
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satisfied,  i.e.,  kQ=2TtAv  (k^  being  the  wavepacket  incident  wavenumber  and  w  being  the  slit  width);  (2)  Due  to 
the  finite  thickness  of  the  slit  (in  our  simulation,  the  aspect  ratio  of  the  slit  is  equal  to  OJ),  the  main  diffracted 
lobe  is  found  to  be  fragmented  in  the  direction  peiwndicular  to  the  slit  aperture  as  a  result  of  the  multiple 
reflections  through  the  slit  of  finite  length;  (3)  No  r  raunhoffer-like  diffraction  pattern  is  seen  far  from  the  slit 
even  when  the  slit  width  is  comparable  to  the  electron  De  Broglie  wavelength. 

One  should  emphasize  that  the  Fresnel-like  diffraction  pattern  observed  in  our  simulations  is  an  artifact  of  the 
use  of  infinitely  sharp  corners  in  the  potential  energy  profile.  A  more  realistic  simulation  would  require 
self-consistent  calculations  of  the  piotential  energy  profile  below  a  split-gate  such  as  reported  by  Laux  and 
coworkers  (1988).  Even  with  appropriate  collimaUon  of  a  nearly  monochromatic  inciaent  electron  beam,  the 
rounding  of  the  comers  below  the  split-gate  on  the  scale  of  the  wavelength  of  the  incident  electron  could  be  the 
dominant  source  of  suppression  of  the  diffraction  pattern  at  low  temperatures. 


3,000X 


l,500X 


Figure  1.  Time  evolution  of  a  Gaussian 
wavepacket  moving  through  a  slit  (in 
the  direction  Indicated  by  the  arrow 
k^) .  The  incident  kinetic  energy  of 

the  wavepacket  is  10  meV.  No 
diffraction  pattern  with  different 
lobes  is  seen  beyond  the  slit.  This 
figure  should  be  compared  with  Figure  2 
in  ihc  elect.- or*  iftsorOyAxc 

wavelength  is  approximately  equal  to 
the  slit  width. 
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»  =  0.15  ps 


t  =  0.25  ps 


t  =  0.35  ps 


Figure  2.  Same  as  Figure  1  for  a  wavepacket  with  an  incident  energy  equal 
to  45  roeV.  Note  the  formation  of  diffraction  lobes  after  crossing  the  slit 
which  can  clearly  be  seen  at  time  t  -  O.lSps.  At  time  t  “  0.35ps,  the  main 
transmitted  lobe  of  the  diffracted  wavepacket  is  clearly  segmented  into 
various  fragments.  These  are  due  to  multiple  reflections  of  the  wavepacket 
inside  the  slit  which  has  finite  length.  As  a  result,  different  parts  of  the 
wavepacket  are  delayed  further  in  time  leading  to  successive  transmitted  lobes 
Spurious  effect  due  to  reflection  at  the  boundaries  of  the  simulation  domain 
are  only  seen  for  t  -  0.5ps. 


t  =  O  05  ps 


t  =  0.30  ps 


1,500  A 


3,000 A 


l,500A 


Figure  3.  Same  as  Figure  1  while  shooting  at  an  angle  of  AS 
on  the  slit  The  incident  kinetic  energy  is  10  meV. 
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ABSTRACT 

Mathenatlcally  well  posed  problems  have  been 
formulated  and  numerical  schemes  have  been 
implemented  for  solving  the  equation  of  motion  of 
the  density  matrix.  Equivalent  first  order  systems 
of  equations  have  been  derived  for  the  equation  of 
motion  of  the  density  matrix  and  these  systems  are 
solved  rather  than  the  second  order  equation  of 
motion.  The  advantage  of  the  approach  Is  that 
suitable  forms  of  boundary  conditions  can  be 
Inferred  from  the  characteristic  directions  of  the 
first  order  system.  Further,  the  first  order  system 
can  be  solved  rapidly  by  a  'method  of 
characteristics'  based  numerical  scheme  that  has 
been  developed.  Typical  solutions,  which  arc 
discussed,  require  about  20  seconds  of  Cray-XKP  CPU 
time  and  demonstrate  the  accuracy  and  efficiency  of 
the  solution  procedure. 

IMTRODUCTION 

Frensley  [1]  was  the  first  to  suggest  that  solutions 
to  the  equation  of  motion  of  the  density  matrix  In 
the  coordinate  representation  could  be  an  Important 
cool  In  understanding  the  physics  of  quantum  well 
aiul  quantum  barrier  devices.  In  [2]  we  described 
typical  results  obtained  from  a  solution  algorithm 
we  have  developed  for  Che  equation  of  motion  of  the 
density  matrix.  In  this  paper  we  describe  essential 
features  of  the  solution  algorithm. 

SOUmOK  PROCEDURE 

The  equation  of  motion  of  the  density  macrxx  Is  (!]: 

a  1*  fa*  a*  ' 

—  — - p(x,x',t) 

at  2ia  lax*  ax'* 

(1) 

i 

- Iv(x,t)-v(x',t)]p(x,x',t) 

A 

where  p  Is  the  density  matrix,  m  Is  Che  effective 
mass,  and  v  is  Che  potential  energy,  Ve  focus 
attention  here  on  obtaining  time -Independent 
solutions  of  equation  (1). 

For  convenience  of  solution  and  determining  suitable 
forsM  of  boundary  conditions,  equation  (1)  Is 
written  In  the  form  of  a  coupled  first  order  system 
of  equations: 

iA  fap  ap  1 

U  (X,X')  +  —  I —  +  -  •«  0  (2) 

2a  lax  ax'J 


ap  rau  au  i  i 

—  + - +  -  [v(x)-v(x')]p-0  (3) 

*at  lax  ax'J  a 
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Equation  (2)  Is  the  definition  of  u(x,x')  and 
equation  (3)  Is  equation  (1)  written  In  terms  of  u 
and  p.  The  characteristic  directions  for  the  system 
of  equations  (2)  and  (3)  ace 

q  •=  (X  +  x')/2=constzmt  (4) 

f  —  (X  -  x')/2«constant  (5) 

In  terms  of  the  characteristic  directions  q  and  , 
equations  (2)  and  (3)  can  be  written  as 

iA  ap 

U(x,x')  + - -  0  (6) 

2n  at! 

ap  au  1 

—  +  —  +  -  [v(x)-v(x')]p*=o  (7) 

at  af  A 


Suitable  boundary  conditions  for  equations  (6)  and 
(7)  are  the  specification  of  p  and  u  along  the 
boundary  x'-O  and  the  specification  of  u  along  the 
boundary  x-L,  where  L  Is  Che  length  of  the  device. 
Along  the  boundary  x-0,  p  Is  specified  as  Che 
complex  conjugate  of  p(x,0),  since  p  Is  hermlclan, 
and  u  Is  computed  from  the  outgoing  characteristic 
equation  (7).  Along  the  boundary  x-L,  p  is  computed 
from  the  outgoing  characteristic  equation  (6). 

An  alternate  system  of  first  order  equations  can  be 
formulated  In  terms  of  the  current  matrix. 


lA  ap 

j(x,x') 

+ - -  0 

(8) 

2m  ar 

ap  aj 

i 

—  +  — 

+  -(v(x)-v(x') ]p=0 

(9) 

at  aq 

ft 

Equations  (8)  and  (9)  have  the  same  characteristic 
directions  q  and  f  as  equations  (2)  and  (3). 
Suitable  boundary  conditions  for  equations  (8)  and 
(9)  are  the  specification  of  p  and  J  along  the 
boundary  x'-O  and  the  specification  of  p  along  the 
boundary  x-L.  Along  the  boundary  x-0,  J  Is 
specified  as  the  complex  conjugate  of  J(x,0)  and  p 
Is  computed  from  the  outgoing  characteristic 
equation  (8).  Along  the  boundary  x-L,  j  Is  computed 
from  the  outgoing  characteristic  equation  (9).  Both 
sets  of  the  first  order  system  of  equations, 
equations  <6)- (7)  and  equations  (8) -(9),  are  useful 
In  applications  since  they  allow  different  forms  of 
boundary  conditions.  Both  sets  of  equations  can  be 
solved  by  the  same  numerical  procedure. 

The  overall  solution  procedure  consists  of  solving 
the  first  order  system  of  equations  as  an  initial 
boundary-value  problem  s  .artlng  from  conditions 
along  the  line  x'-O  and  marching  to  the  line  x'-L 
using  the  method  of  characteristics.  A  character¬ 
istic  net  for  the  equation  of  motion  of  the  density 
matrix  can  be  constructed  a-prlorl  from  grid  points 
of  a  uniform  square  grid.  A  discrete  form  of 
equations  <6)  and  (7)  can  be  written  on  this  grid  as 
(Figure  1]. 


P(i, j)-p(i-l, j-1) 


iA 


tu] 


Ai; 


0  (10) 


apl  U(i+1, j-l)-U(i, j) 

-  +  - 

atjav  Af 

i 

+  -  (V{x)-v(x') )av[p]av  “  0 

A 


(11) 


where  [‘Jav  represents  an  average  over  the  grid 
cell.  Depending  upon  the  fom  of  averaging  chosen, 
equations  (10)  and  (11)  fom  a  systen  of  2x2  block 
tridiagonal  or  block  diagonal  algebraic  equations 
that  can  be  easily  solved  at  x'-J  from  known  values 
at  X*— J-1.  Thus,  the  solution  procedure  can  be 
marched  from  boundary  conditions  at  x'-O,  In  steps 
along  x' ,  to  x'-L.  A  similar  procedure  can  be 
utilized  for  equations  (8)  and  (9).  Self- 
consistency  Is  Included  in  the  analysis  by  iterating 
the  solution  of  the  density  matrix  equation  with  the 
solution  of  Poisson's  eqtiatlon  to  convergence. 

RESULTS 

The  equilibrium  distribution  of  potential  and 
density  for  a  double  barrier  structure  of  (3]  were 
computed.  Computations  were  carried  out  on  a 
uniform  grid  of  300  x  300  points.  Maxwellian 
distribution  of  density  was  assumed  at  the  boundary. 
Flgiire  2  shows  the  self -consistent  potential 
distribution  in  the  structure.  Figure  2  also  shows 
the  density  distribution  with  local  accuaKila^ion  of 
charge  In  the  quantum  well  due  to  tunneling.  Figure 
3  Is  a  plot  of  the  density  matrix  on  the  (x,x') 
plane,  showing  Maxwellian  distribution  near  the 
boundaries,  depletion  of  charge  In  the  vicinity  of 
the  barriers,  and  accumulation  in  the  qvtantum  well. 
For  brevity,  only  one  typical  result  baa  been 
presented  here.  The  solution  procedure  has  been 
used  to  compute  distributions  In  a  variety  of 
structures. 


CONCLUSIONS 

The  ability  to  compute  solutions  of  the  equation  of 
motion  of  the  density  matrix  provides  a 
computational  tool  to  examine  fundamental  issues 
associated  with  sMsoscoplc  structures.  Of 
partl''ular  relevance  Is  that  the  time  Independent 
computations  took  about  20  CPU  seconds  on  a  Cray 
X-MP  and  thereby  provide  a  new  efficient  means  of 
exploring  the  physics  of  quantum  structures. 
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Figure  1.  Discretization  of  the  Equations 


Figure  2.  Potential  and  Density  vs.  Distance 


Flgture  3.  Density  Matrix  on  the  (x,x')  Plane 
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I .  INTRODUCTION 


In  the  Effective -Mass  Superlattice  (EMSL)  proposed 
by  Sasaki  [1,2],  the  effective-mass  of  electrons  (or 
holes)  is  changed  periodically  and  the  conduction 
(for  electrons)  or  valence  band  (for  holes)  is 
supposedly  aligned.  This  eliminates  the  potential 
discontinuities  between  the  respective  superlattice 
layers  which  are  present  in  doping  and  compositional 
superlattices .  Previous  work  [3]  neglecting  space - 
charge  effects  has  shown  the  threshold  voltage  for 
negative  differential  resistance  (NDR)  in  EMSL  to  be 
much  lower  (==10mV)  and  the  current  density 
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(10“^ -10^  A/cm^)  much  higher  than  in  resonant 
tunneling  diodes  (RTD) .  Hereafter,  we  report 
preliminary  self-consistent  calculations  of 
space -charge  effects  in  EMSL  and  compositional 
superlattices  under  zero  bias  condition. 


II.  THEORY 


We  use  the  technique  described  in  ref.  [4]  to 
calculate  quantura-raechanically  the  electron  charge 
density  profile  in  various  types  of  superlattices 
(see  fig.  1).  Under  the  assumption  of  ballistic 
transport,  the  Schroedinger  equation  throughout  the 
entire  structure  is  solved  using  the  scattering- 
matrix  approach  [5]  for  each  electron  impinging  from 
the  contacts,  with  the  usual  boundary  conditions  for 
plane-wave  solutions,  requiring  everywhere 
continuity  of  the  wave  function  and  its  first 
derivative  divided  by  the  electron  effective  mass. 
Once  the  Schroedinger  equation  is  solved  for 
electrons  impinging  from  the  contacts,  the  electron 
density  is  then  calculated  for  the  two  streams  of 
electrons  incident  from  the  left  and  right  contacts. 
The  calculation  for  electron  density  and 
electrostatic  potential  are  then  performed 
iteratively,  until  the  electrostatic  potential 
converges  to  a  final  solution. 

M - SUPERLATTICE - ► 

I - 1  I  I  I  I  I - 1  I - < 

I  I  I  I  I  I  I  )(  I 

I  I  I  I  I  I  I  I  I  ■ 

. I  I  I  I 

_ I  I _ I  I - 1  I _ I  I _ • 

-  I  —  .——.  . — 

^  Buffer ^l^pacer^  jmjjm,  j  012*011  {  ^2  |  oi,  joi2joi, 

( m, )  I  ( 012  j  I  ,  I  I  j  I  !  (  {  '  ( '”2 )  i  ( ) 


Fig.  1.  Spatial  variation  of  the  conduction 
band  edges  of  the  device  studied  in  the  text.  The 
dashed  (full)  curve  describes  a  compositional  (EMSL) 
super lattice  respectively.  Only  the  buffer  regions 
are  assumed  to  be  doped  (lO^^cra"^). 
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III.  NUMERICAL  EXAMPLES 


The  technique  outlined  above  was  applied  to 
calculate  self-consistently  the  electron  density  and 
conduction  band  energy  profile  of  two  types  of 
superlattices  whose  configuration  are  shown  in 
figure  1.  The  first  superlattice  consists  of  five 
regions  of  100,72^30,28^^0.86^0.14  with  width 
41A  and  effective  mass  ra]^=0.039mo  separated  by 
InP  regions  of  width  29.34A  and  effective  mass 
m2'0.073mQ.  The  current-voltage  characteristic 
of  this  EMSL  was  calculated  in  ref.  [3]  neglecting 
space-charge  effects.  As  shown  in  ref.  [1],  this 
special  choice  of  layer  thicknesses  leads  to  equal 
energy  gaps  between  the  EMSL  energy  subbands. 

This  EMSL  is  sandwiched  between  two  undoped  In? 
spacer  layers  (50A)  contacting  two  highly  doped 
(lO^^cm"^)  InP  buffer  regions  (500A).  The 
second  superlattice  is  a  compositional  superlattice 
obtained  by  replacing  the  Ino,  72^30, 28^^0. 86^0. 
regions  in  the  EMSL  described  above  by 
Alo .  3Gao ,  yAs/GaAs  respectively  (m2^=0.09i9mo  and 
m2=0.067mo).  The  conduction  band  discontinuity 
between  AIq, 3Gao , 7As/GaAs  is  taken  to  be  0.209  eV 
(see  fig.  1) . 

Figure  2a)  shows  the  self-consistent  charge 
density  profile  in  the  two  different  superlattices 
under  zero  bias  condition.  Also  shown  for 
comparison  is  the  charge  density  profile  obtained 
without  iteration,  i.e. ,  using  the  conduction  band 
energy  profile  such  as  shown  in  figure  1.  As  a 
result  of  the  doping  gradients,  internal  contact 
potentials  are  important  (while  different)  in  both 
types  of  superlattices  as  can  be  seen  in  figure 
2b) .  This  leads  to  an  upward  shift  of  the 
conduction  band  energy  profile  in  both  superlattices 
resulting  in  a  substantial  decrease  of  the  charge 
density,  when  calculated  self-consistently,  from 
their  non  self-consistent  values.  As  shown  in  fig. 
2a) ,  for  identical  doping  concentration  in  the 
buffer  regions,  the  charge  density  is  about  one 
order  of  magnitude  higher  in  the  EMSL  than  in  the 
compositional  superlattice. 
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Fig.  2.  Self-consistent  and  non  self-consistent 
(a)  charge  density  profiles,  (b)  conduction  band 
energy  profile  in  an  EMSL  and  compositional 
superlattice  (see  text  and  fig.  1)  with  identical 
doping  concentration  in  the  buffer  regions.  The 
dashed-dotted  curves  show  the  self -consistent 
results  for  a  structure  in  which  the  compositional 
superlattice  has  been  replaced  by  undoped  GaAs 
(n'^-n-n'^)  . 
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IV.  CONCLUSIONS 


Our  preliminary  (zero  bias)  self-consistent 
calculations  have  shown  that,  in  superlattices  of 
similar  dimensions,  space-charge  effects  are  more 
important  in  EMSL  than  in  compositional  super¬ 
lattices  .  The  resulting  increased  overall 
capacitance  of  EMSL  should  therefore  have  a 
detrimental  effect  on  their  potentiality  for 
high-speed  device  applications  [3].  Indeed, 
inclusion  of  space -charge  effects  in  RTD  has  been 
shown  to  reduce  (up  to  501)  their  calculated  peak 
current  and  peak- to -valley  ratio  and  to  shift  to  a 
higher  voltage  their  NDR  [4]  (to  an  extent 
increasing  with  the  importance  of  space  charge) . 

A  careful  investigation  of  the  superiority  of  EMSL 
over  RTD  for  achieving  low-power  and  ultrafast 
bistable  switches  must  therefore  await  a  detailed 
self-consistent  calculation  of  their  current -voltage 
characteristic . 
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INTRODUCTION 

Frendqr  [1]  w«  the  first  to  suggest  that  solutions  to  the  equatioa  of  motioa  of  the 
densi^  ma^  in  the  coordinate  representation,  p(xXt)  as  an  important  tool  in 
understandiiu  of  the  piques  of  quantum  well  and  quantum  bamer  devices.  A 
mathematicaUy  well  pos^  problem  has  been  implemented  here  for  solving  tbi< 
equation. 

In  one  dimension,  the  equation  of  motion  of  the  density  matrix  is: 

-(1/3)  tv(*,t)-v(xst)  J/»(x,x^t) 

We  have  solved  the  ab^eouation  using  two  separate  solution  procedures.  The 
first  procedure,  for  udii A  only  time  indq)endent  results  are  disqrsser^  invdhies 
the  method  of  characteristics  m  wind  we  obtain  equivalent  first  order  tystems. 
For  example,  in  <me  case,  the  coupled  first  order  sj^tem  h  obtr^ied  after  defining 
the  oraeitt  matrix: 

(2)  'i(X,X')  +  (l9/2n)l3p/ax-3p/3x'}-0 
vdiich  transforms  the  density  matrix  equation  into: 

(3)  3^/3t+3  j/3x+3  j/3x^+(i/*)  (V(X)-V(X')  ]/>-0 


w  cnir^  along  the  boon^ries  z* -i  0, 0  <x  <  L,  L  is  the  lengdi  of  the 

deviy,tly.q)cciBcatiooofcuxTeotaloqgx»«(1^0<xr<L,artdt^n)edtotionof 
density  along  x*>L,0<x*<L.  Seif-consisterxyisitKhidedfatfaeanatyris. 

RESULTS 


Equlllbrima  xcro  bias  calculations:  The  first  set  of  calculations  represent 
cquilibriom  self-consistent  ealailatiotK  with  Maxwidlian  hmnvfaiyf  in 

- .  — *■  -  ^  J • ^ A  .  r^_  _  ■%  ^  A  _  ^ 


the  indicated  densities.  The  point  to  note  about  this  strocture  is 


*J“t  B  conasteirt  with  |2]  arad  fi 

the  intfOdncHon  of  the  gnanfitm  pfitential  in  fa]-  Itis 


of  the  basis  for 
tbatasfbekngdiof 
le  potential  increases 


bvxierb  centrally  placed  within  the  structure  is  shown  in  Figore  2.  Fortibu 
calculation  the  barrier  height  is  approximatety  lOmev  l«s  tnan  that  of  figure  1, 
Md  suggests  that  simple  estimates  of  the  bamer  in  mesoscopic  stroctures  are 


f  lira  C,  T*  n  riTTTTiT^O  ^  H  i%M-l  {,2  1 


flouble  warier  structure  of  (4]  who  obtain^  solotiotB  using  the  VWgner 
fommlatiotL  The  agremaent  for  this  zero  bias  calculation  b  rerrumable.  Note 
^  for  the  structme  with  doping  up  to  the  ^>acer  btyers,  there  is  a  r^on  local 

<h^e  accumulatioa  The  energy  of  the  carriers  is  near  and  about  the  resonant 


Flg.3.  Potcatu!  and  Density  vs.  Dotanoe. 


Fig.  4.  Potential  and  Density  vs.  Distance. 
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The  calculations  of  Figures  1  through  4,  assume  Maxwellian  boundary  conditions 
on  both  the  up-and  downstream  regions.  If  the  distribution  of  quantum  states  is 
different  for  me  two  boundaries,  it  is  expected  that  the  corresponding  distribution 
of  carriers  will  be  difference.  Several  cmculadons  of  this  type  were  performed 
under  the  assumption  of  Maxwellian  distributions  at  different  temperatures.  No 
physical  significance  is  attributed  to  this  temperature  difference;  rather  the 
calculations  are  displayed  to  shown  the  significance  of  the  boundary  conditions. 
Figure  S  shows  the  re^t  in  which  the  anode  boundary  temperature  is  higher  than 
that  of  the  u]»treain  contact  Fi^e  6  corresponds  to  the  case  when  the 
temperature  is  lower.  Note  that  in  one  case  there  is  local  accumulation  at  the 
downstream  region.  In  the  second  case  there  is  local  depletion. 

Finite  bias  calculations:  In  a  companion  paper,  15],  the  moments  of  the  Wimer 
function  are  used  to  examine  tran^rt  in  a  oouble  barrier  beterostnicture.  The 
barrier  height  is  21mev,  and  the  barrier  widths  are  a  narrow  S  A.  The  bias  across 
the  device^n  in  the  study  is  Sinev.  The  same  structure  in  the  structure  also 

examined  using  the  density  matrix.  The  equilSrrinm  distribution  of  carrier  is 
displayed  in  Figure  7  and  the  nonequilibrium  distribution  is  displwed  in  Hgure  8. 
Note  that  in  Figure  8,  current  int^on  with  an  entering  velocity  of  4x10*  cmfi  was 
assumed.  The  results  are  in  qualitative  agreement  with  that  of  [S]particularfy 
with  respect  to  the  almost  lin^  variation  of  potential  across  the  structure, 
although  we  note  the  bi^  carrier  density  at  tne  downstream  boundary.  As 
indicated  earlier  t^  value  is  determined  by  the  ^>ectra  of  energies  at  the 
downstream  contact 

Time  dependent  calculations:  A  non  self-consistent  time  dependent  problem  was 
solved  a  different  algorithm.  The  structure  was  a  650A  long  GaAs  element 
with  two  Si&thidc  barriers  separated  tty  a  50Ar^on.  At  time  t«»0^  a  uniform 
density  of  carriers  was  createa  While  this  initial  distribution  is  not  pltysical,  we 

E lint  out  that  within  the  first  IQfs  of  the  «dculation,  the  carriers  are  removed 
om  the  barriers  and  settle  within  the  well  and  outside  of  the  barriers. 
Subsequent  development  shows  that  as  the  carriers  outside  of  the  double 
barrier  move  toward  the  contacts  there  is  a  reduction  of  charge  within  the  well 
Upon  refiecdon  at  the  contacts  the  diatge  in  the  well  increases. 

CONCLUSIONS 

The  above  calculation  demonstrates  the  ability  of  the  density  matrix  to  examine 
fundamental  issues  associated  with  mesoscopic  structures.  Of  particular 
relevance  is  that  the  time  indqiendent  olo^tions  took  less  tl^  10  CPU  seconds 
and  thereby  provide  a  new  emdent  means  of  exploring  the  physics  of  quantum 
structures. 
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Abstract.  Barrier  calculations  based  upon  solutions  of  the  Liouville  equation  in 
the  coordinate  representation  reveal  a  complicated  spatial  dependence  of  the 
quantum  distribution  function  near  and  within  the  barriers.  Within  the  framework 
of  classical  transport  this  spatial  dependence  suggests  equilibrium  electron 
temperature  values  that  differ  from  the  ambient.  The  prospect  of  quantum  heating 
and  cooling  under  equilibrium  conditions  is  examined  and  dispelled  in  favour  of 
an  interpretation  that  includes  density-gradient  contributions. 


1.  Introduction 

Calculations  based  upon  solutions  of  the  Liouville 
equation  in  a  density  matrix  formulation  yield  a  com¬ 
plicated  spatial  distribution  with  a  mean  kinetic  energy 
in  equilibrium  that  differs  significantly  from  the  classical 
result  In  particular,  where  classical  physics  teaches  that 
the  energy  per  particle  is  T/l  per  degree  of  freedom  for 
a  Boltzmann  distribution,  quantum  physics,  as  pointed 
out  by  Wigner  [1],  teaches  otherwise.  The  origin  of  this 
difference  lies  in  the  presence  of  quantum  mechanical 
forces  arising  from  gradients  in  density  (Ancona  and 
lafrate,  [2]),  and  are  suggestive  of  a  spatially  dependent 
local  temperature  in  both  equilibrium  and  non¬ 
equilibrium  cases,  although  spatial  dependent  carrier 
temperature  in  equilibrium  introduces  interpretive  dif¬ 
ficulties.  To  avoid  this  difficulty  one  either  abandons  the 
spatial-dependent  temperature  description,  or  retains  it 
for  non-equilibrium  studies  and  seeks  another  descrip¬ 
tion  for  equilibrium.  But  in  either  case,  it  is  necessary  to 
demonstrate  its  origin.  This  is  provided  below  for 
equilibrium  conditions. 


2.  Energy  and  temperature 

Classical  physics  indicates  that  the  mean  kinetic  energy 
in  equilibrium  for  carriers  obeying  Boltzmann  statistics  is 

=  (2^  ^  (I) 

sps 

where  f(x,  p)  is  the  classical  distribution  function,  and  c, 
the  mean  kinetic  energy  per  particle,  is  independent  of 


position,  as  is  the  consequent  electron  temperature.  For 
quantum  structures.  In  which  quantum  distribution 
functions  are  required,  e  is  generally  spatially  dependent 
[1],  and  on  the  basis  of  equation  (1)  suggests  a  spatially 
dependent  equilibrium  carrier  temperature.  Because  of 
the  significance  of  carrier  temperature  in  interpreting 
hot-carrier  phenomena,  the  spatial  dependence  of  the 
mean  energy  per  particle,  and  the  origin  of  this  de¬ 
pendence  is  discussed  through  solutions  to  the  Liouville 
equation 

\hdp„Jdt  =  [H,  p,p]  (2) 

which  in  the  coordinate  representation  is  a  differential 
equation  for  p{x,  x',  t) 

dp/dt  -y  (ft/2miXV^  —  V'^)p 

-(l/ift)[K(x,t)- K(x',f)]p  =  0.  (3) 

To  expose  the  essential  features  of  this  discussion,  we 
assume  Boltzmann  statistics,  spatial  variations  only 
along  the  x  direction,  and  free  particle  behaviour  along 
the  y  and  z  directions.  Transforming  equation  (3)  to 
centre  of  mass  and  non-local  coordinates,  r=(x-t-x')/2, 
C=(x-x')/2,  p=f-p(r-f  C.r-O,  we  find 

p, + (hf2mi)p,^  ~  (l/ift)[  V(r  +  C,t)-  V{r  -  C,  t)]p  =  0.  (4) 

In  the  above  equation  subscripts  denote  differentiation. 
The  potential  V  in  equations  (3)  and  (4)  is  the  sum  of  all 
heterostructure  contributions,  ki(x),  and  contributions 
from  Poisson’s  equation: 

d/dxlt(x)dVJdx]  =  -e^lp{x,  x)-po(^)]  (5) 

where  the  subscript  ‘sc’  denotes  self-consistent;  Poix)  is 
the  background  ‘jellium’  doping  distribution.  The  dia¬ 
gonal  components  of  solutions  to  equation  (4)  (along  the 
diagonal  r  —  x  and  C  =  0)  provide  the  density,  while  the 
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expectation  value  of  energy  <£)  is  obtained  from  the 
diagonal  components  of  the  kinetic  energy  density  matrix 
[3] 

E{r  +  C,r  -  0  =  -h^/Smp...  (6) 

An  approximate  form  of  the  expectation  values  of  the 
density  and  the  energy  density  [1,4,5]  for  one  degree  of 
freedom  is: 

p{x)  =  p(x,  x)  =  Po  exp[ -(V  +  Q/2)/k^ 7]  (7) 

£(x)  =  [(kp  r/2)  +  (ft  V24mkB  T)V,Jp{x).  (8) 

In  equation  (7),  is  a  reference  density,  and  Q{x)  is  the 
Bohm  quantum  potential  (see,  e.g.,  [6]): 

Q(x)=-(h^l2mlp^i\Jp^'\  (9) 

The  second  term  of  equation  (8)  is  referred  to  as  the 
Wigner  contribution.  In  equilibrium  the  spatial  de¬ 
pendence  of  the  energy  per  particle,  e,  as  given  by 
equation  (8)  is  second  order  in  ft.  To  this  order,  if  the 
potential  appearing  in  equation  (8)  is  represented  by  the 
Boltzmann  relation  between  density  and  potential 
energy,  p(x)  =  p- er.p[  -  V(x)/kB7^,  it  is  seen  that  the 
spatial  dependence  of  e  is  a  direct  consequence  of  the 
spatial  derivatives  of  density.  In  this  context  the  origin  of 
the  quantum  correction  to  e  is  the  same  as  the  origin  of 
the  quantum  potential. 


3.  Calculations 

The  spatial  dependence  of  c  and  the  origin  of  the 
quantum  contributions  to  transport  arise  from  gradients 
in  the  carrier  density.  These  features  are  illustrated 
through  solutions  to  the  Liouville  equation  for  two 
equilibrium  solutions  using  Maxwellian  boundary  con¬ 
dition  as  discussed  in  [7].  Two  cases  are  considered.  For 
the  first  calculation  a  single  barrier  characterized  by  a 
potential 

Ko(x)  =  300(meV)exp[-(x/12.5A)"]  (10) 

is  placed  within  a  uniform,  1500  A  long  structure  doped 
to  10‘*cm"^.  The  two-dimensional  density  matrix, 
p(x,x')  as  obtained  from  the  Liouville  and  Poisson 
equations  is  displayed  in  figure  l(<i).  In  equilibrium  the 
density  matrix  is  real  and  symmetric,  p(x,x')  =  p(x',x), 
and  the  solution  is  completely  represented  by  one-half  of 
the  matrix  on  either  side  of  the  diagonal,  x  =  x',  as 
displayed  in  figure  1(a).  The  charge  density  p(x)  =  p(x,x) 
is  displayed  as  a  line  plot  in  figure  l(ft),  where  since  most 
of  the  structure  in  the  solution  is  contained  within  a 
range  of  250  A,  about  the  centre,  only  500  A  of  the  results 
are  displayed.  Figure  1(b)  displays  a  significant  reduction 
of  charge  within  the  barrier,  as  well  as  charge  accumula¬ 
tion  on  either  side  of  the  barrier.  While  the  reduction  of 
charge  within  the  barrier  is  a  consequence  of  the  presence 
of  the  barrier,  the  excess  charge  adjacent  to  the  barrier  is 
a  consequence  of  both  self-consistency  in  the  calculation 
and  wavcfunction  (or  density  matrix)  continuity  acro.ss 
the  h.irrier  The  spatial  deriendenre  of  the  charee  is 


consistent  with  the  condition  of  global  charge  neutrality. 
Figure  1(b)  also  displays  two  additional  plots.  The  curve 
reaching  the  lowest  value  of  density  within  the  barrier  is 
obtained  from  the  classical  Boltzmann  relation  betweeii 
density  and  potential  energy;  the  curve  reaching  inter¬ 
mediate  values  of  density  within  the  barrier  is  obtained 
from  equation  (7).  Note  that  away  from  the  barrier  the 
density  from  equation  (7)  approaches  a  value  that  is  less 
than  the  classical  value,  a  result  that  is  a  consequence  of  a 
change  in  curvature  of  the  potential  as  the  boundaries  are 
approached.  Neither  approximate  solution  can  be  regar¬ 
ded  as  an  adequate  representation  of  the  complete 
solution,  although  the  quantum-corrected  solution 
possesses  the  general  features  of  a  density  that  is  higher 
(than  classical)  within  the  barrier  and  lower  (than  clas¬ 
sical)  adjacent  to  the  barrier.  Figure  1(c)  displays  the 
potential  distribution.  The  lowering  of  the  potential 
adjacent  to  the  barrier  ( ~  20  meV)  is  a  consequence  of  the 
excess  charge  and  self-consistency.  Figure  !(</)  displays 
the  quantum  potential;  note  that  its  value  is  greater  than 
—  300  me V  in  the  centre  of  the  barrier.  The  energy 
density  matrix  represents  the  curvature  of  the  density 
matrix  in  the  non-local  direction.  As  seen  in  figure  1(a), 
the  curvature  is  steeper  where  there  is  excess  charge  and 
changes  sign  within  the  barrier.  The  mean  kinetic  energy 
per  particle,  e,  obtained  from  the  density  matrix  is 
displayed  in  figure  l(e\  along  with  the  Wigner  con¬ 
tribution  as  obtained  from  equation  (8).  It  is  apparent 
that  the  main  origin  of  the  structure  leading  to  the 
Wigner  contribution  i.'^  the  quantum  potential.  The 
negative  value  of  s  within  the  barrier  and  the  positive 
excess  value  of  energy  adjacent  to  the  barrier  suggest 
that  the  Wigner  contribution  is  not  a  correction,  but 
represents  a  dominant  effect,  and  that  temperature  con¬ 
cepts  (which  must  include  negative  values)  are  not  likely 
to  be  germane  within  the  context  of  equilibrium 
transport. 

The  spatial  dependence  of  the  mean  kinetic  energy 
per  particle  is  also  of  significance  in  multiple  barrier 
structures.  This  is  examined  for  a  double  barrier 
structure  with 

Fo(x)  =  (300meV){exp[-(x-75  A)/12.5  A]^ 

-t-exp[-(x-f  75  A)/12.5  A]^}  (11) 

The  barriers  are  centrally  placed  within  a  1500  A 
n*n”n*  structure  with  adjacent  10‘*cm“^n^  regions, 
and  a  centrally  placed  500  A,  10^’  cm  “  ^  region.  The  two- 
dimensional  density  matrix  is  displayed  in  figure  2(a), 
obtained  from  the  Liouville  and  Poisson  equations. 
There  is  excess  charge  between  the  barriers,  a  modest 
increase  in  curvature  between  the  barriers  and  a  change 
in  sign  of  the  curvature  within  the  barriers.  The  line  plot 
of  density  is  shown  in  figure  2(b)  over  a  reduced  range  of 
600  A.  The  density  as  obtained  from  equation  (7)  displays 
a  significantly  lower  charge  density  within  the  barrier  but 
order  of  magnitude  agreement  within  the  quantum  well. 
The  classical  solution  for  density  is  completely  unac¬ 
ceptable.  The  potential  distribution,  shown  in  figure  2(c), 
reaches  flat-band  bevond  4(X)  A  on  either  side  of  the 
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Figure  1.  (a)  Density  matrix  for  a  single-barrier  structure.  The  physical  dimension  of  the  structure  is  1500  A,  requiring  that 
the  density  matrix,  which  is  calculated  over  a  square  matrix,  is  of  side  1500A/y2.  The  centre  of  mass  and  non-local 
directions  are  indicated;  (b)  diagonal  component  of  the  density  matrix  ( - ),  from  equation  (7)  ( - ),  classical  relation 

(  )'■  (c)  potential  energy  V'{x):  {d)  quantum  potential;  (e)  energy  per  particle  from  density  matrix  (-  -  -),  from 

equation  (8)  ( - ). 


origin;  its  increase  arises  from  self-consistency  and  the 
reduction  of  charge  in  the  low-doped  region  compared 
with  the  bounding  charge  density.  The  quantum  potent¬ 
ial  displayed  in  figure  2(d)  is  positive  within  the  quantum 
well,  and  emphasizes  the  reduction  in  charge  density 
compared  with  the  classical  value;  it  is  negative  within 
the  barriers,  as  in  the  case  of  the  single-barrier  structure, 
and  positive  outside  of  the  barriers.  The  jxrsitive  value 
outside  of  the  barriers  is  a  consequence  of  wavefunction 
and  density  matrix  continuity  within  the  classically 
accessible  region.  Note  again  that  the  structure  of  the 
quantum  potential  is  apparently  the  main  origin  of  the 
structure  leading  to  the  Wigner  contribution  to  the 


energy  per  particle,  as  seen  in  figure  2(e).  As  in  the  case  of 
the  single  barrier  of  figure  1  the  calculations  suggest  that 
the  Wigner  contribution  is  not  a  correction  but  repres¬ 
ents  a  dominant  effect,  and  that  tempierature  concepts  are 
not  germane  in  the  context  of  equilibrium  transport. 

4.  Conclusions 

The  calculations  of  figures  1  and  2  reveal  significant 
spatial  variations  in  energy  associated  with  density  gra¬ 
dients.  Mathematically,  these  energy  variations,  which 
arc  a  consequence  of  wavefunction  continuity  as  repre- 
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Figure  2.  (a)  Density  matrix  for  a  double-barrier  structure;  (/»)  diagonal  component  of  the  density  matrix  ( — ),  from 

equation  (7)  ( - ),  classical  relation  ( - );  (c)  potential  energy  V'(x);  {cf)  quantum  potential;  (e)  energy  per  particle 

from  density  matrix  ( — ),  from  equation  (8)  ( - ). 


sented  by  ‘curvature’  in  the  density  matrix,  are  suggestive 
of  quantum  heating  and/or  cooling.  Physically  these 
energy  variations  represent  the  influence  of  local 
quantum  mechanical  density  dependent  forces  on  the 
carriers.  While  their  magnitudes  indicate  that  they  must 
be  accounted  for  in  all  quantum  mechanical  treatments 
of  transport  in  mesoscopic  structures,  an  interpretation 
in  terms  of  heating  or  cooling  in  equilibrium  is 
problematic. 
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Abstract.  The  time-dependent  Uouville  equation  for  the  density  matrix  in  the 
coordinate  representation,  incorporating  scattering  ejects  through  a  quasi- Fermi 
level,  and  Poisson's  equation  have  been  solved  numerically  for  a  single-  and 
double-barrier  structure  using  algorithms  based  on  'characteristics',  and  show 
significant  charge  accumulation  on  the  emitter  side  of  the  barrier,  as  well  as 
significant  charge  screening  in  multiple-barrier  structures. 


1.  Introduction 

Barrier  calculations  based  upon  solutions  to  the 
Liouville  equation  with  a  density  matrix  formulation 
display  signibcaut  charge  screening.  In  particular,  the 
Liouville  equation  which  in  the  coordinate  represen¬ 
tation  and  in  the  absence  of  scattering  is  the  differential 
equation 

dp/dt  +  (/i/2mi)(V^-V'^)p-(l/ifi)[F(jt,t)-  K(x',0]p  =  0 

(i; 

yields  solutions  for  p{x,x',t)  that  contain  significant 
quantum  departures  from  the  classical  solution  when 
there  are  strong  gradients  in  the  carrier  density.  We 
illustrate  these  features  for  Boltzmann  statistics,  spatial 
variations  only  along  the  x  direction,  and  free  particle 
behaviour  along  the  y  and  z  directions.  We  also  trans¬ 
form  equation  (1)  to  centre  of  mass  and  non-local 
coordinates,  r  =  (x  +  x")/!, 

C  =  (x  -  x')/2,p=>  p{r  +  C,r-n 

p,  +  (h/2mi)p^  -  (lM)[K(r  -F  C,f)  -  K(r  -  C,t)]p  =  0. 

(2) 

In  the  above  equation  subscripts  denote  differentiation, 
and  the  potential  V  includes  all  heterostructure  con¬ 
tributions,  lo(x),  as  well  as  contributions  from  Poisson’s 
equation: 

(d/dx)[£(x)dKc/^x]  =  -e^[p(x,x)  -  Po(x)]  (3) 

The  subscript  ‘sc’  denotes  self-consistent,  and  Po(x)  is  the 
background  ‘jellium’  doping  distribution  Note:  along  the 
diagonal  r  =  x  and  C  =  0- 

Significant  quantum  effects  which  are  revealed 
through  numerical  solutions  are  also  revealed  through 
approximate  solutions  when  the  potential  in  equation  (2) 
is  expanded  in  a  Taylor  scries  and  only  the  first  two  non¬ 


zero  terms  are  retained 

K(r  -F  0  -  K(r  -  0  «  2C f;  -F  C'  K,„/3  (4) 

To  second  order  in  ft,  a  solution  to  equation  (2)  is  [Ij 

p(r  -F  C,  r  -  0  =  Poexp[-/?(K+  Q/3) 

-F(CM)H1  +  ^A^K„/6)]  (5) 

where  p  =  l/k^T,  =  ph^/2m  and 

Q  =  (AV2XK„  -  (6) 

With  Q  regarded  as  a  quantum  correction,  inserting  the 
Boltzmann  relation  between  density  and  potential 
iPV  =  ~\n p),  Q  is  transformed  to 

C(x)=  -(ftV2mXp’'\>’'^  (7) 

which  is  a  generalization  of  the  Bohm  quantum  potential 
(e.g.,  see  [2])  for  a  multiparticle  system  and  teaches  that 
screening  is  significant  when  there  are  strong  gradients  in 
carrier  density. 

While  the  above  discussion  ignores  scattering,  when 
relaxation  effects  similar  to  Fokker- Planck  dissipation 
are  introduced  the  diagonal  component  of  density  is 
approximately  given  by 

p{x)  =  p{x,  x)  =  PoCxpC  -  P(  K(x)  -  (^(x)  -F  Q(x)/3)]  (8) 

with  a  velocity  flux  density  proportional  to  —pd^/dx, 
and  suggests  that  there  arc  situations  in  which  the  effect 
of  scattering  on  the  density  can  be  qualitatively  represen¬ 
ted  by  introducing  a  scattering  operator,  part  of  which 
consists  of  an  algebraic  contribution,  (^(x).  As  a  rationale, 
consider  an  n^n  n*  structure  in  which  the  heavily 
doped  regions  are  long  enough  for  the  carriers  to  relax  to 
background.  Under  bias  assume  a  quasi-Fermi  level 
exists  whose  value  is  such  that  the  density  is  constant  at 
the  upstream  and  downstream  boundaries  This  as¬ 
sumption  IS  implicit  in  most  device  analyses  (e  p  .  see  [  Ul- 
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with  the  conbincd  effects  of  an  applied  bias  and  scattering 
elsewhere  in  the  structure  determining  the  remaining 
charge  distribution. 

In  the  absence  of  a  detailed  description  of  scattering 
in  structures,  how  may  one  expect  the  density  to  vary? 
This  is  explored  through  the  assumption  that  an  alge¬ 
braic  scattering  potential,  <f(x),  qualitatively  represents 
the  effects  of  scattering  on  the  density  within  the 
structure,  c^nd  approaches  the  quasi-Fermi  level,  E,, 
within  the  boundary  regions,  ^(x)  is  taken  as  constant 
and  equal  to  the  value  of  £,  at  the  emitter  boundary  £f(c), 
until  the  quantum  structure  is  reached.  At  the  end  of  the 
quantum  structure  (^(x)  is  set  equal  to  the  value  of  E,  at 
the  collector  boundary,  Ef{c).  The  intermediate  values 
within  the  quantum  region  are  indicated  below.  In  terms 
of  equation  (2),  incorporating  if(x)  amounts  to  replacing 
T(x)  with  K(x)  —  ^(x). 


2.  Calculations 

The  Liouvillo  equation  was  solved  for  1500  A  n^nn^ 
structures  under  bias.  Displaced  Maxwellian  conditions 
are  assuuied  and  a  flat  band  occurs  at  the  boundaries.  In 
one  case  a  single  barrier  50  A  wide  and  300  meV  high  was 
centrally  located  (figure  1  calculations);  in  the  second  case 
two  barriers  each  50  A  wide  and  300  meV  high  were 
centrally  placed  and  separated  y  50  A  (figure  2  cal¬ 
culations).  For  the  single-barrier  calculations  i{x)  was 
chosen  to  ensure  that  the  density  at  the  emitter  and 
collector  regions  was  equal  to  the  background  of 
10‘®cm‘^.  The  values  of  the  algebraic  scattering  potent¬ 
ial  across  the  barrier  varied  linearly  from  the  value  E^e) 
to  Ef(c).  The  small  field  that  forms  in  the  heavily  doped 
regions  under  bias  was  ignored  and  the  current  through 
the  device  was  not  accurately  represented.  Higher  current 
values  that  would  normally  flow  were  computed  from  a 
simple  circuit  equation  and  yielded  density  distributions 
similar  to  that  obtained  when  the  current  contribution 
was  ignored.  The  real  part  of  the  two-dimensional 
density  matrix  pfx,  x')  as  obtained  from  the  Liouville  and 
Poisson  equations  is  displayed  for  a  bias  range 
0  <  Kppiicd  <  -  300meV.  p(x,x^  is  symmetric  and 
completely  represented  by  one  half  of  the  matrix  on 
either  side  of  the  diagonal,  x  =  x'.  The  density  matrix 
shows  a  build-up  of  charge  on  the  emitter  side  of  the 
barrier  as  well  a  3  a  broad  depletion  region  on  the 
collector  side.  Linear  plots  of  the  density  and  potential 
are  displayed  in  figure  1(e)  and  !(/)  respectively.  As  the 
bias  increases  there  is  a  lowering  of  the  barrier,  and  a 
build-up  of  a  ‘notch’  p>otcntial  at  the  emitter  side  of  the 
barrier,  signifying  the  development  of  a  region  of  charge 
accumulation;  the  collector  side  of  the  barrier  shows 
nearly  linear  variation  in  potential  suggestive  of  a  broad 
region  of  charge  depletion.  Both  accumulation  and  de¬ 
piction  regions  for  single-barrier  structures  have  been 
discussed  by  Laves  ei  al  ^4].  Increased  charge  accumula¬ 
tion  on  the  emitter  side  of  the  barrier  tends  to  reduce  the 
relative  change  in  potential  on  the  emitter  side  of  the 
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Note  the  calculations  indicate  that  global  charge  neu¬ 
trality  occurs  (including  the  double-bairier  calculations); 
i.e.  {dx(p  -  Pj)  =  0.  It  is  important  to  note  that  tiie 
quantitative  value  of  the  charge  adjacent  to  the  barrier  is 
dependent  upon  the  value  of  (^(x).  For  example  if  (^(x) 
were  suddenly  reduced  prior  to  the  barrier,  the  density 
prior  to  the  barrier  would  also  be  reduced.  Global  charge 
conservation  would  require  that  local  accumulation 
regions  would  form  elsewhere. 

For  the  double-barrier  calculations  the  low  doped 
region  extends  over  a  distance  of  600  A.  Two  sets  of 
simulations  were  performed.  In  the  first  case  the  bias  on 
the  collector  ranged  from  0  <  ^  —  300meV. 

Here,  the  intermediate  value  of  (^(x)  was  equal  to 
E^e)  to  the  beginning  of  the  second  barrier  and  varied 
'inearly  to  £f(c)  at  the  end  of  the  second  barrier.  For  the 
second  case  the  bias  on  the  collector  ranged  from 

-  lOOmeV  <  ^  -  500  me V  and  the  linear  varia¬ 

tion  of  s(x)  started  at  the  beginning  of  the  emitter  barrier 
and  ended  at  the  end  of  the  collector  barrier.  Figures  2(a- 
d)  represent  the  low-bias  range,  while  figures  2(e~h'j  arc 
for  the  high-bias  range.  Figures  2(a,  b  and  c)  display 
density,  potential  and  Q  +  V  fwithin  select  regions  of  the 
structure).  Figure  2(d)  displays  the  real  part  of  the  density 
matrix  for  an  applied  voltage  of  —  300  me V.  The  potent¬ 
ial  variation  shows  the  presence  of  the  formation  of  a 
notch  region  at  the  emitter  side  of  the  first  barrier,  a 
smaller  fraction  of  potential  falling  across  the  emitter 
barrier,  compared  to  the  collector  barrier,  as  the  bias  is 
increased.  Indeed,  further  increases  in  bias  for  this  varia¬ 
tion  of  ^(x)  result  in  small  variations  in  potential  drop 
across  the  emitter  barrier  region  [5].  The  density  distri¬ 
bution  shows  the  build  up  of  charge  in  the  quantum  well, 
and  an  increased  accumulation  on  the  emitter  side  cf  the 
barrier.  At  higher  bias  levels  the  charge  in  the  well  begins 
to  screen  the  emitter  barrier  from  the  collector  and  the 
potential  drop  across  the  emitter  region  is  smaller  than 
the  corresponding  case  of  the  single  barrier.  The  term 
Q  ■¥  V  within  the  quantum  well  is  displayed  because  its 
value  is  approximately  equal  to  the  value  of  the  quasi¬ 
bound  state  energy  and  moves  down  as  the  bias  is 
increased.  Successive  increases  in  bias,  reduce  the  move¬ 
ment  of  the  emitter  barrier  and  delay  or  inhibit  the 
movement  of  the  quasi-bound  state  to  values  below  the 
£f(e),  particularly  when  the  n  “  regions  are  very  narrow. 

The  reduction  of  charge  in  the  well  occurs  with  the 
variation  in  ij(x)  as  indicated  for  the  second  range  of  bias 
levels  in  figures  2(e-h).  (Figure  2(h)  is  the  real  part  of  the 
density  matrix  for  a  value  of  applied  potential  equal  to 

—  300  me V  and  is  directly  compared  to  figure  (</)).  The 
result  is  a  distribution  of  charge  and  potential  in  the 
vicinity  of  the  emitter  barrier  more  closely  related  to  that 
across  the  single-barrier  structure.  There  is  a  more 
equitable  potential  drop  across  the  double  barriers  [5], 
and  a  skewing  of  the  Q  +  V  terms.  A  comparison  of  the 
potential  distribution  and  charge  distribution  at 
Kppiitd  =  —300  meV  indicates  that  a  larger  potential 
drop  falls  across  the  double  barrier  in  the  absence  of 
charge  in  the  well,  which  as  a  consequence  will  lead  to 
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Figure  1 .  (a-d)  Density  matrix  for  a  single-barrier  structure.  The  physical  dimension  of  the 
structure  is  1500  A.  requiring  that  the  density  matrix,  which  is  calculated  over  square  matrix, 
is  of  side  1500A/.^2.  The  centre  of  mass  and  non-local  directions  are  indicated,  (a) 

»^„o...d  =  0.0meV,  (6)  - lOOmeV.  (c)  =  - 300 meV.  (d) 

^•ppii»d=  —  BOOmeV.  (e)  Diagonal  component  of  density  matrix:  ( - )  0.0 meV;  ( - ) 

-lOOmeV;  ( - )  -300meV;  ( - )  — 500meV.  {f)  Potential  energy  ( - )  0.0  meV; 

( - )  -lOOmeV;  ( - )  -300meV;  ( - )  -500meV. 


in  the  measured  current  voltage  characteristics  (see,  e.g., 
[4])  is  a  measure  of  the  variations  in  the  charge  distribu¬ 
tion  as  a  function  of  scattering  within  the  structure. 
Indeed  if  we  assume  that  the  range  of  scattering  potential 
where  <J(x)  is  constant  to  the  collector  barrier  is  represen¬ 
tative  of  the  scattering  dynamics  until  the  peak  current  is 
reached  in  resonant  tunnelling  structures,  and  that  scat¬ 


tering  results  in  a  reduction  of  ^(x)  prior  to  the  collector 
barrier,  for  the  next  and  higher  range  of  applied  bias 
levels,  then  the  build  up  of  charge  in  the  well  of  the 
resonant  tunnelling  diode,  and  the  subsequent  redistri¬ 
bution  of  charge  to  a  region  upstream  of  the  emitter 
barrier  at  higher  bias  levels  emerges  from  this 
de.scription. 


0(Mir.ity  iii.iinx  .iiu)  ()ii.intum  moment  studies 


Figure  2.  Double-barrier  structure  under  bias;  (a  and  d)  density,  (b  and  f)  potential  energy 
(c  and  g)  Q+V  aX  select  portions  of  the  structure,  (dand  h)  density  matrix  at  —  300meV. 
Differences  in  result  depend  upon  algebraic  scattering  potential  (see  text),  (a),  (6)  and  (c) 

( — )  O.Ome'/;  ( - )  — lOOmeV:  ( - )  — SOOmeV.  (e),  {f)  and  (g)  ( — )  O.OmeV;  ( - ) 

-300meV:  ( - )  —  500meV, 


3.  Summary 

An  appropriate  description  of  transport  requires,  at  least 
at  the  boundaries,  that  the  carriers  approach  their 
anticipated  backpround  values  This  is  achieved  Ihrouph 
the  introduction  of  scattermp  events.  I  he  presence  of 
scattering  events  alters  tin*  I  loiiville  couation.  with  a 


consequent  change  in  the  distribution  of  density  and  in 
current.  In  the  absence  of  a  detailed  description  of 
scattering  its  has  been  introduced  through  the 

incorporation  of  an  algebraic  scattering  potential.  While 
the  results  presented  here  are  dejaendent  upon  the  specific 
representation  of  c(\)  they  are  representative  of  the  fact 
the  fact  tli.it  a  simolc  nicliirc  of  Iransnorl  as  nrovided  hv 
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standard  kinetic  and  potential  energy  concepts,  albeit 
quantum  mechanical,  is  not  likely  to  be  adequate  to 
describe  transport  in  multiple-barrier  systems. 
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1.1  INTRODUCTION 

It  is  arguable  that  the  history  of  gallium  arsenide  semiconductor  devices, 
from  the  early  1960s  to  the  present  time,  falls  into  three  groups.  First,  there 
was  the  experimental  work  of  Gunn  [1],  demonstrating  the  generation  of 
sustained  oscillations  upon  application  of  a  sufficiently  large  dc  bias.  This 
work  opened  up  the  possibility  of  fabricating  bulk  microwave  and  mil¬ 
limeter-wave  devices,  and  hastened  additional  and  intense  studies  of  the 
properties  of  compound  semiconductor  devices.  Second,  there  was  the  study 
of  Ruch  [2],  whose  results  suggested  that  the  transient,  or  nonsteady-state, 
aspects  of  semiconductor  transport  would  improve  the  speed  of  devices  by 
almost  an  order  of  magnitude.  This,  of  course,  is  the  argument  behind  much 
of  the  move  toward  submicron  and  ultrasubmicron  structures.  The  third  era, 
the  one  we  are  presently  in,  involves  the  incorporation  of  gallium  arsenide 
into  material-engineered  highly  complex  structures,  some  of  which  have 
provided  remarkable  millimeter  wave  characteristics,  such  as  the  pseudo- 
morphic  HEMT  [3].  Much  of  this  book  is  concerned  with  this  third  era,  and 
thus  this  chapter  will  only  briefly  touch  upon  it.  Rather,  this  section  will 
present  a  road  map  of  the  consequences  of  using  compound  semicon¬ 
ductors  for  device  applications,  using  gallium  arsenide  as  the  paradigm 
example. 

The  band  structure  of  gallium  arsenide  is  familiar  to  most  and  is  displayed 
in  Fig.  1.  [4].  It  is  a  direct  bandgap  material.  Hie  minimum  in  the 
conduction  band  is  at  F  with  relevant  subsidiary  conduction  band  minima  at 
L  and  X.  The  curvature  at  F  is  such  that  the  effective  mass  of  the  F-valley  is 
lower  than  that  of  the  next  two  adjacent  subsidiary  L-  and  A^valleys.  For 
the  valence  band,  the  two  valleys  of  significance  are  those  associated  with 
the  light  and  heavy  holes.  We  will  concentrate  on  transport  contributions 
from  these  five  valleys. 

In  equilibrium,  the  relative  population  of  electrons  in  the  valleys  is 
dependent  on  the  density  of  available  states  and  the  energy  separation,  for 
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Figure  1.1  Band  structure  of  the  semiconductor  gallium  arsenide  [4],  (Reprinted, 
with  permission,  from  Physical  Review  [Section^  B:  Solid  State.) 


example: 


where  and  denote  the  equilibrium  density  of  the  F-  and  L-valley 
carriers,  respectively,  and  are  their  effective  masses,  and  A  is  the 
F-L  energy  separation.  Thus,  in  equilibrium  virtually  all  of  the  electrons  of 
interest  are  in  the  F-valley.  For  the  holes,  the  valleys  are  degenerate. 

Gallium  arsenide  is  a  compound  semiconductor.  At  low  values  of  electric 
field,  apart  from  carrier-carrier  scattering,  there  are  three  important  scatter¬ 
ing  mechanisms:  polar  optical  phonon  scattering,  acoustic  phonon  scatter¬ 
ing,  and  impurity  scattering.  For  F-valley  electrons,  the  contribution  to  the 
momentum  scattering  rate  from  polar  optical  phonons  is  approximately  two 
orders  of  magnitude  larger  than  that  of  the  acoustic  phonon.  Since,  with 
regard  to  mobility,  scattering  rates  are  additive,  the  polar  optical  phonon  is 
the  dominant  scatterer.  Ideal  room-temperature  electron  mobilities  are  in 
the  range  of  8000  to  9000cm^/V*s.  For  the  subsidiary  valleys,  the  effective 
masses  of  the  carriers  are  much  larger  than  that  of  the  F-valley,  and  the 
relative  contribution  of  the  acoustic  phonon  increases.  Nevertheless,  the 
polar  phonon  dominates  the  transport.  For  holes,  the  situation  is  mixed. 
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with  the  dominant  scattering  being  polar  and  nonpolar  deformation  poten¬ 
tial  coupling.  For  momentum  scattering,  the  nonpolar  deformation  potential 
scattering  dominates. 

At  high  values  of  electric  field  and  for  electrons,  nonpolar  phonons  enter 
the  picture,  intervalley  transfer  from  F  to  L  takes  place,  and  the  situation 
becomes  complex.  For  example,  the  spatially  uniform,  field-dependent 
velocity  characteristics  of  gallium  arsenide,  ignoring  electron-hole  interac¬ 
tion,  display  a  region  of  negative  differential  mobility,  as  shown  in  Fig.  1.2 
[5],  where  at  values  of  field  in  excess  of  3  kV/cm  the  mean  carrier  velocity 
begins  to  decrease  with  increasing  electric  field.  This  is  an  unusual  situation 
and  it  is  perhaps  important  to  recognize  that  the  mean  electron  velocity  of  a 
given  species  of  carrier,  assuming  a  parabolic  band,  is  not  decreasing  with 
increasing  electric  field.  Rather,  the  numbers  of  high-mobility  electrons  are 
decreasing,  due  to  transfer  to  the  subsidiary  larger  effective  mass  valleys. 

The  situation  with  holes  is  different.  Here,  the  dominant  transport  is 
through  the  heavy  hole.  Interband  hole  scattering  is  always  present  even  at 
very  low  fields,  however  the  relative  population  is  fixed  through  the  ratio  of 
the  effective  masses,  and  the  existence  of  a  dc  negative  conductance  for 
holes,  on  the  basis  of  available  data,  is  ruled  out.  The  field  dependence  of 
the  mean  hole  velocity,  ignoring  interaction  with  the  electrons,  is  displayed 
in  Fig.  1.3  [6],  and  there  are  two  important  features  of  note.  First,  there  is 
the  extremely  low  mobility  of  the  holes  at  low-field  values.  Second,  there  is 
the  saturated  drift  velocity,  which  is  expected  to  be  higher  than  that  of 
electrons  at  high  fields.  We  note  there  is  no  hard  data  on  the  high-field 
carrier  velocity  of  holes  in  gallium  arsenide. 

Calculations  of  the  type  displayed  in  Figs.  1.2  and  1.3  have  been 
described  by  many  workers  and  are  routinely  incorporated  into  simulation 
codes.  Of  more  recent  interest,  because  of  mixed  conduction  heterostruc¬ 
ture  devices,  are  the  modifications  that  may  be  expected  when  electron-hole 


FIELD  (kv/cm) 

Figure  1.2  Field-dependent  electron  mean  velocity  for  gallium  arsenide  [5). 


INTRODUCTION  TO  THE  PHYSICS  OF  GALLIUM  ARSENIDE  DEVICES 


E  (V/cm) 

Figure  1.3  Field-dependent  hole  mean  velocity  for  gallium  arsenide  [6].  (Re¬ 
printed,  with  permission,  from  Journal  of  Applied  Physics.) 


Figure  1.4  Field-dependent  electron  mean  velocity  for  gallium  indium  arsenide 
assuming  an  interaction  with  heavy  holes  [8].  (Reprinted,  with  permission,  from 
Applied  Physics  Letters.) 


THE  ROLE  OF  BAND  STRUCTURE  ON  DEVICE  OPERATION  5 


scattering  occurs.  However,  because  of  the  limited  number  of  studies  with 
GaAs  and  because  of  similarities  with  other  compound  semiconductors, 
results  of  InGaAs  studies  are  presented.  Additionally,  because  of  ex¬ 
perimental  work  in  InGaAs  [7],  the  role  of  carrier-carrier  scattering  has 
been  most  extensively  studied  for  that  material.  Monte  Carlo  calculations 
incorporating  electron-hole  scattering  are  displayed  in  Fig.  1.4  [8].  The 
results  require  some  detailed  discussion  and  are  considered  later.  Here,  we 
simply  note  that  the  presence  of  holes  leads  to  reduction  in  the  low-field 
mobility  but  an  increased  peak  carrier  velocity.  These  intriguing  results  are 
also  anticipated  for  GaAs. 

The  question  of  interest  is  how  may  we  expect  the  role  of  the  complicated 
compound  semiconductor  band  structure  to  affect  the  performance  of 
devices.  This  is  considered  next. 


1.2  THE  ROLE  OF  BAND  STRUCTURE  ON  THE  OPERATION 
OF  ELECTRON  DEVICES 

In  examining  the  role  of  band  structure  on  the  operation  of  electron  devices, 
there  are  several  items  of  immediate  interest:  the  effective  mass,  the 
low-field  mobility,  and  the  direct  bandgap  energy  of  the  binary  III-V 
materials  (see  Table  1.1  [9]).  Additionally,  the  energy  separation  of  the 
conduction  band  minima  to  the  subsidiary  valleys  is  listed  in  Table  1.2  [10]. 
Note  that  of  the  seven  binary  materials  listed,  five  are  direct  bandgap 
materials,  and  two,  GaP  and  AlAs,  are  indirect  materials.  The  indirect 


TABLE  1.1  Critical  Parameters  of  Select  Compound 
Semiconductors 


Compound 

Effective  Mass* 

Electron 

Low-field 

Mobility 

(cm^/V-s) 

Direct  Energy 
Bandgap 
(eV) 

GaAs 

0.063" 

9,200 

1.424 

GaP' 

0.25V0.91' 

160 

2.78 

GaSb 

0.042 

3,750 

0.75 

InAs 

0.0219 

33,000 

0.354 

InP 

0.079 

5,370 

1.344 

InSb 

0.0136 

77,000 

0.230 

AlAs** 

0.71" 

300 

2.98 

Source:  Ref.  9.  Reprinted  with  permission  of  Springer-Verlag. 
'Multiples  of  free  electron  mass  at  the  conduction  band  minima. 
"DOS. 

'The  minima  in  the  conduction  ba;id  are  at  A-axis  near  zone  boundary. 
‘‘The  minimum  in  the  conduction  band  is  at  X. 
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TABLE  1.2  Intervalley  Energy  Separation 


Compound 

r-L 

(eV) 

r-A' 

(eV) 

L-X 

(eV) 

GaAs 

0.34 

0.48 

0.14 

Gap 

-0.27 

-0.39 

-0.37 

GaSb 

0.08 

0.37 

0.23 

InAs 

1.27 

1.60 

0.33 

InP 

0-63 

0.73 

0.10 

InSb 

0.41 

0.97 

0.56 

AlAs 

-0.15 

-0.79 

-0.64 

Source:  Ref.  10.  Reprinted  with  permission  from  Oxford 
University  Press. 


bandgap  materials  have  the  highest  effective  masses  of  the  group  and  also 
the  lowest  mobility.  Of  these  materials,  GaAs,  InAs,  InP,  and  InSb  possess 
regions  of  negative  differential  mobility.  GaSb,  GaP,  and  ALAs  do  not.  It  is 
perhaps  not  surprising  that  the  first  four  mentioned  materials  possess  a 
region  of  negative  differential  mobility,  nor  that  the  last  two  materials  do 
not.  In  the  latter  case,  the  minima  in  energy  is  associated  with  a  large 
effective  mass,  high  density-of-states  energy  level.  The  situation  with  GaSb 
is  peculiar.  But  here,  while  the  effective  mass  of  the  F-valley  is  the  smallest 
of  the  three,  its  closeness  in  energy  to  that  of  the  subsidiary  L- valley  is  such 
that  at  low  values  of  field  conduction,  contributions  arise  from  both  the  T- 
and  L- valley,  effectively  suppressing  the  contributions  of  intervalley  transfer 
to  negative  differential  conductivity. 

The  presence  of  a  region  of  bulk  negative  differential  mobility  has,  as  a 
major  consequence,  the  possibility  of  electrical  instabilities.  These  in¬ 
stabilities  manifest  themselves  either  as  large-signal  dipole-dominated  oscil¬ 
lations,  often  referred  to  as  the  Gunn  effect,  or  as  circuit-controlled 
oscillations,  where  the  semiconductor  behaves  electrically  as  a  van  der  Pol 
oscillator.  The  binary  semiconductors  GaAs,  InP,  and  InAs  have  exhibited 
electrical  instabilities  associated  with  bulk  negative  differential  mobility. 
While  InSb  has  also  sustained  electrical  instabilities,  the  interpretation  of 
the  instability  is  complicated  by  the  small  direct  bandgap  and  the  possibility 
of  avalanching  at  low  bias  levels. 

An  additional  feature  of  importance  is  the  intrinsic  carrier  concentrations 
of  some  of  these  materials,  as  shown  in  Table  1.3  [9].  It  is  clear  that  the 
intrinsic  concentration  of  InAs  and  InSb  make  them  unsuitable  for  a 
unipolar  source.  Indeed,  all  transport  calculations  using  these  latter  materi¬ 
als  must  necessarily  include  multispecies  transport. 

In  choosing  materials  for  electron  devices,  particularly  as  power  sources, 
a  figure  of  merit  has  been  the  peak  to  saturated  drift  velocity  ratio.  From 
this  point  of  view,  indium  phosphide  is  an  attractive  candidate,  but  this  must 
be  weighed  with  the  fact  that  the  low-field  mobility  of  InP  is  less  than  that  of 
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TABLE  1.3  Intrinsic  Concentration 


Compound 

n  (/cm^) 

GaAs 

2.1  X  10* 

InAs 

1.3  X  10'* 

InSb 

2.0  X  10** 

InP 

1.2  X  10® 

Source:  Ref.  9.  Reprinted  with  permission  of  Springer- 
Verlag. 


gallium  arsenide.  A  recent  study  comparing  these  features  suggests  that  the 
F-valley  mobility  is  the  dominant  material  parameter  of  submicron  struc¬ 
tures,  whereas  the  high-field  saturated  drift  velocity  is  the  dominant  material 
parameter  of  micron-length  structures  [11].^  Additionally,  if  a  choice  for  two 
terminal  sources  is  to  be  made  between,  for  example,  InP  and  GaAs,  other 
issues  emerge.  For  instance,  the  scattering  rates  in  InP  indicate  a  shorter 
energy  relaxation  time  than  that  of  GaAs.  The  consequence  of  this  are 
higher-frequency  operation  for  InP.  Thus,  at  least  with  respect  to  these 
materials,  the  peak  to  valley  ratio  of  the  materials  is  only  one  factor  in  the 
design  of  an  electrical  source. 


V  (drain-source)  (V) 

Figure  1.5  Temperature-dependent  pulsed  data  for  a  GaAs  FET,  with  a  3.0  /im 
gate  length,  a  source  to  drain  separation  of  8.5  ptm,  and  an  epitaxial  thickness  of 
3000  ±  500  A.  Nominal  background  doping  is  10‘’/cm’  [12]. 


^(Editor’s  Note)  See,  however.  Chap.  2,  page  82,  last  paragraph  and  reference  [12] 
cited  within  for  an  additional  view  on  this  issue. 
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There  is  less  to  say  about  the  effects  of  negative  differential  mobility  on 
the  operation  of  avalanche  diodes.  Here,  the  effects  of  negative  differential 
mobility  conductivity  are  present  but  are  overshadowed  by  the  effects  of 
avalanching.  For  example,  recent  simulation  studies  show  the  presence  of 
domains  in  IMPATTs,  whose  presence  is  a  direct  consequence  of  negative 
differential  conductivity.  These  domains  can  complicate  the  actual  transit 
time  of  dipole  layers  associated  with  the  avalanche  generation,  but  the 
negative  differential  mobility  is  a  marginal  issue.  Such  is  not  the  case  with 
three-terminal  devices. 

For  three-terminal  device  observations  of  bias-dependent  white  light  in 
GaAs  FETs,  as  from  either  the  drain  side  of  the  gate  contact  and  the  gate 
side  of  the  drain  contact,  are  consistent  with  numerical  calculations  showing 
the  presence  of  local  high-field  dipole  layers  near  the  gate  and  drain 
contacts.  In  addition,  for  a  range  of  bias,  some  devices  display  a  current 
dropback  consistent  with  bias-dependent  formation  of  high-field  domains 
and  concurrent  current  oscillations.  This  last  result  is  shown  in  Fig.  1.5  [12]. 
Remaining  questions  of  interest  focus  on  the  manner  in  which  transport  in 
these  devices  is  examined.  We  begin  with  the  equilibrium  description  of 
transport. 


1.3  EQUILIBRIUM  DESCRIPTION  OF  TRANSPORT 

The  steady-state  equilibrium  description  of  transport  has  traditionally  pro¬ 
vided  most  details  of  device  behavior.  Nevertheless,  the  description  ignores 
acceleration.  It  assumes  that  the  carrier  velocity  is  determined  by  the  local 
electric  field  and  that  the  total  current  is  governed  by  a  balance  of  a  drift 
component  and  a  diffusive  component.  Typically,  the  continuity  equation  is 
solved  simultaneously  with  the  current  equation,  which  for  electrons  is  of 
the  form 

=  (2) 

and  for  holes; 

=  (3) 

Here,  n  and  p  denote  electron  and  hole  concentration,  respectively,  v 
velocity,  and  D  diffusivity.  The  usual  derivations  of  Eqs.  (2)  and  (3)  proceed 
from  a  linearization  of  the  Boltzmann  transport  equation.  The  assumption  is 
then  made  that  the  equation  is  valid  for  high-field  nonlinear  transport. 
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Typically,  the  field-dependent  velocity  assumed  in  these  equations  is  of  the 
type  displayed  in  Fig.  1.2. 

While  the  use  of  the  field-dependent  velocity  in  these  equations  is 
universal,  the  type  of  diffusivity  coefficient  used  in  these  studies  is  almost  as 
numerous  as  the  numbers  of  workers  involved  in  numerical  studies.  How¬ 
ever,  a  number  of  important  issues  are  at  stake  in  the  description  of  the 
diffusivity.  For  example,  if  the  Einstein  relation 


D  = 


jxkT 

e 


(4) 


is  used,  then,  under  equilibrium  and/or  zero  current  conditions,  the  depen¬ 
dence  of  carrier  density  on  conduction  and  valence  band  energy  is  given  by 
either  the  equilibrium  Boltzmann  or  Fermi  distribution.  However,  under 
nonequilibrium  conditions  (and  near-zero  current  conditions),  the  Einstein 
relation  inadequately  describes  diffusive  transport  [13].  To  correct  for  the 
latter  deficiency,  the  field-dependent  diffusivity  often  used  in  calculations  is 
of  a  form  similar  to  that  shown  below  [14]: 


thl 

e 


-I-  TV 


sat 


(5) 


where  at  high  values  of  electric  field,  the  diffusivity  only  gradually  de¬ 
creases.  While  the  diffusivity  coefficient  of  Eq.  (5)  more  adequately  repre¬ 
sents  high-field  phenomena,  because  its  field  dependence  is  conceptually 
consistent  only  with  the  assumption  of  nonequilibrium  conditions  it  is 
conceptually  inconsistent  with  equilibrium  conditions,  and  will  lead  to 
incorrect  built-in  potentials  [15]. 

While  the  drift  and  diffusion  equations  (DDE)  clearly  offer  conceptual 
difficulties  with  respect  to  consistency  of  physics,  they  nevertheless  offer 
considerable  insight  into  the  physics  of  device  operation  and  are  useful 
providing  their  limitations  are  kept  in  mind.  For  example,  instabilities  in 
long  GaAs  structures  are  known  to  depend  critically  on  conditions  at  the 
contacts.  A  study  in  1969  [16]  demonstrated  that  by  experimentally  creating 
different  conditions  at  the  boundaries  to  the  active  region  of  GaAs,  a  wide 
range  of  different  electrical  instabilities  could  be  obtained.  Corresponding 
numerical  studies  were  performed  through  solutions  to  the  above  drift  and 
diffusion  equations,  in  which  a  value  for  the  electric  field  was  specified  at  the 
cathode  (and  anode)  boundary.  It  was  found  that  the  boundary-dependent 
electrical  behavior  could  be  broken  into  three  categories,  as  summarized  in 
Fig.  1.6  [16].  The  key  conclusion  of  the  study  was  that  the  electrical 
behavior  of  compound  semiconductors  devices  was  dependent  in  a  detailed 
way  on  contact  conditions.  This  same  critical  result  has  reappeared  numer¬ 
ous  times  in  a  variety  of  different  types  of  structures. 


I 
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1.4  NONEQUILIBRIUM  DESCRIPTIONS  OF  TRANSPORT 

The  situation  of  most  interest  lies  in  nonequilibrium  transport.  The  most 
critical  area  of  interest  is  the  incorporation  of  acceleration  into  the  govern¬ 
ing  equations. 

In  examining  nonequilibrium  transport,  several  approaches  have  been 
used.  One  is  the  Monte  Carlo  method,  where  the  trajectory  of  a  particle  is 
followed  through  its  acceleration  and  subsequent  scattering  events.  In  the 
discussion  below,  results  of  Monte  Carlo  calculations  will  be  presented,  but 
we  first  concentrate  on  nonequilibrium  phenomena  as  described  by  the 
moments  of  the  Boltzmann  transport  equation.  These  equations,  in  their 
simplest  form  for  parabolic  bands,  a  position-dependent  conduction 
band,  and  a  position-dependent  effective  mass,  take  the  form  shown 
below  [17]: 

Carrier  balance: 


—  +  V,  •  /i - 

dt  m 


2 

(Itt)'  J  dt 


d\ 

coll 


Momentum  balance: 


(6) 


—  nhk.  -1-  V,  • 
dt  ^ 


m 


=  -nV^E^  +  qny  x  B  —  V^nkT 


+ 


\  2m  2  /  m  J  ^tlcoii  ^  ' 


Energy  balance: 


dt 


ttV  J  dt 


rk.-k 


+ 1  ytr 


2m 


^^k*k 


coll  2m 


d\ 


In  the  above,  hkj  is  the  mean  momentum  of  the  carriers,  T  is  the  carrier 
temperature,  and,  for  electrons,  is  the  position-dependent  conduction 
band  energy.  B  is  an  applied  magnetic  field.  The  terms  on  the  right  side 
represent  scattering  and/or  electron-hole  interaction,  as  through  avalanch¬ 
ing.  For  example,  the  right  side  of  Eq.  (6)  represents  intervallcy  scattering. 
If  avalanching  occurs,  generation  is  expressed  through  an  erergy-dependent 
ionization  coefficient  [18].  If  a  carrier  temperature  model  is  assumed,  then 
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carrier  generation  is  given  by 

na{T)  (9) 

where  a(T)  is  the  ionization  coefficient.  In  the  absence  of  a  first-principle 
determination  of  oi(T),  the  following  relation  can  be  assumed  as  a  starting 
point; 

a(T)  =  a*(F)v(F)  (10) 

where  oc*(F)  and  v(F)  are  the  equilibrium  ionization  rates  and  field- 
dependent  velocities,  respectively  and  the  relation  between  T  and  F  is 
determined  from  the  equilibrium  solution.  While  the  Eq.  (10)  relation  is 
uncertain,  it  has  the  conceptual  advantage  of  relating  ionization  to  energy, 
rather  than  field. 

But,  perhaps  the  most  significant  feature  of  these  equations  is  the 
presence  of  acceleration,  both  spatial  and  temporal  in  the  momentum 
balance  equation.  These  acceleration  terms  are  absent  from  the  drift  and 
diffusion  equations.  Additionally,  under  equilibrium  conditions,  and  hence, 
zero  current  (i.e.,  n^k^/m  =  0),  the  electron  temperature  model  teaches 
that  for  any  spatially  nonuniform  structure,  such  as  a  p-n  junction,  the 
electron  temperature  is  everywhere  constant  and  equal  to  the  ambient. 
Thus,  CO-  ceptual  problems  arising  from  the  form  of  the  diffusion  contribu¬ 
tion  to  the  drift  and  diffusion  equations  do  not  enter  here.  Note  that  a 
generalized  drift  and  diffusion  current  term  is  obtained  when  the  left  side  of 
Eq.  (7)  is  set  to  zero. 

Equations  of  the  type  shown  above  provide  a  considerable  amount  of 
information  with  regard  to  transport.  For  example,  with  a  T-L-X  orienta¬ 
tion  in  GaAs  the  distribution  of  carriers  as  a  function  of  field  is  shown  in 
Fig.  1.7  [5].  Here  the  relative  distribution  of  carriers  in  each  of  the  valleys  is 
determined  by  the  distribution  of  temperature  in  each  of  the  valleys,  which 
in  turn  is  driven  by  the  electric  field,  as  shown  in  Fig.  1.8  [5].  Note  that  for 
fields  below  4kV/cm,  the  carriers  reside  in  the  F-valley.  At  fields  above 
6kV/cm,  the  L- valley  population  exceeds  that  of  the  F-valley.  It  should  be 
emphasized,  however,  that  because  of  the  very  low  subsidiary  valley  mobili¬ 
ty,  most  of  the  current,  for  fields  up  to  50  kV/cm,  is  carried  by  the  F-valley 
carriers. 

The  interest  in  nonequilibrium  equations  lies  not  in  the  steady-state 
uniform  field  distribution,  but  in  transients  and  nonuniform  fields.  The 
transient  distribution  of  carrier  density  and  velocity  for  electrons  subject  to  a 
sudden  change  in  electric  field  is  shown  in  Fig.  1.9  [5],  where  the  high  peak 
velocity  can  be  noted. 

The  high  peak  velocity  in  Fig.  1. 9a  is  primarily  associated  with  F-valley 
transport.  Indeed,  the  F-valley  velocity  at  2.8  ps  is  near  4  x  lO’  cm/s.  The 
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Figure  1.7  Steady-state  uniform  field  carrier  distribution  for  T-L-X  orientation  in 
GaAs  [5]. 


Figure  1.8  Temperature  dependence  of  electric  field  for  a  T-L-X  orientation  in 
GaAs.  Electric  field  is  the  independent  variable  [5]. 


14  INTRODUCTION  TO  THE  PHYSICS  OF  GALLIUM  ARSENIDE  DEVICES 


T/(.28ps) 

Figure  1.9  Transient  overshoot  for  a  field  of  9.6kV/cm  [5].  T-L-X  orientation, 
and  an  applied  field  of  (a)  F  and  mean  velocity,  (b)  F-  and  L-valley  carrier  density. 

mean  velocity 

i;  =  (rZpUp  + (11) 

where  n  denotes  net  population  of  the  F-,  Z.-,  and  A'-levels  and  Nq  is  the 
total  carrier  density,  is  also  shown  in  Fig.  1.9.  We  note  that  the  significant 
drop  in  mean  velocity  is  a  consequence  of  electron  transfer  from  the  central 
to  the  satellite  valleys  (see  Fig.  1.9b).  Also  note  a  transient  decrease  in 
F-valley  velocity.  This  is  a  consequence  of  the  difference  between  the  energy 
(longer)  and  momentum  (shorter)  relaxation  times  in  GaAs.  The  time- 
independent  spatially  nonuniform  situation  also  displays  overshoot  effects. 

The  situation  when,  in  one  dimension,  space  charge  effects  are  intro¬ 
duced  is  displayed  in  Figs.  1.10  and  1.11  [5],  where  for  a  gallium  arsenide 
device  of  different  lengths  we  show  the  field  and  space  charge  distribution  of 
the  F-valley  electrons  and  the  current-voltage  characteristics.  The  feature  to 
note  is  that  as  the  device  length  decreases  the  current  drive  increases.  Note 
that  in  all  cases  the  field  is  nonuniform  and  increases  toward  the  anode. 
Electron  transfer  exists  for  all  four  structures,  with  the  greatest  amount  of 
transfer  occurring  for  the  longest  device.  Additionally,  since  high-field 
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DISTANCE /L 

Figure  1.10  Electric  field  and  F-valley  carrier  distribution  for  GaAs  (with  a 
two-level  model)  two-terminal  devices  of  different  lengths  and  a  mean  field  of 
5kV/cm  [5], 


F  (kV/cm) 
o 

Figure  1.11  Current-voltage  characteristics,  as  a  function  of  lengths,  for  the  struc¬ 
tures  of  Fig.  1.10  [5]. 
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Figure  1.12  Dimensions  and  doping  level  of  the  simu 
lated  PBT  [20], 

values  are  synonymous  with  carrier  accumulation,  we  see  that  electron 
transfer  here  is  accompanied  by  an  accumulation  of  L-valley  carriers. 
Saturation  in  the  current  density  occurs  at  high  bias,  and  even  for  the 
shortest  device  there  is  electron  transfer  at  the  anode  side  of  the  device. 

The  role  of  nonequilibrium  transport  on  two-dimensional  simulations  is 
discussed  for  the  vertical  three-terminal  GaAs  permeable  base  transistor 
(19],  one  cell  of  which  is  displayed  in  Fig.  1.12  [20].  One  important 
advantage  of  this  structure  is  the  parallel  placement  of  the  source  and  drain 
contacts*  and  the  absence  of  any  substrate  through  which  current  can  flow 
and  reduce  the  transfer  characteristics  of  the  device.  The  simulations  were 
performed  for  a  1  fim  long  source  to  drain  region  and  a  200  A  gate.  Also, 
for  comparison,  results  of  the  drift  and  diffusion  equation  simulations  were 
included. 

The  computed  I-V  characteristics  of  the  device  shown  in  Fig.  1,12  are 
presented  in  Fig.  1.13  [20].  The  moments  of  the  Boltzmann  transport 
equation  results  are  extrapolated  to  the  origin,  as  indicated  by  the  long 
broken  lines.  The  shorter  dashed  curves  show  the  results  for  the  DDE.  The 
comparison  shows  that  the  predicted  current  levels  are  significantly  higher 
for  the  moment  equation  solutions,  a  result  consistent  with  FET  calculations 
performed  by  Cook  and  Frey  [21],  who  used  a  highly  simplified  momentum- 
energy  transport  model.  TTie  present  calculation  results  also  indicate  a 
region  of  negative  differential  forward  conductivity  at  =  0.6  V.  The 
origin  of  this  phenomena  is  a  <»nsequence  of  electron  transfer.^  The 
presence  of  a  dc  negative  forward  conductance  is  also  a  feature  of  PBT 
measurements  but  is  clearly  absent  from  DDE  simulations  [20]. 

‘(Editor’s  Note)  They  are  also  called  emitter  and  collector  in  PBT  terminology. 

’(Editor’s  Note)  See  also  Chap.  2,  Section  2.2.4  for  additional  views  on  the  issue  of 
negative  differential  conductance. 


EMITTER 
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Figure  1.13  Collector  current  versus  collector-emitter  voltage  for  different  values 
of  base-emitter  voltage  [20]. 

A  comparison  between  the  total  carrier  density  distribution  along  the 
center  of  the  channel  for  drift  and  diffusion  and  MBTE  solutions  is  shown  in 
Fig.  1.14  [20]  for  Fee  =  1.0  V  and  Vbe=0.4V.  The  moment  equation 
prediction  for  the  F-viley  carrier  density  is  also  shown. 

As  seen  in  Fig.  1.14  for  the  DDE  simulations,  the  carrier  density  reaches 
a  maximum  between  base  contacts  and  displays  a  significant  dipole  layer. 
Here,  with  the  velocity  in  saturation  and  the  cross-sectional  area  at  a 
minimum,  the  carrier  density  must  increase  to  maintain  current  continuity. 
In  the  moment  equation  simulation,  the  constraints  of  current  continuity  are 
more  complex.  First,  a  decrease  in  the  cross-sectional  areas  is,  as  in  the 
DDE,  accompanied  by  an  increase  in  field  along  the  channel.  The  field 
increase  under  both  equilibrium  and  nonequilibrium  conditions  is  qualita¬ 
tively  similar,  as  may  be  observed  from  Fig.  1.15  [20],  which  shows  the 
potential  distribution  along  the  center  of  the  PBT  channel.  However, 
consequent  changes  in  electron  temperature,  both  increasing  and  decreas¬ 
ing,  lag  behind  the  equiUbrium  state.  This  leads  to  velocity  overshoot  (Fig. 
1.16  [20])  and  a  delay  in  electron  transfer.  As  a  result,  for  nearly  the  first 
half  of  the  device,  transport  is  almost  exclusively  F-valley  transport.  The 
implication  is  that  if  the  F-valley  carrier  velocity  increases  with  increasing 
field,  then  the  product  of  density  and  cross-sectional  area  normal  to  current 
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Figure  1,16  F-valley  velocity  along  the  center  of  the  PBT  channel  [20]. 


tions  the  L- valley  carriers  make  a  negligible  contribution  to  current.  Thus,  a 
decrease  in  carrier  velocity  results  in  a  net  increase  in  carrier  concentration. 
However,  initially,  the  decrease  in  field  is  not  accompanied  by  a  correspond¬ 
ing  temperature  decrease  (as  experienced  in  the  uniform  field  calculations). 
Thus,  the  high  F-valley  temperature  results  in  transfer  to  the  L-valley, 
giving  rise  to  the  second  minimum  in  the  F-valley  carrier  density  shown  in 
Fig.  1.14.  Further  toward  the  drain,  the  field  decreases.  However,  relaxation 
is  incomplete  and  the  field  at  the  collector  is  not  equal  to  the  field  at  the 
emitter.  Also  note  that  the  moment  equation  potential  distribution  gives  rise 
to  a  slightly  higher  field  upstream  of  the  base,  and  a  lower  field,  over  a 
longer  distance,  downstream  of  the  base  compared  with  the  drift  and 
diffusion  result.  More  significantly,  the  electron  temperature  at  the  collector 
exceeds  that  at  the  emitter.  It  is  noted  that  as  the  field  relaxes,  the  electrons 
transfer  back  to  the  central  valley. 


1.5  NONEQUILffiRIUM  ELECTRON-HOLE  TRANSPORT 

Additional  nonequilibrium  studies  were  mentioned  at  the  beginning  of  this 
chapter.  This  concerns  nonequilibrium  electron-hole  transport,  which  for 
specificity  was  discussed  for  InGaAs.  The  details  are  considered  below. 
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In  this  recent  study,  ensemble  Monte  Carlo  studies  were  performed  in 
which  electrons  were  injected  into  p-type  InGaAs  [8].  In  one  case,  the 
acceptor  doping  was  10*’  cm~^  and  in  the  second  case  5  x  10*’  cm~  .  The 
calculations  were  at  300  K  and  the  ratio  of  the  injected  electrons  to  the 
majority  holes  was  taken  to  be  1:5.  (Note:  The  ensemble  Monte  Carlo 
avoids  any  assumptions  on  the  magnitudes  of  the  energy  and  momentum 
exchange  in  an  electron-hole  [e-h]  scattering  process  and  the  evolution  of 
the  electron  and  hole  distribution  functions.)  Tlie  electrons  and  holes  were 
assumed  to  be  in  equilibrium  with  the  lattice  when  the  electric  field  was 
switched  on,  and  the  band  model  consisted  of  three  nonparabolic  valleys  for 
the  conduction  band  and  two  parabolic  light-  and  heavy-hole  bands.  The 
role  of  the  light  holes  was  suppressed  in  this  study.  The  model  includes  the 
elastic  acoustic  phonon,  impurity  scattering  using  Ridley’s  model,  alloy 
scattering,  deformation  potential,  and  intervalley  and  intravalley  phonon 
scattering  process.  The  e-h  and  screened  carrier-phonon  scattering  are 
calculated  from  the  expressions  given  in  Osman  and  Ferry  [22],  using  a 
self-consistent  screening  model.  There  is  only  one  LO  phonon  mode  in  this 


Figure  1.17  (a)  Electron  drift  velocity  in  p-type  InGaAs,  (b-c)  percentage  of 
upper-valley  electrons  in  p-type  InGaAs  [fci  (Reprinted,  with  permission,  from  I 
Applied  Physics  Letters.) 
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calculation.  The  electron  transport  parameters  for  In^  sjOa^  47  As  are  the 
same  as  those  reported  in  Ahmed  et  al.  (23).  However,  for  hole  effective 
masses  and  deformation  potential  constants,  appropriate  to  GaAs  were 
used.  The  interaction  between  the  L-valley  electrons  and  the  heavy  holes 
was  ignored. 

The  drift  velocities  of  the  electrons  injected  into  the  p-type  InGaAs  as  a 
function  of  the  applied  electric  fields  are  plotted  in  Figs.  1.4  and  1.17a  [8], 
for  doping  levels  of  5  x  10*’  and  lO”  cm"^  respectively.  The  curves  connect¬ 
ing  the  open  circles  in  these  figures  correspond  to  situations  in  which  the 
interaction  with  the  mobile  holes  is  ignored  and  only  the  interaction  with  the 
ionized  acceptor  impurities  is  taken  into  account. 

When  the  interaction  between  the  minority  electrons  and  the  mobile 
majority  holes  is  taken  into  account,  the  electron  velocities  are  lower  for 
fields  below  4  kV/cm  compared  with  majority  electrons,  as  can  be  seen  from 
these  figures.  At  these  low  fields,  the  electron  transfer  to  the  upper  valleys  is 
negligible,  as  shown  in  Fig.  1.17b  [8],  and  the  energy  loss  through  e-h 
interaction  is  not  significant  because  the  rate  at  which  the  electrons  gain 
energy  from  the  electric  field  is  small,  as  can  be  seen  in  Fig.  1.18  [8]. 
However,  the  e-h  scattering,  which  has  the  same  angular  distribution  as  the 
impurity  scattering,  has  the  same  effect  on  the  electron  mobility  as  doubling 
the  doping  level.  Consequently,  the  mobility  of  the  electrons  is  reduced, 
leading  to  lower  velocities.  As  the  electric  field  is  increased,  the  fraction  of 
electrons  with  enough  energy  to  transfer  to  the  upper  valleys  increases. 


Figure  1.18  Electron  energy  in  fvtype  InGaAs  as  a  function  of  field  [8]. 
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For  instance,  from  th  equation  of  motion  of  the  density  matrix,  for  a 
system  including  mobile  carriers  and  scattering  centers,  the  first  three 
moment  equations  have  the  following  form  [27]: 

«?'“’»  +  iP'"}  =  I  «[H„  (12) 

^mTx  S  «("■’  -^"’1^ 
where  the  {(  ))  denote  the  quantum  ensemble  averages,  and  using  Dirac 


notation,  the  operators  of  interest  are  of  the  form 

/><“'  =  |a:„  X  j:„|  (15) 

/><‘'  =  (|)(P|jr„Xx„|  +  |j:„XAr„|P)  (16) 

pm  =  (  j  X  +  2/.|x„  X  xji>  +  |x„  X  x„|P")  (17) 

pm  =  (  J  )=(/>=|x„  X  x„|  +  3P|x<,  X  x„|  P"  +  3P"|x„  X  x„lP 

+  |x„xxo|P’)  (18) 


where  P  is  the  momentum  operator.  Note  that  the  terms  involving 
incorporate  dissipation.  In  a  diagonal  representation,  the  ensemble  average 
of  the  first  three  operators  breaks  down  into  the  following  form: 

(/p(‘’)))=  2  p„n..(Xo)  =  n(xo) 

=  n{xo)mv^ 

iP^^'^))=  E  ^  E  Pi./*/  ^  In  tii 

=  +  m\v\  -  E  P,.n,  In  n^ 

where  p^^  is  the  diagonal  element  of  the  density  matrix,  and 

=  E  PunM^i  - 


(19) 

(20) 

(21) 

(22) 


It  is  clear  that  with  the  exception  of  the  third  operator,  which  contains  a 
term  involving  Planck’s  constant,  the  equations  are  of  a  classical  form.  Thus, 
in  the  simplest  approximation,  there  appears  to  be  a  close  similarity  between 
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the  classical  moment  equations  and  that  obtained  quantum  mechanically. 
The  difficulty  is,  of  course,  in  solving  these  equations. 

There  is,  however,  an  interesting  situation  to  consider:  that  in  which 
p„  =  1  /N  the  system.  In  this  case,  the  first  two  moment  equations,  including 
dissipation  in  momentum,  reduce  to 


+  div{nv)  =  0 


+  —  nv.  +  —  = 


\  dx 


dV ^  dQ\  n 
dx  dx  /  m 


where 


Q  =  - 


^  1  d  Vn 

2m  'dH 


The  quantity  V  represents  an  imposed  barrier  and  the  self-consistent  energy 
associated  with  Poisson’s  equation.  The  potential  Q  [24]  is  density-depen- 
dent  and  tends  to  become  significant  near  strong  barriers,  where  the 
curvature  of  V«  will  either  enhance  or  diminish  the  imposed  barrier. 
Tunneling  and  resonance  arise  from  Q. 

In  mr'tiple  dimensions,  these  equations  are  subject  to  the  constraint 

^  mv '  dx  —  nh  (26) 

or  in  a  gauge  that  includes  a  vector  potential,  the  constraint 

(f  (mv  +  -  A)  •  dx  =  nh  (27) 


Contact 


Contact 


Figure  1.20  Configuration  s-  itable  for  (he  Aharonov-Bohm  constraint  [29].  (Re 
printed,  with  permission,  from  Applied  Physics  Letters.) 
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Presently,  device  systems  are  being  constructed  which  are  influenced  by  the 
constraint  of  Eq.  (27),  often  called  the  Aharonov-Bohm  condition  [28],  In 
particular,  structures  are  being  constructed  in  which  the  path  of  an  incident 
beam  of  electrons  is  split  and  then  recombined.  The  split  path  lengths  of  the 
original  beam  are  different,  and  under  coherent  reconstruction  in  which  Eq. 
(27)  is  satisfied,  conduction  oscillations  are  anticipated.  A  structure  original¬ 
ly  proposed  to  deal  with  this  is  displayed  in  Fig.  1.20  [29]. 
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As  device  dimensions  in  nanoscale  structures  and  mesoscopic  systems 
are  reduced,  the  characteristics  and  interactions  of  dimensionally-confined 
longitudinal-optical  (LO)  phonons  deviate  substantially  from  those  of  bulk  polar 
semiconductors.  This  account  emphasizes  LO-phonon  effects  arising  in  three 
separate  systems:  Short-period  AlAs-GaAs  superlattices;  rectangular  GaAs 
quantum  wires  embedded  in  A1  As;  and  metal-semiconductor  interfaces  such  as 
those  in  de  Broglie  wave  quantum-interference  devices. 

I.  INTRODUCTION 

In  nanoscale  and  mesoscopic  systems,  the  effects  of  confinement  on 
carriers  have  been  studied  extensively.  However,  to  properly  model  carrier 
energy  loss  in  nanoscale  and  mesoscopic  systems,  it  is  essential  that  calculations 
of  carrier  scattering  by  longitudinal-optical  phonons  take  into  account  the  fact 
that  confinement  also  changes  the  strength  and  spatial  properties  of  longitudinal- 
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optical  phonons.  As  a  first  example  of  phonon-confinement  effects,  we  shall 
emphasize  the  dominant  role  of  interface-optical  phonons  on  electron  scattering 
in  heterostructures  with  confinement  lengths  less  than  about  40  Angstroms. 

Electron  interactions  with  longitudinal-optical  (LO)  phonon  modes  in 
heterostructures  are  strongly  affected  by  the  changes  in  theFrohlich  Hamiltonian 
caused  by  phonon  confinement  and  localization,  as  well  as  by  the  changes  in  the 
electronic  wave  function  due  to  the  confining  potential.  The  presence  of 
heterointerfaces  produces  large  changes  in  the  dielectric  constant  and  gives  rise 
to  the  confinement  of  optical  phonons  in  each  layer  (i.e.,  confines  modes)  and  the 
localization  in  the  vicinity  of  interfaces  (i.e.,  interface  modes  or  surface-optical 
(SO)  modes).  These  confined  and  interface  modes  arise  because  the  abrupt 
changes  in  dielectric  constant  near  heterointerfaces  make  it  impossible  to  satisfy 
the  normal  bulk  dispersion  relation  which  normally  restricts  LO  phonon  frequen¬ 
cies  to  those  where  the  dielectric  constant  vanishes.  There  have  been  suggestions 
that  interactions  by  the  interface  modes  can  be  significant  and  that  the  scattering 
rate  due  to  the  confined  modes  can  be  considerably  reduced  in  some  structures 
compared  to  the  bulk  LO-phonon  scattering  rate.  1-3  Therefore,  an  appropriate 
treatment  of  the  optical-phonon  modes  in  quantized  systems  is  essential  for 
understanding  electron  transport  in  heterostructures.  Recently ,  both  macroscopic 
and  microscopic  approaches  to  electron-optical-phonon  interactions  in 
heterostructures  i-n  have  been  applied  in  theoretical  treatments.  Enhanced 
electron-SO-phonon  (i.e.,  electron-interface-phonon)  scattering  in  polar  semi¬ 
conductors  with  confining  dimensions  less  than  about  50  Angstroms  has  been 
indicated  recently  12.13  on  the  basis  of  scalar-potential  modes.  Physically,  this 
enhanced  scattering  rate  is  due  to  the  interface  phonon  potential  which  increases 
approximately  exponentially  near  heterointerfaces,  where  there  is  an  abrupt 
change  in  the  dielectric  constant.  The  appropriateness  of  using  a  scalar  potential 
to  model  the  electron-SO-phonon  interaction  has  been  questioned.i4  Herein,  we 
model  the  potential  to  calculate  the  ratio  of  electron-SO-phonon  scattering  to 
electron-confined-LO-phonon  scattering  in  a  short-period  GaAs/AlAs 
superlattices.  15  The  results  provide  a  criterion  which,  in  conjunction  with 
experiments,  can  be  used  to  examine  the  validity  of  the  scalar-potential  model. 

n.  MODEL  AND  PROCEDURE 

Figure  1  shows  a  schematic  drawing  of  the  GaAs/AlAs  short-period 
superlattice  considered  in  this  study.  Each  GaAs  layer  has  width,  a,  and  each 
AlAs  layer  has  width,  b.  We  will  consider  the  case  where  a  uniform  electron 
current,  J,  is  incident  parallel  to  the  superlattice  heterojunctions.  The  ratio  of 
electron-SO-phonon  scattering  to  electron -LO-phonon  scattering  in  the  superlattice 
of  Fig.  1  is  calculated  through  the  Fermi  golden  rule  by  treating  the  electron-Lo- 
phonon  and  electron-SO-phonon  interaction  Hamiltonians  as  perturbation 
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FIG  1.  Short-period  GaAs/AIAs  superlattice  with  a  uniform  electron  curent 
parallel  to  the  heterojunction  interfaces. 

Hamiltonians  and  by  taking  the  electronic  wavefunction  in  the  z-direction  as  that 
given  by  Cho  and  Prucnal  for  the  maximum  edge  of  the  first  subband.16 

To  obtain  an  order-of-magnitude  estimate  of  the  ratio  of  electron-SO- 
phonon  scattering  to  electron-LO-phonon  scattering  in  theGaAs/Al  As  superlattice 
of  Fig.  1,  we  take  the  ratio  of  the  Fermi  golden  rule  transition  probability  for 
electron  scattering  from  the  symmetric  AlAs-likc  interface  mode  to  the  corre¬ 
sponding  transition  rate  for  electron  scattering  from  the  lowest  order  confined 
phonon  mode. 

In  Fig.  2,  this  ratio  is  plotted  as  a  function  of  the  parallel  phonon 
wavevector  times  the  quantum  well  width,  qa,  for  a  GaAs/AIAs  short-period 
superlattice  with  a=b.  For  qa  <  it,  the  ratio  is  larger  than  unity.  A  principal  factor 
contributing  to  the  rapid  increase  for  small  values  of  aq  (=X)  is  the  exponential 
increase  in  the  strength  of  the  SO-phonon  mode  near  the  heterojunction.  In 
agreement  with  our  results,  a  surge  in  the  strength  of  the  S+  mode  with  a 
decreasing  well-width,  a,  is  observed  in  a  recent  low-temperature  experiment,!^ 
which  clearly  demonstrates  the  dominance  of  the  SO-phonon  modes  over  iheLO- 
phonon  modes.  These  results  are  based  on  the  implicit  assumption  that  the 
uniform  current,  J,  is  used  to  prepare  the  electronic  state  given  by  Ref.  16.  This 
assumption  is  reasonable  if  the  electrons  have  incident  energies  signiHcantly  less 
than  the  GaAs/AIAs  conduction  band  discontinuity  and  if  the  superlattice  b^er 
widths  are  roughly  50  Angstroms  or  less.  Furthermore,  the  frequency  of  the 
symmetric  AlAs-like  mode  foraGaAs/AlAs  heterojunction  varies  from  about45 
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FIG.  2.  Approximate  ratio  of  electron-SO-phonon  (S+  mode)  scattering  to  elec- 
tron-LO-phonon  (m=1  confined  mode)  scattering. 

meV  for  qa  >  6  to  50  meV  near  qa=0  and  the  frequency  of  the  lowest-order 
confined  mode  is  approximately  36.2  meV.  Thus,  the  scaling  illustrated  in  Fig. 
2  is  expected  to  hold  for  a  <  50  Angstroms  and  for  incident  electron  energies  in 
the  range  from  about  50  meV  to  very  roughly  150  meV.  In  this  eno-gy  range,  the 
phonon  momentum,  q,  will  be  approximately  equal  to  the  incident  electT''n 
momentum.  Thus,  we  conclude  that  enhanced  SO-phonon  scattering  will  oc  r 
in  the  regime  of  low-voltage  (also,  low-temperature  in  general)  operation 
characteristic  of  mesoscale  devices  currently  being  considered  by  the  electron 
device  community. 


As  a  second  example,  of  the  nature  of  Lo-phonon  confinement  effects 
in  polar  semiconductor  nanostructures,  we  consider  Golden  rule  scattering  rates 
for  emission  and  absorption  in  rectangular  GaAs  quantum  wires  embedded  in 
AlAs.  In  these  calculations,  the  carriers  are  taken  to  be  free  in  the  x-direction 
and  bound  in  the  ground  states  of  finitely  deep  quantum  wells  in  the  y-  and  z- 
directions.  For  a  40- Angstrom-by-40- Angstrom  GaAs  quantum  wire  of  infinite 
length  in  the  x-direction,  we  find  strong  one-dimensional  dcnsity-of-states  effects 
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in  the  rales  for  surface-optical  (SO)  phonon  emission.  The  emission  peaks  for  the 
SO-phonons,  also  known  as  interface-optical  phonons,  occur  near  the  AlAs-like 
and  GaAs-like  LO-phonon  frequencies.  The  confined  and  interface  phonon 
modes  used  in  obtaining  these  results  are  those  considered  previously  by  the 
authors.6.9.12.21  pig.  3  depicts  the  Golden  rule  scattering  rates  for  the  quantum 
wire.  The  strong  SO-phonon  emission  rates  near  the  AlAs-like  and  GaAs-like 
LO-phonon  frequencies  are  more  than  an  order  of  magnitude  larger  than  the 
corresponding  bulk  LO  phonon  emission  rates.9.21  However,  if  the  lateral 
dimensions  of  the  quantum  wire  are  varied  sinusoidally  in  amplitude  by  only  ten 
percent,  the  narrow  one-dimensional  density-of-staies  emission  peaks  have 
averaged  values  close  to  those  of  bulk  carrier — ^LO-phonon  interaction  rates.21 
In  this  case,  we  learn  that  carrier — ^LO-phonon  scattering  rates  may  deviate 


FIG.3.  Scattering  rates  at  300  K  as  a  function  of  electron  energy  in  the 
longitudinal  direction.  The  size  of  quantum  wire  is  40  Angstroms  by  40  Angstroms. 
The  SO-phonon  scattering  rate  is  denoted  as  SO,  f  forthe  free-standing  quantum 
wire,  and  as  SO  for  the  GaAs/AIAs  quantum  wire,  respectively,  while  LO 
represents  the  rate  of  electron  scattering  with  confined  LO-phonon  modes. 
Phonon  emission  (absorption)  process  Is  represented  as  emi  (abs),  respectively. 
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substantially  from  the  corresprinding  wires;  however,  in  the  more  realistic  case 
of  variable  quantum  wire  dimensions,  one-dimensional  effects  arc  suppressed. 

As  a  third  and  final  example  of  how  optical  phonon  properties  in 
nanoscale  structures  may  deviate  substantially  from  those  of  bulk  optical 
phonons,  we  consider  a  metal-semiconductor  interface  for  the  case  where  the 
semiconductor  thickness  is  restricted  to  nanoscale  dimensions.  In  this  analysis, 
we  find  that  the  nature  of  the  interface  phonon  modes  at  the  metal-semiconductor 
interface  are  modified  substantially  for  two  reasons:  (1)  at  the  metal-semicon¬ 
ductor  interface,  only  those  phonon  modes  with  Frohlich  potentials  which 
approach  zero  at  the  interface  are  allowed;  and  (2)  the  magnitude  of  the  Frohlich 
potential  for  the  surviving  interface  phonon  mode  decreases  exponentially  with 
the  product  of  the  phonon  wavevector  and  thickness  of  the  semiconductor  slab. 
As  demonstrated  previously ,'*.6  the  components  of  the  optical  phonon  polariza¬ 
tion  vectors  for  dimensionally-confined  modes  have  opposite  parities  for  polar¬ 
ization  components  normal  to  the  semiconductor  interface  and  for  polarization 
components  parallel  to  the  interface.  Since  the  average  tangential  components 
of  the  optical-phonon  electric  field  must  vanish  at  the  metal  interface  and  since 
the  phonon  polarization  vector  is  proportional  to  the  phonon  electric  field,  it 
follows  that  only  optical  phonon  modes  with  even  polarization  vectors  normal  to 
the  metal-semiconductor  interface  will  survive.  Furthermore,  since  the  laplacian 
of  the  Frohlich  potential  for  a  given  phonon  mode  is  proportional  to  the 
divergence  of  the  polarization  vector  for  that  mode,  it  follows  that  only  optical 
phonon  modes  with  odd  Frohlich  potentials  satisfy  the  metal-semiconductor 
interface.  In  particular,  of  the  well-known  symmetric  and  antisymmetric 
interface  modes,2i  only  the  antisymmetric  mode  will  survive.  In  the  limit  of  no 
dispersion,  the  antisymmetric  interface  mode  has  a  Frohlich  potential  which 
scales  as  one  divided  by  the  square  root  of  cosh(qd)*sinh(qd)  where  q  is  the 
magnitude  of  the  phonon  wavevector  parallel  to  the  metal-semiconductor 
interface  and  d  is  the  thickness  of  the  semiconductor  slab.  For  small  values  of  qd, 
the  Frohlich  potential  of  the  antisymmetric  interface  mode  then  decreases 
exponentially  with  qd.  Hence,  we  conclude  that  the  carrier-optical-phonon 
scattering  rate  decreases  to  zero  at  the  metal -semiconductor  in  terface  and  that  the 
strength  of  the  Frohlich  potential  for  the  antisymmetric  interface  mode  in  the 
semiconductor  slab  decreases  exponentially  with  qd.  (It  is  important  to  note  that 
these  results  depend  upon  having  a  perfect  metal-semiconductor  interface  since 
even  a  thin  oxide  at  the  metal-semiconductor  interface  is  sufficient  to  isolate  the 
semiconductor  optical-phonon  modes  from  the  metal.)  The  impact  of  these 
observations  on  the  physics  of  nanoscale  structures  such  as  the  de  Broglie  wave 
interference  device  is  clear;  carrier-optical-phonon  scattering  will  be  reduced 
near  metal-semiconductor  interfaces.  Additionally,  Schottky  barrier  heights  at 
metal-semiconductor  interfaces  in  nanoscale  structures  may  deviate  from  those 
of  metal-semiconductor  structures  with  semi-infinite  metal  and  semiconductor 
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regions;  the  extent  of  this  deviation  will  depend  on  the  relative  importance  of 
phonon  scattering  and  other  factors  such  as  interface  states  which  determine  the 
equilibrium  Schottky  barrier  height. 

in.  CONCLUSION 

In  conclusion,  the  characteristics  and  interactions  of  dimensionally- 
confined  optical  phonon  modes  deviate  substantially  from  those  of  bulk  polar 
semiconductors.  In  particular,  theory  and  experiment  point  to  greatly  enhanced 
carrier — LO-phonon  scattering  rates  when  dimensional  constraints  optimize 
interface  mode  strengths.  For  thecase  of  enhanced  carrier — LO-phonon  emission 
rates  in  rectangular  quantum  wires,  our  analysis  indicates  the  one-dimensional 
densily-of-states  effects  may  be  greatly  reduced  as  a  result  of  minor  fluctuations 
in  the  laterial  dimensions  of  the  quantum  wire.  Finally,  we  conclude  that 
necessary  restrictions  in  phonon  symmetries  at  metal-semiconductor  interfaces 
will  modify  carrier-LO-phonon  scattering  rates  in  nanostructures  where  charge 
transport  occurs  in  thin  polar  semiconductor  regions  in  intimate  contact  with 
metals. 
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Abstract 

Novel  quantum-effect  polar-semiconductor  structures  underlie  technologies 
portending  dramatic  enhancements  in  the  capability  to  process  information  orders 
of  magnitude  faster  than  is  possible  currently.  In  many  embodiments  of  these 
quantum-effect  structures,  charges  are  transported  in  quasi-one-dimensional 
quantum  wires  which  must  support  the  transport  of  charges  at  high  mobilities. 
However,  it  has  recently  been  demonstrated  that  the  longitudinal-optical  (LO) 
phonons  established  at  quantum-wire  interfaces  lead  to  dramatic  enhancements  in 
carrier-phonon  interactions  and  concomitant  degradation  in  carrier  mobility.  This 
letter  demonstrates  that  phonon  modes  may  be  tailored  through  the  judicious  use  of 
metal-semiconductor  interfaces  in  such  a  way  as  to  dramatically  reduce  unwanted 
emission  of  interface  LO  phonons  and,  consequently,  to  lead  to  the  achievement  of 
high  quantum-wire  mobility. 


Modem  fabrication  techniques  for  making  semiconductor  structures  with 
nanometre-scale  dimensional  features  have  been  essential  in  leading  to  fundamental 
discoveries  (Esaki  1974,  Stein,  von  Klitzing  and  Weimann  1983)  as  well  as  in 
stimulating  concepts  for  future  information-processing  systems  based  on  the  exploit¬ 
ation  of  quantum  effects  occurring  in  such  structures  (Capasso  and  Datta  1990).  As 
originally  proposed  by  Sakaki  (1980),  the  predicted  high  mobilities  of  quasi-one- 
dimensional  wire-like  regions  of  a  semiconducting  material  underlie  many  proposed 
quantum-wire  system  concepts  (Luryi  and  Capasso  1985,  Sakaki  1989)  and  have 
resulted  in  pioneering  efforts  leading  to  the  actual  fabrication  of  quantum  wires  (Watt, 
Sotomayor  Torres,  Araot  and  Beaumont  1990).  Furthermore,  quantum-wire  struc¬ 
tures  have  been  essential  to  seminal  studies  underlying  quantum-wire  laser  concepts 
(Tsuchiya  et  al.  1989)  and  quantum-coupled  electronic  systems  (Reed  et  al.  1988). 
Recently,  however,  theoretical  studies  of  the  interaction  between  longitudinal-optical 
(LO)  phonons  and  carriers  in  a  polar-semiconductor  quantum  wire  (Stroscio  1989) 
have  revealed  the  presence  of  discrete  LO  phonon  modes  similar  to  those  identified  for 
polar-semiconductor  quantum  wells  (Kliewer  and  Fuchs  1966,  Licari  and  Evrard  1977, 
Mori  and  Ando  1989,  Kim  and  Stroscio  1990).  As  for  the  case  of  quantum  wells, 
interface  LO  phonons  are  established  at  the  semiconductor-semiconductor  bound¬ 
aries  of  quantum  wires  (Kim  et  a/.  1991).  In  addition,  it  has  been  shown  that,  for  carrier 
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energies  in  excess  of  the  interface  LO-phonon  energy,  the  inelastic  scattering  caused  by 
carrier-interface-phonon  interactions  dominates  over  other  scattering  mechanisms 
when  confinement  occurs  on  a  scale  for  about  40  A  or  less  (Stroscio  et  al.  1991).  In  this 
letter,  it  is  demonstrated  that  establishing  a  metal-semiconductor  interface  at  the 
lateral  boundaries  of  polar-semiconductor  quantum  wires  introduces  a  set  of  boundary 
conditions  that  dramatically  reduces  or  eliminates  unwanted  carrier  energy  loss  caused 
by  the  interactions  with  interface  LO-phonon  modes. 

Motivated  by  the  recently  demonstrated  technology  for  the  epitaxial  growth  of 
metals  in  intimate  contact  with  polar  semiconductors  (Harbison  et  al.  1988,  1989, 
Guivarc’h  et  al.  1989),  the  dielectric  continuum  model  of  interface-phonon  modes  has 
been  applied  to  determine  the  carrier-interface-phonon  interaction  Hamiltonian  near 
a  semi-infinite  metal-polar-semiconductor  interlace  (Stroscio  et  al.  1992).  Detailed 
microscopic  calculations  of  interface  modes  in  polar  semiconductors  (Rucher,  Molihari 
and  Lugli  1991)  indicate  that  the  dielectric  continuum  model  provides  an  accurate 
formalism  for  modelling  carrier-interface-phonon  interactions.  As  summarized  by 
Mori  and  Ando  (1989),  the  components  of  the  optical  polarization  vectors  for  confined 
and  interface  phonons  have  opposite  parities  for  polarization  components  normal  to 
the  semiconductor  interface  and  lo*  polarization  components  parallel  to  the  interface. 
Since  the  average  tangential  components  of  the  optical-phonon  electric  field  must 
vanish  at  the  metal  interface  (it  may  survive  at  a  semiconductor-semiconductor 
interface)  and  since  the  phonon  polarization  vector  is  proportional  to  the  phonon 
electric  field,  it  follows  that  only  optical-phonon  modes  with  even  polarization  vectors 
normal  to  the  metal-semiconductor  interface  will  survive.  Furthermore,  since  the 
Laplacian  of  the  Frohlich  potential  for  a  given  phonon  mode  is  proportional  to  the 
divergence  of  the  polarization  vector  for  that  mode,  it  follows  that  only  optical-phonon 
modes  with  odd  Frohlich  potentials  satisfy  the  metal-semiconductor  interface 
conditions.  Hence,  only  those  confined  and  interface  modes  of  odd  potential  satisfy  the 
correct  boundary  conditions  at  the  metal-semiconductor  interface.  In  particular,  of  the 
well  known  symmetric  and  antisymmetric  interface  modes,  only  the  antisymmetric 
modes  will  survive;  then  straightforwardly  from  the  formulation  of  Kim  and  Stroscio 
(1990),  it  follows  that 

H  =y  (  ftexp^L-^ 

^  q  \eo(5/d(u)[Ki(cu)  coth  (qd/2)  +  K2(m)] / 

1  w  ♦  .  sinh(gz)  ^  d 

where  denotes  the  antisymmetric  electron-interface-phonon  interaction  Hamil¬ 
tonian  in  a  polar-semiconductor  slab  of  dielectric  function  Kj(cu),  bounded  for  z<0 
with  a  semi-infinite  metal  and  for  z>d/2  with  a  semi-infinite  region  of  dielectric 
function  Kjfcu).  In  this  result,  q  is  the  component  of  the  interface-phonon  wave-vector 
parallel  to  the  semiconductor  heterojunctions,  p  is  the  pair  of  spatial  coordinates 
parallel  to  the  planes  of  the  semiconductor  heterojunctions,  z  is  the  coordinate 
orthogonal  to  p,  L  is  the  normalization  length,  a^  and  at.,  are  the  usual  phonon 
annihilation  and  creation  operators,  ard  Eq  is  the  permittivity  of  vacuum.  In  the 
dispersion  regions  where  q  does  not  :..pend  strongly  on  to,  1/^  scales  with  qd/l  as 
[cosh  {qd/2)  sinh  {qdl2)]  "  or  equivalently  as  2exp(— ^d/2);  thus,  in  the  dispersionless 
region  (away  from  ^  =  0),  the  strength  of  the  electron-interface-phonon  interaction 
decreases  exponentially  as  dl2  becomes  large.  The  impact  of  these  observations  on  the 
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physics  of  nanostructures  is  clear;  carrier-optical-phonon  scattering  will  be  reduced 
near  metal-semiconductor  interfaces. 

For  rectangular  polar-semiconductor  quantum  wires  embedded  in  a  metal,  the  only 
confined  and  interface-LO-phonon  modes  satisfying  the  metal-semiconductor  bound¬ 
ary  conditions  are  the  modes  with  odd  potentials  at  all  boundaries  of  the  quantum  wire. 
For  the  case  of  confined  modes,  the  electron-LO-phonon  Hamiltonian  is  precisely  that 
given  previously  by  Stroscio  (1989)  since  all  the  modes  derived  earlier  are  odd  at  the 
quantum-wire  interfaces;  thus  all  results  of  Stroscio  11989)  hold  for  the  case  of  a 
rectangular  polar-semiconductor  quantum  wire  embedded  in  a  metal.  For  the  case  of 
interface  modes,  there  are  no  allowed  solutions  since  it  is  impossible  to  satisfy  the 
metal-semiconductor  boundary  conditions  simultaneously  at  two  of  the  quantum- 
wire  metal-semiconductor  interfaces.  While  the  demonstration  holds  rigorously  only 
for  quantum  wires  with  rectangular  cross-sections,  it  is  clear  that  the  same  conclusion 
will  be  substantially  correct  for  a  quantum  wire  of  arbitrary  cross-section  since  the 
interface  modes  with  long-range  Coulomb  character  must  be  reduced  severely  by 
invoking  metal-semiconductor  boundary  conditions  at  two  opposing  locations  on  the 
quantum-wire  cross-section. 

The  conclusion  that  carrier-interface-phonon  interactions  will  be  absent  in 
quantum  wires  surrounded  by  metal  is  critical  in  maintaining  high-mobility  transport 
in  quasi-one-dimensional  quantum-wire  structures.  Indeed,  the  recent  theoretical  and 
experimental  findings  (Tsen  et  al.  1991)  indicate  that  the  strength  of  the  carrier- 
interface-phonon  scales  approximately  inversely  with  the  thickness  of  the  confining 
region.  However,  it  is  precisely  the  regime  of  the  smallest  dimensional  scales  which  has 
portended  the  greatest  advances  in  the  novel  application  and  discovery  of  quantum 
effects  in  semiconductor  structures.  As  demonstrated  in  this  letter,  this  highly 
unsatisfactory  situation  for  polar-semiconductor  quantum  wires  can  be  circumvented 
by  eliminating  interface-optical-phonon  modes  through  the  judidous  use  of  metal- 
semiconductor  interfaces. 
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Abstract  ~Thu  paper  presents  a  nets  hydrodynamic  transport  model  toilh  non- 
parabolie  conduction  bands  and  quantum  correction  t‘rms.  For  the  first  time 
solutions  for  the  full  quantum  balance  equations,  applied  to  an  ultrasmail  t,lee- 
Iron  device,  .re  presented. 

1.  Introduction 

Requirements  for  faster  electronics  have  produced  smaller  devices  influenced  by 
quantum  interference  efTects  which  can  not  be  modeled  by  classical  theory.  These 
devices  operate  under  nonequilibrium  conditions  where  the  average  carrier  energy 
reaches  many  times  its  equilibrium  value.  To  address  these  problems  several  ap¬ 
proaches  (1]  have  been  considered.  This  paper  describes  a  preliminary  investigation 
into  solutions  for  electron  density,  average  velocity  and  average  energy  for  an  ultra- 
small  electronic  structure  under  the  conditions  of  noiiclassical  electron  transport, 
using  an  approach  previously  described  as  quantum  hydrodynamics  [2]. 

2.  The  Nonparabolic  Model 

First,  we  present  a  brief  summary  of  c  unique  form  of  the  hydrodynamic  transport 
model  applicable  to  nonparabolic  conduction  bands.  A  complete  description  of  this 
model  will  be  formafued  elsewhere  (3).  Our  hydrodynamic  model  was  developed 
by  studying  moments  of  the  Boltzmann  transport  eeuation  [4].  This  process  was 
achieved  using  the  moment  operators  ♦i  =r  1.  ♦*  =  u(k)  =  and  ♦a  = 

jm(k)u(k)  •  u(k)  =:  £^e(k).  Here,  the  Kane  dispersion  relation,  =  Ec{l  + 

oEf).  has  been  used  to  define  u(k)  and  m(k)  =  m*  y/i  +  These  moment 

operators  were  chosen  because  they  lead  to  a  form  which  can  be  manipulated 
more  easily  and  one  in  which  simplifying  approximations  can  be  seen  more  clearly. 

Performing  the  moment  process  and  simplying  to  first  order  yields  immediately  the 
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results; 


Sn 

St 


=  -V,  (nv) 


—  =  -V  -  V,v  +  — - ;•  V,  •  (Pv 

at  TTi*  nm 

^  = -V  •  V,  tu  +  F  .  V  - -V,  ■(v-[Pwl  +  q) 
at  n 


(1) 

(2) 

(3) 


where  m**.  [Py].  (Pw]  *>*<1  q  have  integral  definitions  [3]  which  depend  on  /. 
At  this  point  additional  assumptions  or  relations  are  necessary  to  dose  the  mathe¬ 
matical  system  .  A  familiar  approach  is  to  use  a  displaced  Maxwellian  distribution 
function  [7],  However,  to  resolve  the  dilemma  of  two  distinctly  different  effec¬ 
tive  pressure  tensors  Py  and  Pw.  (for  parabofic  bands  and  displaced  Maxwellian, 
[Py]  =  [Pw]  =  nili£T[I]),  we  use  stationary  Monte  Carlo  calculations  and  physical 
intuition  to  suggest  the  following  distribution  function 


(4) 


as  a  constitutive  relation  to  dose  the  moment  equations.  Here  T^,  has  been  chosen 
to  replace  T  because  jlcgT^  approximates  the  effective  thermal  energy  wdt  for  the 
stationary  transport  case,  and  tn«(7U)  has  been  introduced  because  we  expect  a 
nonconstant  effective  mass  strongly  dependent  on  T^.  Using  Equation  (4)  in  the 
definitions  for  the  transport  parameters,  and  limiting  the  analysis  to  first  and  second 
orders  in  3^.  we  arrive  at  the  the  supplemental  relations:  w  =  | -b  - -  J y .y . 
m»*(T.)  =  m*(l  [Py]  =  ^nJfesT.p].  (Pw)  =  n*j»r.(I]  and 

rj  =  6a(^r)*(fcsl^)*nar.  These  relations  can  be  used  to  dose  Equations  (1)  - 
(3)  expressing  them  in  terms  of  electron  density  n,  average  vdodty  «  and  average 
energy  to. 


3.  The  Quantum  Correction  Terms 


Since  we  desire  a  general  model  suitable  for  studyir<g  ultrasmall  devices,  an  ap¬ 
proach  is  illustrated  to  develop  some  quantum  corrections  for  our  semi-dassicai 
nonparabolic  model  equations.  Grubin  and  Kveskovsky  [2]  have  previously  pre¬ 
sented  a  detailed  derivation  which  will  achieve  th^  goal.  However,  for  our  initial 
investigations  we  use  a  slightly  more  restricted  version  of  their  quantum  hydrody¬ 
namic  equations.  Spedfically,  we  neglect  any  spatial  or  nonparabolic  effects  on  the 
electron  effective  mass  and  ignore  any  defidency  of  the  general  quantum  distribution 
function  to  agree  with  Fermi  statistics  (low  temperature  effects).  We  begin  with 
the  general  moment  equations  of  Strosdo  [5].  The  collisionless  one-dimensional 
forms  of  the  first  three  moments  are: 


Sn 


dpi 

m 


1  s 

nm*  dx 


(n{(p  -  p^))) 


(5) 


(6) 
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^Up  -  pj’))  =  -  p-)*))  -  -  P-)’» 

2  S  ,  „  .  Pd((p  -  Pd)*) 

+  ^ 

where  n  is  the  electron  density,  p  is  the  single  electron  momentum,  is  the 
classical  average  momentum  and  the  total  effective  electric  potential.  If  we 

follow  Ref.  2  and  use  the  momentum  displaced  nonequilibrium  Wigner  distribution 
function, 

of  Ancona  and  lafrate  (6],  we  obtain  n((p  —  p^)*)  =  *^(1  — 

n((p  —  p^)*)  =  0  where  a  =  0  =  5^  and  7  =  These  results  can 

be  used  with  Equations  (5)  -  (7)  to  ^eid  the  collisionless  quantum  hydrodynamic 


equations; 


dn  _ 

~di~~  8x 

1  S  Q.  1 


dw  dta  8  1 8(nvkBT)  .  A*  8  ,1 8n^8v  ,  ^ 

wb«f<  «=  |ipT+  }»•••'-  •"<1  <?= 

ing  the  above  quantum  equations  to  the  semidasskal  nonparatwlk  hydrodynamic 
Equations  (1)  -  (3),  approximating  7^  by  T  and  using  the  supplemental  relations 
of  Section  2  we  arrive  at  the  quantum  corrected  collisionless  nonparaboGc  hydro- 
dynamic  transport  model  in  one  dimension, 

8n  8  . 

^  =  (12) 


vw  /rr  . 

-  +  3  ^  n^T 


1  8(nkBT) 
m*  8z 

A*  8  ,li 


A^  8  I  £»v  ^  p  » 

nm* 


_  8v  ^  ^  (If  2  8 

81  ^  8x  m*u8x  3  Znm*  8x 


2  8  m* 

—  — (Tur(w-Wf-  — 


8w  8v$  8  Q.  2  8  m*  j  8v 

with  nonparaboGc  correction  terms  fi(n,v,w)  =  1-1-  2arw(«,w)(w  —  w,  — 
s/(fs,  V.  w)  =  «(*,  w)  =  An,  V,  w)  =  - 

S^u’)  and  quantum  corrections  to,  =  -  Q  =  -iSF;^VJr- 

Q  b  the  standard  form  of  the  quantum  potential  and  w,  b  a  purely  quantum 
mechanical  contribution  as  first  noted  by  Wigner. 


4.  Results  amd  Conclusions 


An  initial  study  has  been  made  of  a  double-barrier  heterostructure,  Al^^Gai-g^At 
barriers  with  Xc  —  0.3,  doped  to  10*'*cm~*  with  source  and  drain  regions  of  length 
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200  A  and  doped  to  10**cfn“*.  This  structure  is  shown  in  Figure  2  with  the  source 
and  drain  regions  excluded.  In  the  self-consistent  simulations,  the  heterostructure 
was  treated  using  the  Anderson  rule;  =  —q^tSf  — t(*).  where  x(*)  = 

0.697zc)  is  a  position  dependent  electron  affinity  and  the  applied  potential  'itjf 
was  determined  from  solving  Poisson's  equation.  The  dissipation  mechanisms  were 
treated  in  a  classical  manner  [7]  and  we  assumed  TuttUt  =  SOO/T.  The  profile  result 
of  10.  and  io,,  parameters  intrinsic  to  the  energy  balance  equation,  are  displayed 
in  Rgure  1  for  =  8mV.  Figure  2  shows  the  corresponding  profile  result 

for  Q  along  with  the  conduction  band  Ec  (E^  =  This  study  is  the  first 

to  use  the  three  quantum  hydrodynamic  (QHD)  equations  to  study  transport  in 
heterostructure  devices.  The  key  additional  feature  of  the  QHD  equations  is  the 
incorporation  of  the  effects  of  density  gradients  through  the  quantum  potential. 
Insofar  as  S  classical  solution  does  not  exist  in  the  presence  of  barriers,  Q  must  at 
least  cancel  the  effects  of  the  barriers  and  permit  carrier  transport.  Indeed  ^  does 
approximately  balance  the  barrier  potential  a  result  consistent  with  [2]. 


Figure  1.  Average  energy  (io)  and  quantum  correction  (w,)  at  =  8mV. 


mppiitd  =  8mV . 
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Abstract 

Numerical  simulations  of  a  submiaon  GaAs  FET  have  been  perfonned  using 
both  Monte  Carlo  (MQ  methods  and  the  Moments  of  Bwtzmann  Thmqiort 

Equation  (MBT^.  The  I>V  characteristics  as  well  as  details  of  the  internal 
distribotion  of  carriers  and  potential  were  obtained.  The  MC  calculations  show 
no  regions  or  negative  forward  conductance.  However,  the  MBTB  results  show 
that  nerative  forward  conductance  can  be  present  or  al»ent  dq)ending  on  the 
value  of  the  thermal  conductivity. 


Introduction 

Recently,  both  MC  and  MBTE  or  lydrodynamicipCTOadiet  to  device  rimnlatfam 
have  be^  receiving  greater  attention.  Tins  is  a  re^  of  the  need  to  investigate, 
and  accurately  rimulite,  electron  dynamics  in  submicroo  derices.  IXfitfathis 
reduction  in  oevice  size  nonequflilwium  effects  objected  in  the  more 
conventiooal  drift  and  diffiKhwiapptoadi  must  be  considered.  This  rives  impetus 
to  the  use  of  the  more  complex  but  more  ftmdamental  methods  menbooed  above. 

While  attention  has  focused  on  the  inylcmentatioa  of  MC  and  Mii'iEsiinulation 
procedures,  a  direct  conqrarison  of  tb^  methods,  qmlied  to  the  same  problem, 
has  not  been  made.  This  paper  presents  the  results  msudi  a  coimnrative  study. 
While  the  results  are  prelimmary  in  nature,  they  do  provide  signincant  Insight  and 
direction  for  future  work. 


Analysis 

The  device  considered  is  a  sittple  GaAs  FET  with  a  (L6  microo  sooree-drain 
n>acingandaqrnuQetricalIyplaced(lL2tnlcxDOga^  Thederieedqpthistakenas 
(Ll  miaonwitb  auar^fotmo9[dngof  Iz  lO^Vcm^. 

The  MC  procednre  empkiyed  follows  doseljr  that  ootlined  by  Hodmqr  and 
Eastwooo[ll.  The  drange  in  particle  momentum  and  position  are  advanced  in 
time  under  toe  influenoe  of  an  applied  field  assuming  a  non-parrirolk  band 
structure.  At  the  end  of  the  free  tli^peruxl,  the  partides  are  scattered.  This 
process  b  rq>eated  for  all  parddes  nntu  a  qpedfie  tnerement  far  time  is  readied. 
At  that  time  a  qnce  diarge  distribution  b  cooputed  ftmu  the  psrtide  distrfbution 
using  a  doud-mcell  algorithm  and  Poisson's  eqoatioa  b  srived  mxbting  the 
potential  Thb  process  b  repeated  until  a  statistical  steady  state  badiined.  The 
MC  pfXKedure  considers  transfer  the  between  r,  X  sod  L  valleys. 

The  MBTB  approadi  solves  continni^,  momentum  and  energy  equatitms  for  r 
and  Lvall^  electrons.  These  equations  are  expressed  as 


(1)  anj/at  -  -  aniY/ax  -  ni/r^  + 
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(2)  an^Vi/at  -  -  an^Vj^Vj/ax  -  l/m  an^Tik^j/ax  + 
nj^e/aj^  a<^/ax  +  fi  a’V/ax*  -  njy/r2 


(3)  an^Ti/at  -  -  an^ViTi/ax  -  2/3  n^Ti  av/ax  + 


2/3Ko  a/ax  (xaTi/ax)  -  tijT^/Ts  +  njTj/rg  + 


®lVl*/3l«0  ("2/^2  “  "lAl)  +  2niBiViV31Co  *^3 


forr  (subsd^t  »  l)vallqrcaiiieis.  A  similar  ^stem  with  subsist  2  applies  to 
theLvaH^cairieis.  The  equations  are  expressM  here  in  ooe-dimensioiial  form. 
The  extension  the  two-dimeasional  case  considered  here  is  strain  foivwud.  b 
Eg.  (2),  the  mot.^'entum  equation,  the  coeffident  of  the  fourth  tenn  on  the  RHS> 
p,  controls  diss^ntioo.  This  tenn  arise  Grom  the  nonmherical  nature  of  the 
distribution  fiinctioiL  Ingaiera],weassutDethiscoemdenttobezeroorvew 
small.  If  non-gero  and  agnlficant,  the  slope  of  die  1-V  characteristics  at  low  nelds 
win  be  signtficantfy  in  error  and  below  tfairt  anticipated  finom  the  low  field  mobili^ 
of  the  nudetial  in  question.  The  coefficient,  <,mpearing  in  the  third  term  on  the 
RHS  of  the  energy  equation,  is  the  thermal  conmidivi^  of  the  electron  gas. 

While  its  precise  value  is  nnlmown,  it  is  inqrortant  and  Its  role  in  the  belmvior 
the  solntioQS  is  investigated  here.  If  his  zero,  there  wfll  be  no  heat  transfer!^ 
conduction  vASe  if  h  iqiproadies  infini^  the  energy  equation  reduces  to 

(4)  v*T  -  0 


The  quantities  indicated  by  r  in  Eqs.  (1-3)  are  relaxation  times  determined  fiom 
a  prior  evahration  the  coUision  int^rals  ane  are  taken  as  energy 
dqiendent  in  the  present  study  p). 

Equations  (1-3)  and  the  LvaUqr  counter  parts  are  cormled  to  Pohsoo^  equation 
andtheresnlting^stemissolvediisinganalgoti^basedontiieLBI 
(linenhed,  Blodc,  InqiBchltedmiqne  devefaped  by  Brilqr  and  MdXnald  p] 
and  paralldrng  that  nswbyl&eshwm  and  Grubtn  [4]  fm  solving  die  drift  aria 
diffusion  equmioos.  The  ugorithm  is  highly  efficient  and  a  Ugh  degree  of 
parallelism  IS  present  making  h  ideally  smtM  for  vector  processing.  Abrief 
discussion  of  ns  rqndication  to  the  prcmrrt  system  may  be  found  in  [:q. 

Results 


Simulations  of  the  FET  shown  in  Fig.  1  were  performed  for  gate  bias  levels  of  0.1, 
OJ  and  0l5  volts  with  the  drain  bias  varying  from  0  to  2J)  volts.  Both  the  MCand 
the  MBTE  shnulations  used  an  11  z  61  gnd  point  mesh  yidding  square  cells  with 
lOOk^ndr^  Three .'~oosaiid (3000) partk^ (an average ofS per c^ were 
used  in  the  MC  calculations.  The  MBlBdmnlations  requited  qmiosbnately  25 
minutes  of  Gray  XMP  CPU  time  per  biaspdnt  The  Mdsunnlanons  reqiAed 
^jproximatelyOmtnutesperMaspoint  Results  for  MCsunulations  using  9000 

pAidea  showed  no  dgnifiramiKfer»giev*hntfeqirfr«rfrign?firMflygriii^ 

time. 


Figure  2  shows  the  l-Vdmractetrstics  of  the  device  as  determined  fiom  the  two 
simnlationappioadies.  At  a  gate  bias  of  0.1  vdts  his  observed  that  both  the  MC 
and  MBTBr^ts  show  goodagreement  At  this  bias  level,  the  rvaUeyetecOons 
do  not  becwne  too  hot  (maximum  local  termrerature  of  1^0*  la  and  efectioa 
transfetisamaxiinnmOfMproximate^^9arvaIlqfand38%Lva^.  Astbe 
gate  bias  is  increased  to  03  vohs  we  see  that  the  agreement  b  good  up  to 
approximately  0.6  volts  Vjy.  Above  thb  value  of  orain  bias  the  MBTjB  result 
rxniUtsiiMative  forward  oooductanoe.  No  sndi  effect  b  observed  for  tlieMC 
result  Atv^fc  «  2f>  volts,  the  ituiximnmrvaflqr  electron  temperature  fixxm  the 
MBTE  resoureacfaes  2400*  K 1^  the  population  of  caniers  in  the  r  and  L 
valleys  b  20%  and  80%  reflectively  at  the  drain;  the  point  of  maximum  transfer. 
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The  trend  ofwder  disagreement  between  MC  and  MBTBpredictioos  is  scy  to 
grow  as  the  gate  bias  ismereased.  AtV«  »  03  volts  we  observe  that  negative 
ccmductanoe  lypeats  in  ib&  MBTE  simmation  at  about  Vj;  03v^ts.  Here  the 
r  valley  teiiq>etatures  readies  a  maximum  of 292S*  K  ana& the  drain  the 
population  of  carriers  in  the  r  vall^  is  o^  10%  of  the  totaL  AtV^js**  LOvolts, 
the  mairinniin  r  vall»  temperature  is  only  1860*  K,  and  the  V  valley  population  at 
the  drain  is  25%.  It  snould  be  noted  the  population  of  carrier  in  the  Xvall^  from 
the  MC  results  was  ne^gible  at  all  bias  l^els. 

The  reason  for  the  negadve  conductance  is  strai^t  forward.  On«  the  bias  levels 
become  high  enou^  mat  sufficient  energy  is  supplied  to  the  carriers  to  cause 
them  to  transfer,  the  possibility  exists  that  the  carriers  will  transfer  at  a  rate  vidiich 
is  greater  than  the  r  valley  velod^  increases.  Smee  the  current  is  medmmnantly 
through  r  vaUw  transport,  if  this  occurs  the  current  win  decrease.  This  is 
predselyvriiat  B  hajmening  here  and  further  si^ests  two  possible  souroes  of  the 
disctepaacy.  First,  the  scattering  rate  constants.  r{,  which  are  fimcdoiisttf^ 
eleotroa  ten^erature,  could  have  an  iruppropriate  hi^  ten^erature  variatioiL 
However,  uiuer  uniform  field  condidons  the  results  obtained  from  both  MC  and 
MBTE  amulations  yirid  reasonable  velocity  field  carves  wdudi  ate  in  general 

Mtti  <nrwrim/^t«l  Alta  This  suggcsts  that  file  ptobleffl  Ees  elsewiiere, 
and  it  associated  wnth  spatial  nonnnifotmitics  in  the  solatia  Examination  M  the 
contimiity  and  momentum  equatioos  does  not  reveal  aiy  possible  souroes  of  the 
ptoUem.  (As  diseased  eariier;  the  coefficient  of  the  dtsspation  term, /t,  affects 
the  low  field  portion  of  the  I-V  curve,  but  has  little  effect  at  U^biasV  utbe 
energy  equation,  the  onfy  term  vriiidi  can  be  r^atxied  witii  so^pidoa  Is  the  heat 
oonoudioatetnL  This  profited  ns  to  exunine  the  result  if  the  thermal 
oonductiviurwereinfinue.  Inis  result  is  shown  for  Vk  *■  05  volts  as  file  curve 


conductance.  In  fac^  we  have  been  aUe  to  obtain  results  wfakhmatdi  the  MC 
lestfitshyvaningthevi^ofe.  However,  the  requbed  value  of  a  is  hii^  bias 
dqiendent  ForexanH)le,atV«:  -  (LSandV(«s  »  LOavahteofafivetiniratbat 


2J)avahie6f«  100tiaiestbecarve*A*vafaiegaveacaneiitfiiatwasttill20% 
bdow  the  MC  result  One  other  possibility  is  that  the  scattering  rates  should  be 
based  on  a  iiistorkally  averaged  tenperature,  rather  than  the  load  tenpentnre  as 
done  at  present  Oe^,  more  research  needs  to  be  done  here. 


Finally,  we  oonpare  the  distribation  of  potential  and  carrier  den^wifiiin  the 
device,  as  determined  tom  the  two  procedures,  in  F^  3  and  4.  ni^3riiowstbe 
potential  efistribation  at  ^  ■<  03andV,w  LO.  Here  the  I-V  carves  are  far 

good  agreement  and  in  Fig.  3  we  observe  mat  the  potential  distributions  are 
almostidenticaL  The  MC  result  shows  a  sli^  amount  of  noise  which  b  to  be 
eipected  since  the  ^noe  charge  must  be  derived  from  the  discrete  particle 
distribution.  Hbwew,fiie  similarity  in  the  results  ineScatessimilan^  in  the  yace 


simulation.  Here  the  gate  lies  along  fiie  centerline  of  the  figure.  The  similarity  in 
the  dqdetion  re^on  b  erident  here. 


Conclusions 


andMBTBipproadies.  The  two  approaches  were  found  to  give  good  agreement 
in  the  I-V  draracteristics  at  low  values  of  gate  bias.  As  file  g^  bias  vras 
increased,  the  MBTE  rmult  exhibited  negative  forward  conductance  whereas  only 
positive  forward  cooducrance  was  observed  in  the  MC  results.  The  presents  of 
negative  conductance  was  found  to  be  dependent  on  the  magnitude  of  the  thermal 
conductivity  coefficient  in  the  r  valley  energy  equation,  with  negative 
conductance  disappearing  altogether  as  the  thermal  conductivity  ^proacbed 
infinity.  At  bias  leveb  where  good  agreement  was  present  in  the  jn^cted  I-V 
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result,  good  agreement  was  also  observed  in  the  potential  and  carrier  distributln.w 

within  the  dew*.  Overall  the  results  indicate  that  both  MC  and  MBTE 
approaches  are  viable  for  submiaon  device  amnlation,  however,  fiirther  reseaidi 
is  needed  to  provide  proper  dosure  of  tte  MBTE  approad^  ^pc^cally  wift^* 
regard  to  the  role  of  the  neat  flux  terms  m  the  eoer^  equation,  and  the 
evmuadon  of  scattering  rates  based  on  historic  rather  than  local  quantities. 
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ABSTRACT:  Solutions  to  the  first  three  moments  of  the  Boltzmann  transport  equation  and  Poisson's 
equation  are  obtained  for  a  permeable  base  transistor  (PBHr)  using  linearized,  block  implicit  (LBI)  and 
ADI  techniques.  Two  level  electron  transfer  is  considered.  The  results  of  the  simulations  are  compared  to 
results  obtained  from  the  drift  and  diffusion  equations.  The  comparison  indicates  that  nonequilibrium 
transport  and  velodty  overshoot  are  important  in  the  PBT.  The  prediaed  l-V  characteristics  of  the  device 
show  substantially  higher  current  levels  arxf'  a  higher  cutoff  frequency  arc  obtained  with  the  moment 
equations. 


INTRODUCTION 

Numerical  and  experimental  studies  of  the  PBT  were  first  performed  by  Bozler  and 
Alley  [Ij.  Their  results  suggested  that  a  PBT  could  obtain  cutoff  frequendes  as  high 
as  300  GHz.  This  result  has  not  yet  been  obtained.  In  some  of  their  drift  and 
diffusion  computations  Bozler  and  Alley  {IJ  used  a  two  piece  velodty-field  curve  in 
an  attempt  to  account  for  nonequtlibrium  effects.  The  velodty  varied  linearly  with 
field  up  to  a  value  of  4  kV/cm  after  which  it  was  held  constant  at  2  x  lO’cm/sec. 
This  approach  will  yield  higher  current  levels,  but  it  offers  no  improvements  in 
either  the  qualitative  or  quantitative  consequences  of  nonequilibrium  transport. 
Nonequilibrium  effects,  in  the  presence  of  nonuniform  fields,  result  in  carrier  and 
velocity  distributions  that  arc  significantly  different  than  those  offered  by  excess 
saturation  velodty  coupled  to  the  drift  and  diffusion  equations. 

Nonequilibrium  effects  in  a  PBT  have  been  investigated  using  Monte  Carlo 
methods  by  Hwang  et  al.  (2J.  Their  results  showed  significant  differences  compared 
to  the  drift  and  diffusion  results;  the  peak  velocity  along  the  center  of  the  channel 
e.xceeded  4  x  10^  cm/scc  and  the  density  distribution  was  also  significantly  different. 
The  Monte  Carlo  results  also  indicated  a  cutoff  frequency  approximately  60%  higher 
than  the  corresponding  drift  and  diffusion  prediction. 

In  the  present  note,  the  performance  of  the  PBT  structure  shown  in  Fig.  1  is 
analyzed  using  both  the  drift  and  diffusion  and  the  moments  of  the  Boltzmann 
transport  equation.  (MBTE).  The  base  penetration  of  the  present  device  structure  is 
only  one  half  of  the  channel  width.  This  in  itself  results  in  a  higher/,-,  compared  to 
devices  with  equal  base  penetration  and  channel  width  (3J.  However,  further  and 
substantial  increases  in  fy  are  observed  from  the  MBTE  predictions  as  well  as  higher 
current  levels. 
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Figure  I:  Dimensions  and  doping  level  of  the  simulated  PBT. 
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ANALYSIS  AND  NUMERICAL  METHOD 


The  equations  to  be  solved  are  obtained  by  taking  the  first  three  moments  of  the 
Boltzmann  transport  equation.  The  resulting  equations  represent  conservation  of 
mass,  momentum,  and  energy,  and  may  be  expressed  for  the  central  valley  carriers 
as 


any, 

at 


an,  an,V\ 

at  ax,  •  " 

an,vy,  kp  an,T,  _  J_ 
aXf  m,  ax,  m. 


+  «2^12 


axf  m,  ax,  ' 


'r 


(1) 

(2) 


an,T,  an,V\T,  5  ^  _  2  ^  ^  _L 

at  axj  3  '  axf  3  k„  ax,  3k„  ax^ 

—  n,T,r,^  +  n.,T-,r,^  +  (2n,r],  —  njr*,,  +  .  (3) 

A  similar  system  is  obtained  for  the  satellite  valley  carriers.  These  equations  are 
coupled  to  Poisson’s  equation 


cVV  =  e{n,  +  nj  -  N„)  .  (4) 

In  eqs.  (l)-(4)  n  is  carrier  density,  V  the  velocity  vector,  T  the  carrier  temperature, 
(fr  the  electrostatic  potential,  is  Boltzmann's  constant,  m  is  the  effective  electron 
mass.  are  scattering  rate  constants,  cr’’  is  a  stress  tensor  arising  from  the 
nonspherical  nature  of  the  distribution  function,  and  k  is  a  thermal  conductivity.  The 
scattering  rates  are  determined  a  priori  from  evaluation  of  the  collision  integrals  [4]. 

Equations  (l)-(3),  their  satellite  \'allcy  counterparts,  and  Poisson’s  equation  form 
a  coupled  system  of  nine  nonlinear  PDE's  in  two  space  dimensions  of  the  form 
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+  S{t]>) 


(3) 


whore  <f>  represents  the  vector  of  unknowns.  H  and  S  represent  nonlinear  functions 
of  <!>.  and  D{<f>)  represents  a  nonlinear,  partial  differential  operator.  To  solve  this 
s>'stem  Poisson’s  equation  is  decoupled  from  the  remaining  eight  equations  by 
differencing  the  electric  field  appearing  in  the  momentum  equation  at  the  explicit 
time  level.  The  remaining  equations  are  solved  by  direct  application  of  linearized 
block  implicit  (LBI)  techniques  (5).  The  equations  arc  first  time  differenced 


At 


=piDi4,r"  +  si4>r"]+(\-fi)[D{4>y+s{4>r]. 


(6) 


The  nonlinear  operators  are  then  linearized  using  Taylor  series,  as  for  example; 

Dir*')  =  Dir)  +  ^r^'  +  OiAr)  .  (7) 

dtp  1 
I 

Equation  (6)  may  then  be  expressed  as 


iA  +  dtL)ar*' =^MDir)+sir)]  (s) 


where 


Equation  (8)  is  then  split  following  the  ADI  procedure  of  Douglas  and  Gunn  |6). 

iA  + AiLJA<b*  =  At[Dir)  + Sir)]  (10a) 

(A  +  AtLJdr*  =  A^r  i  10b) 

and 

Ar"  =  Ar* +  OiAr).  (11) 

.As  a  result  of  the  ADI  splitting,  the  number  of  operations  required  to  solve  the 
remaining  coupled  equations  varies  linearly  with  the  total  number  of  grid  points,  and 
due  to  the  linearization  process,  the  solution  can  be  advanced  in  time  w'ithout 
introducing  nonlinear  iteration.  After  solution  of  the  continuity,  momentum,  and 
energy  equations.  Poisson's  equation  is  solved  using  a  scalar  ADI  procedure  to 
complete  a  time  step.  Steady  solutions,  such  as  those  to  be  presented  here,  are 
obtained  from  the  long  time  asymptotic  transient  solution.  Since  transient  accuracy 
is  not  of  interest  under  such  conditions,  the  time  step  can  be  spatially  scaled  to  speed 
convergence  to  steady  state. 

The  solutions  to  the  drift  and  diffusion  equations  reported  here  were  obtained 
using  a  related  procedure,  described  in  detail  in  (7). 
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COMPUTKI)  RKSUI.TS 

The  computed  l-V  characteristics  of  the  device  shown  in  Rg.  1  arc  presented  in  Fig. 
2.  The  results  for  the  MBTE  calculations  were  extrapolated  to  the  origin,  as 
indicated  by  the  long  broken  lines.  The  shorter  dashed  curves  show  the  results  for 
the  DDE  (.II.  The  comparison  shows  that  the  predicted  current  levels  arc  significant¬ 
ly  higher  for  the  MBTE  solutions,  a  result  consistent  with  FET  calculations 
performed  by  the  present  authors,  as  well  as  by  Cook  and  Frey  |8J  who  used  a  highly 
simplified  momentum-energy  transport  model.  The  present  MBTE  results  also 
indicate  a  rCgion  of  negative  differential  forward  conductivity  at  =  0.6  V.  The 
origin  of  this  phenomena  is  believed  to  be  a  consequence  of  electron  transfer.  The 
presence  of  a  dc  negative  forward  conductance  is  also  a  feature  of  PBT  measure¬ 
ments  (9],  but  is  clearly  absent  from  DDE  simulations. 

A  companson  between  the  total  carrier  density  distribution  along  the  center  of  the 
channel  for  drift  and  diffusion  and  MBTE  solutions  is  shown  in  Hg.  3  for 
V(.fr  =  1.0 V  and  Vf,^  =  QAV.  The  MBTE  prediction  for  the  gamma  valley  carrier 
density  is  also  shown.  The  drift  and^diffusion  result  was  obtained  using  an  equilib¬ 
rium  velocity-field  curve,  and  yields  results  qualitatively  similar  to  that  of  Bozicr  and 
Alley. 

As  seen  in  Rg.  3  for  the  DDE  simulations,  the  carrier  density  reaches  a  maximum 
between  base  contacts.  Here,  with  the  velocity  in  saturation  and  the  cross  sectional 
area  at  a  minimum,  the  carrier  density  must  increase  to  maintain  current  continuity. 
In  the  MBTE  simulation  the  constraints  of  current  continuity  are  more  complex. 
Rrst  a  decrease  in  the  cross  sectional  areas  is,  as  in  the  DDE,  accompanied  by  an 
increase  in  field  along  the  channel.  The  field  increase  under  both  equilibrium  and 
nonequilibrium  conditions  is  qualitatively  similar,  as  may  be  observed  from  Hg.  4 
which  shows  the  potential  distribution  along  the  center  of  the  PBT  channel. 


Figure  2:  Colicciar  currcnl  vs.  collector  emitter  voltage  for  different  values  of  hase  etniiicr  xohaje. 
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However,  consequent  changes  in  electron  temperature,  both  increasing  and  decreas- 
in».  lag  behind  the  equilibrium  state.  This  leads  to  velocity  overshoot  and  a  delay  in 
electron  transfer.  As  a  result,  for  pearly  the  first  half  of  the  device  transport  is 
almost  exclusively  gamma  valley  transport.  The  implication  is  that  if  the  gamma 
valley  carrier  velocity  increases  with  increasing  field,  then  the  produa  of  density  and 
cross  sectional  area  normal  to  current  flow  must  decrease  to  maintain  current 
continuity.  Since  the  velocity  increases  faster  than  the  area  decrease  the  carrier 

density  decreases.  .  c  .j  u 

At  moderate  bias  levels  typical  FET  calculations  show  a  decreasing  field  as  the 
gate  region  is  passed.  This  also  occurs  in  the  PBT.  Now,  as  the  cross  sectional  area 
increase  the  gamma  valley  carriers  exhibit  a  decrease  in  velocity.  It  must  be  noted, 
however,  that  for  the  parameters  of  the  calculations  the  L  valley  carriers  make  a 
negligible  contribution  to  current.  Thus  a  decrease  in  carrier  velocity  results  in  a  net 


Figure  4:  Polcnlial  vs.  disiance  along  center  of  channel  for  the  PBT. 
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Figure  5:  Cutoff  frequcnqt  vs.  base  emitter  wltagc. 


increase  in  carrier  concentration.  However  initially,  the  decrease  in  field  is  not 
accompanied  by  a  corresponding  temp>erature  decrease  (as  experienced  in  the 
uniform  field  calculations).  Thus,  the  high  gamma  valley  temperature  results  in 
transfer  to  the  L  valley  giving  rise  to  the  second  minimum  in  the  gamma  valley 
carrier  densitv  ..hown  in  Fig.  3.  Further  toward  the-  drain,  the  field  decreases. 
However,  relaxation  is  incomplete  and  the  field  at  the  collector  is  not  equal  to  the 
field  at  the  emitter.  Also  note,  ttu.  MBTE  potential  distribution  gives  rise  to  a 
slightly  higher  field  upstream  of  the  base,  and  a  lower  field,  over  a  longer  distance, 
downstream  of  the  base  compared  to  the  DDE  result.  More  significantly  the  electron 
temperature  at  the  collector  exceeds  that  at  the  emitter.  It  is  noted  that  as  the  field 
rela.xes.  the  electrons  transfer  back  to  the  central  valley. 

All  of  the  results  presented  thus  far  are  qualitatively  similar  to  the  Monte  Carlo 
results  of  Hwang  et  al.  {2j.  However,  the  results  for  cutoff  frequency,  sho.vn  in  Fig. 
5.  are  qualitatively  different.  At  V'j.^  =  0.5V,  is  seen  to  incriase  with  from 
approximately  100  GHz  to  130  GHz.  value  2.5  io  3  times  greater  than  that  obtained 
using  the  DDE.  However,  for  higher  values  of  fj.  is  shown  to  decrease  with 
increasing  For  V^£  =  1.5V,  an  in  e.xcess  of  200  GHz  is  predicted  for 

+  0.2  V;  while  at  =  0.6  V,  has  dropped  below  40  GHz.  a  result  that  is 
qualitatively  consistent  with  the  presence  of  the  negative  forward  conductance,  and 
one  that  would  find  no  explanation  in  the  DDE  cal'-”-  tions. 
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Electron  velocity  overshoot  and  valley  repopulation  effects  in  diamond 
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Using  ensemble  Monte  Carlo  procedures,  electron  velocity  overshoot  and  transient  valley 
repopulation  in  diamond  have  been  investigated  as  a  function  of  field  strength  and  orientation. 

It  is  foimd  that  the  response  of  electrons  to  the  sudden  application  of  a  50  ItV/ cm  electric  field 
along  <  100>  results  in  a  velocity  transient  with  a  maximum  of  1.9X  10’  cm/s  after  200  fs  and  a 
steady  state  value  of  1.30X  10’  cm/s.  For  a  field  along  (110),  the  corresponding  maximum  and 
steady  state  values  are  2.0X  10’  and  1.4X  10’  cm/s,  respectively.  The  calculated  temporal  and 
spatial  duration  of  velocity  overshoot  in  diamond  is  longer  than  that  of  silicon  but  shorter  than 
that  of  gallium  arsenide. 


The  superior  thermal  conducting  characteristics  of 
semiconducting  diamond,  coupled  to  hole  and  electron  mO' 
bilities  and  saturated  drift  velocities  that  are  superior  to  sili' 
con,  make  it  an  ideal  candidate  for  high-speed,  high-power, 
and  radiation-resistant  device  applications. To  exploit 
this,  much  effort  is  under  way  in  the  growth  and  doping  of 
synthetic  diamond,  while  interest  in  natural  diamond  re¬ 
mains  high.*"’  The  successful  development  of  a  semicon¬ 
ducting  diamond  device  technology  requires  the  determina¬ 
tion  of  a  number  of  key  transport  parameters,  such  as 
mobilities,  transient  velocities,  etc.  The  mobility  and  steady 
state  drift  velocity  of  electrons  and  holes  in  natural  diamond 
have  been  investigated  by  Nava  et  al,^  and  Reggiani  et  al.^  In 
this  letter  we  examine  two  effects  occurring  in  diamond  for 
thefirsttime:  (1)  velocity  overshoot  and  (2)  transient  valley 
repopulation  (as  have  been  investigated  in  Si  and  Ge*).  The 
calculations  reported  below  are  for  electrons  subjected  to 
sudden  changes  of  field  with  values  20  and  50  kV/cm.  Addi¬ 
tionally,  for  fields  of  50  kV/cm,  the  orientation  dependence 
is  studied.  It  is  found  that  velocity  overshoot  occurs  on  tem¬ 
poral  scales  of  the  order  of 400  f^  which  is  longer  than  that  of 
silicon  and  the  same  order  as  that  of  gallium  arsenide.  The 
spatial  duration  of  the  overshoot  is  of  the  order  of 400  A,  less 
than  that  of  galltiim  arsenide’  but  larger  than  that  of  rili- 
con. Thus  nonequilibrium  effects  arc  expected  to  be  impor¬ 
tant  at  featiue  sizes  below  400  A. 

In  the  ensemble  Monte  Carlo  procedure  used  in  this 
study  of  diamond,  the  kinetics  of  electrons  are  determined 
relative  to  the  minimum  in  k  space  at  which  the  valley  is 
located.  Electrons  located  in  valleys  on  opposite  sides  of  the 
same  axis,  e.g.,  the  ( KW)  and  ( 100)  valleys,  are  not  distin¬ 
guishable,  and  symmetry  allows  specific  calculations  to  be 
performed  for  three  ellipsoidal  valleys  (m,  =  1.4,  m, 
—  0.36)  located  along  the  (100),  (010),  and  <(X)1)  crystal¬ 
lographic  orientations.  The  Monte  Carlo  procedure  takes 
into  account  scattering  by  elastic  intravalley  acoustic  phon¬ 
ons,  g-type  intcrvaJley  phonon  scattering  between  two  paral¬ 
lel  valleys,  and /type  intervalley  phonon  scattering  between 
four  perpendicular  valleys.  Following  an  /type  scattering, 
the  final  valley  is  selected  from  the  possible  perpendicular 
valleys  using  a  random  number.  The  interaction  between  the 
electrons  and  ionized  impurities  is  ignored  in  this  investiga¬ 
tion.  The  longitudinal-optical  phonon  modes  assist  interval¬ 
ley  g  scattering  and  have  a  temperature  of  1900  K,  while 


longitudinal-acoustic  and  transverse-optical  phonon  modes 
assist  intervalley / scattering  and  have  temperatures  of  1500 
and  1720  K,  respectively.*  These  energies  of  the  intervaHey 
phonons  are  considerably  larger  than  the  corresponding  val¬ 
ues  in  SL*  A  single  value,  8X10*  eV/cm,  was  iised  for  all  of 
the  intervalley  deformation  potentials,  and  a  value  of  8.7  eV 
was  used  for  the  intravalle^'  acoustic  deformation  potential.^ 
Initially,  the  electrons  were  assumed  to  be  in  thermal  equi¬ 
librium  corresponding  to  a  Maxwellian  distribution  at  room 
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HG.  I.  Drift  velocity  of  electrons  in  diamond  as  ■  function  of  (a)  time  and 
(b)  position  at  300  K.  The  orientations  of  the  electric  field  are  shown  for 
each  curve. 
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temperature.  The  electric  field  was  then  switched  on  and  the 
time  evolution  of  the  electron  ensemble  was  examined  for 
fieldsof  20  and  SO  kV/cm  applied  along  the  (100)  direction. 
To  illustrate  the  orientation  dependence  of  valley  repopula¬ 
tion  and  velocity  overshoot,  the  calculations  were  repeated 
for  a  50  kV/cm  electric  field  along  the  ( 1 10)  direction. 

The  time  evolution  of  the  electron  velocity  during  the 
first  picosecond  after  switching  on  the  electric  field  is  shown 
in  Fig.  1.  For  a  50  kV/cm  electric  field  along  the  (100) 
direction,  the  transient  velocity  reaches  a  maximum  of 
1.9x  10^  cm/s  after  200  fs  and  then  decreases  to  a  steady 
state  value  of  1.30X  10’  cm/s.  When  a  20  kV/cm  field  is 
applied  along  (100),  the  velocity  increases  more  graduaUy 
to  a  maximum  of  1.4  X 10’  cm/s  after  300  fs,  and  then  de¬ 
creases  to  a  steady  value  of  1.0  X 10’  cm/s.  Changing  the 
orientation  of  the  electric  field  at  50  kV/cm  to  (110),  in¬ 
creases  the  maximum  transient  velocity  and  steady  state  ve¬ 
locity  to  2.00x10’  and  1.4x10’  cm/s,  respectively.  The 
steady  state  results  are  consistent  with  the  earlier  experimen¬ 
tal  work  reported  by  Nava  et  al.,  which  showed  higher  drift 
velocities  in  natural  diamond  for  fields  along  (110)  com¬ 
pared  to  that  along  ( 100).’  In  all  of  the  above  situations,  the 
maximum  velocity  is  about  50%  higher  than  the  corre¬ 
sponding  steady  state  velocity.  The  magnitudes  of  the  maxi¬ 
mum  and  the  steady  state  velocities  are  larger  than  those 
reported  for  Si  at  50  kV/cm.*  Additionally,  the  time  dura¬ 
tion  of  velocity  overshoot  is  longer  in  diamond  than  in  Si. 
For  example,  for  a  20  kV/cm  field  along  (100),  the  time 
duration  of  overshoot  is  approximately  250  fs  in  Si*  com¬ 
pared  to  400  fs  in  diamond  as  can  be  seen  from  Fig.  1(a).  The 
longer  overshoot  duration  makos  it  possible  to  design  diamond 
devices  with  larger  dimensions  than  Si  devices,  while  main¬ 
taining  shorter  transit  times  across  the  device.  Figure  1  (b) 
shows  that  the  distance  over  which  overshoot  effects  are  pro¬ 
nounced  is  approximately  400  A  for  diamond.  This  over¬ 
shoot  distance  is  longer  than  the  300  A  in  silicon  at  compara¬ 
ble  field  values.'” 

Figure  2  displays  the  time  evolution  of  the  average  ener¬ 
gy  of  the  carriers  subject  to  fields  of  20  and  50  kV/cm.  For 
the  50  kV/cm  field,  and  during  the  first  200  fs,  the  average 
energy  of  the  carriers  rises  rapidly  and  is  nearly  independent 
of  orientation.  However,  at  longer  times  the  average  energy 
for  the  field  along  ( 1 10)  is  slightly  higher  than  that  along  the 
(100)  direction.  In  response  to  a  20  kV/cm  field  along 
(100),  the  average  energy  increases  gradually  during  the 
first  400  fs  and  reaches  a  value  of  130  mV  which  is  below  the 
threshold  for  the  emission  ofg-type  intervalley  phonons.  To 
illustrate  the  valley  repopulation  effects  and  how  it  in¬ 
fluences  the  velocity  overshoot,  the  mean  energy  and  the 
number  of  electrons  in  individual  valleys  were  determined. 
In  Fig.  3,  the  average  energy  of  the  electrons  in  the  ( 100) 
(dashed  lines)  and  (001 )  valleys  (solid  lines)  is  plotted  for 
20  and  50  kV/cm  electric  fields.  Note  that,  regardless  of  field 
orientation,  e.g.  along  ( 100)  or  (110),  the  component  of  the 
electric  field  along  the  principal  axis  of  the  (001 )  valleys  is 
zero,  while  it  is  finite  along  the  minor  axes  where  transport  is 
determined  by  the  smaller  transverse  effective  mass.  Conse¬ 
quently,  the  electric  field  heats  the  electrons  in  the  (001) 
valleys  at  a  rate  faster  than  that  of  the  ( 100)  valleys.  This 


FIG.  2.  Mean  energy  as  a  function  of  date  forelectroos  in  diamond.  Differ¬ 
ent  curves  refer  to  different  magnitudes  and  orientadoo  of  the  dectric  field. 

results  in  an  initial  rapid  rise  in  average  energy,  as  well  as  a 
higher  average  energy  at  longer  times  for  the  electrons  in  the 
((X)l)  valleys  relative  to  those  in  the  (100)  valleys  (seeing.  3 
for  both  20  and  50  kV/cm  fields). 

The  (001)  and  (010)  valleys  are  equivalent  when  the 
field  is  along  (1(X))  directions.  Hence  the  rapid  heating  of 
electrons  in  these  valleys  results  in  a  significant  fiactimi  of 
the  electrons  undergoing /^type  intervalley  phonon  scatter¬ 
ing  to  the  (  100)  valleys,  whose  transport  is  governed  by  the 
larger  longitudinal  effective  mass. ' '  This  reduces  the  popu¬ 
lation  of  the  electrons  in  the  (001)  vaUeys  while  increasing 
that  of  ( 100)  valleys,  as  can  be  seen  in  I^.  4.  For  example,  in 
response  to  a  50  kV/cm  field  along  (1(X))  [I^  4(a)],  the 
population  of  the  electrons  in  the  ( 100)  valleys  rapidly  in- 
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FIG.  3.  Mean  energy  as  a  function  of  dme  for  dectron  in  diamond  in  (001 ) 
valleys  (solid  lines)  and  ( 100)  valleys  (dashed  lines)  at  300 1C. 
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FIG.  4.  Fraction  of  electrons  in  (100)  valleys  (dashed  lines)  and  (001) 
valleys  (solid  lines)  as  a  function  time  after  switching  on  the  held  at  300  K: 
(a)  /•=50kV/cm  ||(100),  (b)  /•=  50  kV/cm,  |1(110).  and  (c)  £=60 
kV/cm  11(100). 

creases  from  its  equilibrium  value  to  50%  of  the  total  after 
3(X)  fs  and  then  gradually  decreases  to  45%  after 700  fs.  The 
physics  behind  this  repopulation  is  displayed  in  Rg.  3,  which 
shows  that  after  200  fs  the  average  energy  of  the  electrons  in 
the  (100)  valleys  exceeds  the  threshold  for  intervalley 
phonon  emission.  This  leads  to  a  larger  fraction  of  /^type 
intervalley  scattering,  a  reduction  of  the  ( 1(X))  valley  popu¬ 
lation,  and  an  increase  in  (001 )  valley  population. 

When  the  orientation  of  the  field  is  changed  to  { 1 10), 
the  ( 100)  and  (010)  valleys  are  equivalent;  and  transport  in 
these  valleys  is  through  a  combination  of  longitudinal  and 
transverse  efiective  masses.  Several  new  features  enter:  First, 
the  energy  difference  between  the  (001 )  and  ( 100)  valleys  is 
smaller  than  for  the  field  along  ( 100);  second,  because  of  the 
contribution  of  the  transverse  effective  mass,  the  electrons  in 
the  ( 100)  valleys  have  a  higher  average  energy  and  heat  up 
at  a  faster  rate  compared  to  the  situation  when  the  field  is 
along  the  <  100)  axis  as  can  be  seen  in  Fig.  3.  Consequently, 
the  difference  between  the  /-type  intervallcy  scattering  from 
(KX))  to  (001),  and  the  intervalley  scattering  rate  from 
(001 )  to  ( 100)  are  smaller  for  the  <  1 10)  orientation  than  for 


the  ( 100)  orientation,  and  the  population  of  the  electrons  in 
both  (001 )  and  ( 100)  valleys  changes  by  less  than  5%  from 
the  initial  thermal  equilibrium  values.  The  steady  state  ener¬ 
gy  of  the  electrons  in  the  ( 100)  valleys  remains  well  below 
the  f-typc  phonon  emission  threshold  (sec  Rg.  3)  when  a  20 
k  V/cm  is  applied  along  the  ( 100)  axis,  so  that  only  a  smaller 
fraction  of  electrons  undergoes  intervalley  scattering  to 
(010)  and  (001)  valleys.  Here  the  population  of  the  elec¬ 
trons  in  the  ( 100)  valley  gradually  increases  to  50%  after  1 
ps,  while  that  of  the  ((X)l )  valleys  decreases  to  25%  of  the 
total  electron  population  [Fig.  4(c)  ]. 

In  conclusion,  we  have  shown  that  the  electrons  in  dia¬ 
mond  exhibit  orientation-dependent  velocity  overshoot  and 
that  the  temporal  and  spatial  duration  of  the  overshoot  is 
longer  compared  to  Si.  The  maximum  velocities  were  about 
50%  larger  than  the  steady  state  velocities.  The  magnitudes 
of  the  maximum  velocities  following  a  50  kV/cm  field  were 
2.00X10’  cm/s  for  (110)  field  orientation  and  1.9x10’ 
cm/s  for  (100)  orientation.  The  valley  repopuladon  is 
strongly  affected  by  the  orientation  and  the  magnitude  of  the 
electric  fields.  Finally,  the  influence  of  the  ionized  impurity 
scattering  on  the  above  conclusions  remains  to  be  investigat¬ 
ed.  Because  the  known  donor  impurities,  such  as  activated 
lithium,  require  more  than  100  meV  to  ionize  and  release 
free  electrons,  the  concentration  of  free  electrons  and  ionized 
donors  is  strongly  dependent  on  the  electric  field.  This  can 
influence  the  transient  response  of  electrons  to  applied  elec¬ 
tric  field  due  to  the  finite  time  required  to  release  the  elec¬ 
trons  from  the  donor  impurities. 
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Abstract 


The  role  of  the  electron-electron  (e-e) ,  hole-hole  (h-h) ,  and  electron-hole  (e-h) 
interaction  on  ultrafast  cooling  of  carriers  in  GadVs  is  examined  for  excess  excitation 
energies  of  40,  200,  and  300  meV  using  em  Ensemble  Monte  Carlo  (EMC)  approach.  It  is 
found  that  when  the  initial  energy  of  the  czurrier  is  below  the  phonon  emission  threshold 
carrier-carrier  (c-c)  interactions  stimulate  either  the  optical  phonon  emission  or 
absorption  process  depending  on  whether  the  initial  energy  of  the  carrier  is  edjove  or 
below  the  thermal  energy,  respectively.  The  e-h  interaction  role  is  strong  when 
excitation  energy  is  below  the  LO  phonon  emission  threshold. 


Introduction 


In  recent  years,  the  time-resolved  femtosecond  spectroscopy  and  optical  time  of 
flight  measurement  techniques  have  become  an  attractive  tool  for  investigating  the 
ultrafast  dynamics  of  photoexcited  carriers  in  semiconductors^"^ .  Consideraible 
experimental  information  about  the  carrier-carrier  and  carrier-phonon  interactions 
processes  has  been  obtained.  These  processes  show  a  strong  dependence  on  carrier 
concentration,  energy,  electric  field,  temperature,  and  material  properties.  For 
example,  at  high  electric  fields  the  minority  electrons  in  p-type  InGaAs  exhibit  higher 
saturation  velocities  than  majority  carrier  electrons^,  while  in  p-type  Si  they 
exhibit  mobilities  that  are  much  lower  than  majority  carrier  electrons®,  similarly  an 
increase  in  the  electron  energy-loss-rates  in  the  presence  of  a  cold  hole  plasma®,  and 
an  absolute  negative  mobility  for  the  minority  electrons  and  holes  at  low  temperatures 
in  GaAs/AlGaAs  quantum  well  structures  have  been  observed^,  in  all  these  situations, 
the  e-h  interaction  has  been  suggested  as  the  possible  responsible  mechemism.  More 
recently  it  has  also  been  suggested  that  electrons  excited  by  excess  energy  of  0.9  eV 
relax  initially  by  emitting  LO  phonons  and  subsequently  thermalize  through  e-e 
interactions’^.  These  remarkable  experimental  results  clearly  demonstrate  our  limited 
understanding  of  the  manner  in  which  the  carrier-carrier  (e-e,  h-h,  e-h)  interactions 
influence  the  cooling  and  response  to  externally  applied  electric  fields  of  the 
photoexcited  carriers  at  different  excitation  energies  and  levels. 

In  this  paper,  we  examine  the  energy  dependence  of  the  role  of  e-e,  h-h,  e-h,  and 
carrier-phonon  interactions  on  the  cooling  rates  of  electrons  and  holes.  The  EMC 
approach  is  used  to  avoid  making  any  assumptions  on  the  form  of  the  distribution 
functions  for  the  holes  and  electrons  which  has  been  a  standard  practice  in  previous 
theoretical  attempts  to  investigate  the  cooling  process®' 9.  With  its  built-in  time 
evolving  distribution,  the  EMC  approach  makes  it  possible  to  exeunine  the  details  of  the 
cooling  process  on  picosecond  emd  sub-picosecond  time  scales.  In  an  earlier  study,  we 
examined  the  concentration  dependence  of  em  e-h  plasma  excited  by  a  200  femtosecond, 
1.64  ev  laser  pulse^O.  It  was  found  that  the  e-h  interaction  provided  a  significemt 
channel  through  which  the  electron  energy  is  transferred  to  the  lattice  pai-ticularly  at 
high  excitation  levels  where  the  electron-phonon  interaction  is  strongly  screened  and 
the  e-h  scattering  events  are  more  frequent.  The  cooling  rates  were  lower  than  those 
predicted  by  a  model  that  assumed  a  parameterized  distribution  function  with  equal 
electron  and  hole  temperatures  and  includes  the  hot  phonon  effects.  More,  recently  we 
also  examined  the  hot  phonon  and  e-h  interaction  using  EMC  and  showed  that  for  times 
less  than  2  ps  the  e-h  interaction  controlled  the  relaxation  process  while  for  times 
greater  than  2  ps  the  not  phonons  slowed  the  cooling  rates^^. 


In  the  present  study,  situations  where  the  e-h  plasma  is  photoexcited  with  excess 
energies  of  40,  200,  and  300  meV  are  investigated.  These  energies  correspond  to 
situations  where  the  electron  energy  (Eg) ,  the  hole  energy  (E^) ,  and  the  LO  phonon 
energy  (Eg)  satisfy  the  following  conditions: 


(1)  Ee  <  Eo  and  Eh  <  Eq, 

(2)  Ee  >  Eq  and  Eh  <  Eq, 

(3)  Ee  >  Eo  and  Eh  >  Eq, 

respectively.  Hot  phonon  effects  are  not  considered.  The  density  of  the  e-h  plasma  is 

assumed  to  be  5  x  lO^^cm"^,  and  the  finite  duration  of  the  pulse  is  ignored.  From 
the  results  of  this  study,  conclusions  are  drawn  about  the  influence  of  the  initial 
energy  distribution  of  the  photoexcited  carriers  on  the  manner  in  which  c-c  interactions 
achieve  thermalization  of  the  carrier  distributions. 

Theory 

The  EMC  proceedure  used  in  this  study  is  the  same  used  in  investigating  the  role  of 
e-h  interaction^^' ^2 ^  Thus  only  its  main  features  will  be  summarized  for  the  sake  of 
completeness.  The  band  model  consists  of  three  r.onparabolic  valleys  for  the  conduction 
band  and  a  parabolic  heavy  hole  valence  band.  The  light  hole  band  is  ignored,  partly 
because  of  its  small  density  of  states  and  to  reduce  the  computation  time.  Figure  1, 
shows  the  flow  chart  of  the  EMC  program  used  in  this  investigation.  The  first  step 
corresponds  to  specifying  the  material  parameters,  initial  energy  distribution  of  the 
carriers,  pulse  shape  and  duration,  room  temperature,  value  of  the  electric  field,  and 
which  c-c  interactions  are  active.  The  scattering  rates  for  acoustic  phonons, 
deformation  potential,  intervalley  and  intravalley  scattering  processes  for  electrons 
and  holes  are  then  calculated  and  tabulated.  On  the  other  hand,  the  scattering  rates  of 
c-c  and  screened  carrier-phonon  interactions,  which  depend  on  the  evolving  distribution 
function  of  the  carriers  are  calculated  every  50  femtosecond  initially  and  at  longer 
time  intervals  as  the  time  proceeds.  For  finite  pulse  duration,  the  number  of  the 
carriers  is  adjusted  according  to  the  pulse  shape  with  an  initial  energy  determined  from 
the  energy  of  the  pulse  taking  into  account  the  nonparabolicity  of  conduction  band.  The 
dynamics  of  the  electrons  and  holes  are  then  examined  using  the  regular  EMC  approach. 
Coupling  between  the  electron  EMC  and  the  hole  EMC  programs  is  accomplished  through  the 
e-h  interaction  which  depends  on  the  energy  and  momentum  distribution  of  the  electrons 
and  holes.  The  average  energy  of  the  electrons  and  holes,  energy-loss-rates  through  e-h 
and  electron-phonon  interactions  are  then  calculated  at  the  end  of  each  iteration  or  at 
specified  time  intervals. 

The  electrons  and  holes  are  assumed  to  interact  among  themselves  and  with  each  other 
through  an  exponentially  screened  Coulomb  potential.  The  screening  length  is  calculated 
self  consistently  in  the  static  and  long  wavelength  limit  of  the  random  phase 
approximation.  Unlike  previous  models  for  e-h  scattering,  no  a  priori  assumptions  are 
made  on  the  magnitude  of  the  energy  exchange  between  the  electrons  and  holes  in  an  e-h 
scattering  event.  The  expressions  for  the  scattering  rates  for  c-c  and  screened 
carrier-phonon  interactions  are  given  in  reference  12.  Figure  2,  shows  the  scattering 
rates  for  e-e,  e-h,  h-h  interaction  processes,  assuming  a  Maxwellian  distribution  for 
the  electrons  and  holes  at  temperatures  of  1000  K  and  160  K  respectively  as  a  function 
o..'  the  average  wavevector  kz.  From  this  figure,  it  is  clear  that  the  h-h  scattering 
rates  are  much  stronger  than  the  e-e  rates  due  to  the  larger  density  of  states  in  the 
heavy  hole  band.  Also  notice  that  the  interaction  of  an  electron  with  a  hole  plasma  is 
stronger  than  the  rate  at  which  a  hole  interacts  with  a  sea  of  electrons^^.  This  is 
related  to  the  large  density  of  final  states  for  holes  compared  to  that  of  electrons. 

The  screening  of  the  carrier-phonon  interactions  reduces  the  scattering  rates 
significantly,  especially  at  high  carrier  concentrations  and  low  temperatures. 
Additionally,  the  presence  of  the  relatively  cold  hole  plasma,  which  is  the  case  in  a 
photoexcitation  process,  leads  to  stronger  screening  even  when  the  energy  of  the 
electrons  is  sufficiently  high^2.  m  the  present  analysis  the  acoustic  phonon 
scattering  is  assumed  to  be  elastic  which  is  a  good  approximation  at  lattice  temperature 
of  77  K.  However,  for  excitation  by  total  excess  energy  of  40  meV,  the  use  of  an 
inelastic  acoustic  scattering  model  is  more  appropriate  because  the  relaxation  through 
the  optical  phonon  emission  is  severely  diminished.  Furthermore,  the  interaction 
between  the  L-valley  electrons  and  holes  is  neglected,  because  the  population  of  the 
electrons  the  upper  valleys  is  negligible.  The  elastic  impurity  scattering  is  ignored 
in  our  calculation,  since  we  assume  that  all  of  the  carriers  are  generated  by  the 
photoexcitation  process.  In  order  to  examine  the  roles  of  e-e  ,  h-h,  and  e-h 
interactions  more  closely,  the  phonon  system  was  assumed  to  be  in  eguilibrium. 


Results  and  discussion 


The  cooling  rates  of  carriers  photoexcited  with  total  excess  energies  of  40,  200  and 
300  meV,  corresponding  to  photon  energies  of  1.55,  1.71,  and  1.81  eV  respectively,  were 
examined  under  four  different  combinations  of  active  scattering  processes.  First  the 
simulation  was  carried  out  assuming  only  optical  phonon  interactions  were  active.  Curve 
1  in  Figures  3  through  5  corresponds  to  this  situation.  Notice  that  curve  1  in  Figures 
3  and  4b,  which  represents  cooling  of  electrons  and  holes  with  initial  energies  below 
the  phonon  emission  threshold  Eg  (36.4  meV) ,  does  not  change  with  time  since  no  phonon 
emission  is  possible  and  at  77  K  and  the  probaibility  for  optical  phonon  absorption  is 
small.  However,  when  the  initial  energy  of  the  electrons  is  above  Eg  (Figures  2a  and 
3a) ,  the  electrons  emit  integral  multiple  of  LO  phonons  and  end  up  with  energies  equal 
tc  27.9  meV  and  4.25  meV  at  4  picosecond  after  excitation  by  1.71  and  1.81  eV  laser 
pulses,  respectively.  Similarly  the  holes  reach  a  steady  state  energy  of  5.8  meV  which 
corresponds  to  emitting  a  single  optical  phonon  (see  Figure  3b).  After  reaching  these 
energies  no  further  cooling  is  possible  to  reach  equilibrium  with  the  lattice,  since  the 
acoustic  scattering  is  assumed  to  be  elastic. 


The  situation  where  the  e-e  and  h-h  interactions  were  active,  in  addition  to  the 
optical  and  acoustic  phonons,  was  examined  next.  Curve  2  in  these  figures  shows  the 
cooling  process  for  the  electrons  and  holes.  The  most  interesting  phenomena  in  Fig.  1 
is  that,  the  e-e  interaction  enhances  the  cooling  of  the  electron  system  from  an  initial 
energy  of  35  meV,  which  is  just  below  the  phonon  emission  threshold,  to  a  steady  state 
energy  13.5  meV  that  is  slightly  above  the  room  temperature.  This  enhancement  occurs 
because  some  electrons  gain  enough  energy  from  interacting  with  other  electrons  and  then 
cool  down  by  emitting  LO  phonons.  On  the  other  hand,  the  h-h  interaction  coupled  with 
optical  phonon  absorption  drives  the  hole  system  to  thermal  equilibrium  with  the 
lattice.  The  average  energy  per  hole  changes  slowly  from  an  initial  energy  of  5  meV  to 
a  final  energy  8.7  meV  which  is  just  below  Eth*  Iri  this  case  each  hole  gains  an 
average  energy  of  3.7  meV  from  the  lattice  which  means  that  energy  of  a  single  optical 
phonon  is  shared  by  about  10  holes.  This  is  possible  because  the  strength  of  the  h-h 
interaction  spreads  the  distribution  of  the  holes  quickly  in  energy,  so  that  when  a  hole 
absorbs  an  optical  phonon,  it  is  more  likely  to  share  the  gained  energy  with  an  other 
hole  than  emitting  an  optical  phonon.  This  results  in  two  holes  whose  energies  are 
below  Eq/  but  are  relatively  hotter  than  the  rest  of  the  holes,  so  that  through 
further  h-h  scatterings  the  extra  energy  will  be  distributed  among  the  holes.  When  the 
electrons  are  excited  by  1.71  and  1.81  eV  laser  pulses,  the  e-e  interaction  leads  to  a 
slower  cooling  rate.  This  happens  because  the  interaction  among  the  electrons  quickly 
redistributes  the  electrons  into  high  and  low  energy  regions.  Consequently,  many 
electrons  end  up  with  energies  below  the  phonon  emission  threshold,  while  those  which 
end  up  at  higher  energy  regions  essentially  have  the  same  probability  for  emission  of 
optical  phonons.  Those  electrons  which  have  energies  below  the  phonon  emission 
threshold  can  cool  down  only  by  giving  up  some  of  their  energies  to  other  electrons  or 
gaining  enough  energy  to  emit  an  optical  phonon.  Comparing  curve  2  in  Figures  4a  and 
5a,  we  see  that  the  manner  in  which  e-e  interaction  influences  the  cooling  rate  is  also 
dependent  on  whether  the  final  energy  after  emitting  the  LO  phonons  (assuming  no  c-c 
interactions)  is  below  or  above  E^h*  When  the  e-h  plasma  is  excited  by  excess  energy 
of  200  meV,  the  electrons  would  reach  a  steady  state  energy  of  28  meV  by  emitting  LO 
phonons  only,  which  is  not  in  equilibrium  with  the  lattice,  consequently  the  electrons 
exchange  energy  among  themselves  and  emit  LO  phonons  in  a  manner  similar  to  the 
situation  in  Figure  3a.  This  explains  why  the  cooling  curve  2  in  Figure  4a  drops  below 
curve  1  beyond  2  ps.  On  the  other  hand  curve  2  always  remains  above  curve  1  in  Figure 
5a,  and  the  electron  ensemble  reaches  a  steady  state  energy  of  19.4  meV  which  is  above 
Eth-  The  cooling  curve  for  the  holes  in  Figure  4a  behaves  in  the  same  way  the  the 
electrons  do  in  Figure  3a,  because  in  both  cases  the  initial  energies  of  the  carriers 
are  less  than  E  and  greater  than  E.,  .  When  the  initial  energy  of  the  holes  is  above 
Eg,  the  energy  exchange  among  the  holes  leads  to  faster  equilibrium  with  the  lattice 
as  can  be  seen  in  Figure  5b  (curve  2).  The  holes  retain  some  part  of  the  initial  energy 
instead  of  losing  it  to  the  lattice  via  optical  phonon  emission. 


The  case  where  only  electron-phonon  and  e-h  interactions  were  active  was  then 
investigated.  From  examining  curve  3  in  Figure  3b,  we  see  that  the  average  energy  of 
the  holes  increases  rapidly  during  the  first  picosecond  from  5  meV  to  12.5  meV  as  a 
result  of  the  energy  transferred  from  the  electrons.  Also  the  cooling  curve  for  curve  3 
in  Figure  3a  is  below  that  of  curve  2  reflecting  the  fact  that  e-h  interaction  tends 
gradually  shift  the  electron  population  to  lower  energy  states  by  transferring  the 
energy  to  the  holes.  The  cooling  rates  exhibit  the  same  trend  for  all  the  cases 
examined  as  can  be  seen  from  the  relative  positions  of  curves  2  and  3  in  Figures  3 
through  5.  Finally  the  more  realistic  situation  where  all  the  scattering  mechanisms 

e-h,  h-h,  and  carrier-phonon  are  present  was  investigated.  From  curve  4  in  Fiqure 
3a  it  is  obvious  that  the  coupling  of  e-e  and  e-h  interactions  leads  to  a  faster  cooling 
of  the  electrons  to  room  temperature.  The  same  trend  is  true  for  the  holes,  so  that  the 


energy  of  the  holes  is  intermediate  between  the  situations  where  h-h  or  e-h  interaction 
is  active  only.  However,  it  is  clear  form  examining  curves  3  and  4  in  Figures  3  through 
5,  that  the  cooling  rates  are  primarily  determined  by  the  e-h  interactions,  since  the 
cooling  curve  in  this  case  closely  follows  the  situations  v;here  only  e-h  and 
carrier-phonon  interactions  were  active. 


Conclusion 


The  c-c  interactions  are  sensitive  to  the  initial  energy  distribution  of  the  carriers 
and  they  influence  the  carrier-phonon  interactions  so  as  to  achieve  a  faster 
thermal ization  of  the  carrier  distribution.  when  the  initial  energy  of  the  carrier  is 
below  the  phonon  emission  threshold  c-c  interactions  stimulate  the  optical  phonon 
emission  or  absorption  process  depending  on  whether  the  initial  energy  of  the  carrier  is 
above  or  below  the  thermal  energy  respectively.  The  e-h  interaction  plays  a  dominant 
role  in  the  relaxation  process  even  when  excitation  energy  is  below  the  LO  phonon 
emission  threshold. 
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Figure  3.  Time  evolution  of  the  carrier  energy  for  different  active  scattering 
mechan'sms  (a)  electrons,  (b)  holes. 
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Figure  4.  Time  evolution  of  the  carrier  energy  for  different  active  scattering 
mechanisms  (a)  electrons,  (b)  holes. 
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Recently,  several  authors  [1-3]  have  investigated  the  dynamics  of 
electron-hole  plasma  generated  by  picosecond  and  sub-picosecond  laser 
pulses  in  GaAs.  These  measurements  have  indicated  the  importance  of 
electron-hole  interaction  in  the  cooling  process  (especially  in  the 
presence  of  a  cold  hole  plasma)  and  the  existance  of  non-equilibrium 
phonon  distributions  as  a  result  of  cooling  of  photoexcited  electrons  and 
holes.  Previous  theoretical  studies  of  the  cooling  process  have  ignored 
the  presence  of  the  upper  valleys  even  when  the  energy  of  the  excited 
electrons  exceeds  the  energy  separation  between  the  central  and  the  upper 
valley.  Additionally,  assumptions  had  to  be  imposed  on  the  form  of  the 
distribution  function  of  the  electrons  and  phonons.  En.semble  Monte  Carlo 
(EMC)  techniques  avoid  these  asstxmptions  and  have  been  used  to 
investigate  the  effects  of  electron-hole  (e-h)  interaction  (4]  and 
non- equilibrium  phonons  [5]  on  the  relaxation  rate  of  photoexcited 
electrons.  In  this  study  we  present  the  first  EMC  calculation  to  account 
for  both  the  e-h  interaction  and  hot  phonon  effect. 

The  Monte  Carlo  model  takes  into  account  the  e-h  ,  carrier -phonon 
interactions  and  the  non-equilibrium  LO  phonon- electron  interactions. 

The  band  model  consists  of  three  nonparabolic  valleys  for  the  conduction 
band  and  a  parabolic  heavy-hole  valence  band  presence.  The  parameters 
used  in  the  EMC  calculations  are  the  same  as  those  used  in  references 
[4]  and  [5].  Electron-electron,  electron-plasmon  and  the  screening  of 
electron-phonon  interactions  have  been  negelected.  Initially,  to  avoid 
the  complications  arising  from  the  transfer  of  electrons  to  the  upper 
valleys,  we  have  examined  the  roles  of  hot  phonons  and  electron-hole 
interaction  for  excitation  energy  of  1.8  eV.  Under  this  situation,  the 
energy  of  the  excited  electrons  is  250  meV  while  that  of  the  holes  is  40 
meV,  so  that  the  change  in  the  LO  phonon  population  is  mainly  due  the 
emission  of  LO  phonons  by  the  electrons.  The  time  evolution  of  the • 
LO-phonon  distribution  N(q)  is  calculated  as  a  function  of  the  wavevector 
q  from  the  Monte  Carlo  simulation  using  the  approach  outlined  in  [5]. 
Instead  of  recalculating  the  scattering  rates  for  the  LO-phonons  using 
the  perturbed  phonon  distribution,  we  have  used  a  Large  value  for  N(q) 
and  used  the  rejection  technique  to  accept  the  final  state  after  each 
scattering  event  involving  an  LO  phonon. 
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Ihe  relaxation  of  electrons  excited  by  a  1.8  eV  laser  pulse  was 

mnnn.r?  excitation  level  of  5x10^’^  cm'^  to  understand  the 

n  wi  Cl  e-h  interaction  and  hot  phonons  influence  the  coolinp 
process.  From  Fig.  1,  it  is  clear  that  when  the  e-h  interaction  is 
Ignored,  the  hot  phonon  effect  is  significant  only  for  times  longer 
than  1  ps.  Curve  1  in  Fig.  2  shows  that  between  0.6  and  2  ps  the 
energy  lo.,s-rate  through  the  emission  of  LX)-phonons  is  reduced  when 
t  e  phonon  heating  is  taken  into  account  compared  to  the  situation 
where  hot  phonons  are  negelcted  (curve  1).  When  the  e-h  interaction 
is  taken  into  account,  the  cooling  rate  slows  for  times  less  than  2  ps 
as  shown  by  curve  3  in  Fig.  1.  When  both  the  e-h  interaction  and  hot 
phonons  are  considered  (curve  A),  the  cooling  rate  is  slightly  faster 
for  times  below  1  ps  partly  due  the  initial  overestimation  of  the 
enhancement  of  the  phonon  emission  process  inherent  in  using  a  single 
large  value  for  N(q).  However,  it  is  obvious  that  the  cooling  process 
for  times  larger  than  2  ps  is  controlled  by  the  hot  phonons,  which 
lead  to  a  slower  cooling  rate.  Additionally,  no  significant  change  in 
the  energy  of  the  photoexcited  holes  was  noticed,  when  the 
electron-hot-phonon  interaction  was  considered.  The  situation  where 
hole-hot-phonon  interaction  is  taken  into  account  might  lead  to  an 
increase  in  the  hole  energy  as  a  result  of  the  phonon  reabsorption 
process.  This  work  is  underway  and  the  results  will  be  presented 
elsewhere . 


Fig.  1.  Time  evolution  of  the 
electron  energy:  1)  e-ph, 

2)  e-ph  (hot),  3)  e-ph  +  e-h, 
4)  e-ph  (hot)  +  e-h. 


Fig.  2.  Time  dependence  of  the 
electron-energy-loss  rate:  1) 
w/o  hot  phonons,  2)  with  hot 
phonons . 


The  situation  where  the  electrons  were  transferred  to  the  upper 
valleys,  was  investigated  for  the  case  of  excitation  by  a  588  run  laser 
pulse  using  EMC  and  taking  into  account  the  carrier-phonon 
interactions  only.  The  fraction  of  electrons  in  the  central  valley 
is  shown  in  Fig.  3  as  a  function  of  the  time  after  the  application  of 
the  pulse.  From  Fig.  3,  it  is  obvious  that  the  simple  picture  of 
electrons  cascading  down  the  central  valley  by  emitting  LO-phonons  is 
not  complete.  The  electrons  rather  transfer  to  the  upper  valleys  and 
cool  down  in  these  valleys  by  emitting  phonons  through  polar-optical 
and  deformation  potential  interactions  as  shown  in  Fig.  4  (curves  2 
and  4).  The  energy  loss  rate  through  LO  phonon  emission  is  dominated 
by  the  central  valley  electrons  for  times  greater  than  1  ps  and  reach 
a  maximum  after  2  ps.  Note  that  after  2  ps  about  80%  of  the  electrons 
are  in  the  central  valley.  The  energy  of  the  electrons  in  the  central 
valley  (cuirve  1  in  Fig.  5)  remains  above  100  meV  up  to  3  ps  after  the 
excitation  and  an  increase  in  the  average  energy  occurs  at  2  ps  due  to 
the  return  of  the  electrons  to  the  central  valley.  The  initial  rapid 
drop  in  the  average  kinetic  energy  of  the  electron  ensemble  (curve  2 
in  Fig.  5)  is  due  to  the  transfer  of  a  large  fraction  of  the  electrons 
to  the  upper  valleys  as  can  be  seen  in  Fig.  3.  Curve  3  in  Fig.  5 
indicates  that  the  cooling  of  the  electron  in  the  L- valley  is  not 
negligible . 


Fig.  3.  Time  evolution  of  the 
fraction  of  electrons  in  the 
central  valley. 


Fig.  4.  Time  dependence  of  the 
electron-energy- loss  rate; 

1)  LO  (F).  2)  LO  (L).3)  total  LO 
4)  Intervalley  scattering 
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Fig.  5.  Mean  kinetic  energy  of 
the  electrons  in;  1)  central 
valley,  2)  central  +  L- valley, 
3)  L-valley. 


In  conclusion,  have  shown  that  the  cooling  process  of  the 
photoexcited  electrons  is  controlled  by  e-h  interaction  for  times  less 
than  1  ps  and  by  hot  phonons  for  times  greater  by  2  ps.  In  addition, 
we  have  demonstrated  that  the  cooling  of  the  electrons  in  the  upper 
valleys  are  important  for  high  excitation  energies. 
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AbPtract 

Over  the  past  decade  a  variety  of  physical  phenomena  have  emerged  which  are  exercising 
profound  constraints  on  the  speed  of  devices  whose  primary  transient  is  electronic  in 
origin.  For  example,  the  broad  high  frequency  extrapolations  of  high  speed  arising  from 
velocity  overshoot  were  found  to  be  limited  by  the  effects  of  space  charge  and  the 
imposition  of  fields  that  take  a  finite  time  to  reach  their  specified  level.  Much  of  the 
latter  problem  may  be  mlnimited  through  optical  processes,  but  in  its  place  rests  the 
ccnstr.'lnta  of  electron-hole  Interacti''"  and  other  that  may  require  quantum 

coherence  for  iiipl«i«..ntation .  In  the  case  of  electron-hole  Interaction  it  has  been  found 
that  it  provides  an  important  channel  through  which  electron  energy  is  transferred  to  the 
lattice,  and  that  the  electron  cooling  rates  at  high  hole  concentrations  are  higher  when 
the  electron-hole  interaction  is  included,  than  when  it  is  ignored.  At  low  hole 
concentrations  the  reverse  is  true.  The  effect  of  electron-hole  interaction  also 
influences  velocity  overshoot.  In  the  case  of  quantum  coherence,  of  interest  is  the 
effect,  e.g.,  of  temporal  scatterers  (phonons)  on  the  time  dependent  wave  functions  and 
modifications  of  Interference  effects.  The  influence  of  space  charge  on  transient  velocity 
overshoot,  the  effects  of  electron-hole  interaction  on  the  relaxation  of  an  electron-hole 
plasma,  and  the  effects  of  time  dependent  and  time  independent  scattering  on  quantum 
interference  will  be  discussed. 


Introduction 

The  early  suggestions  of  Ruch^  and  subsequent  calculations  of  others,  have  suggested 
that  high  mobility  semiconductors  were  capable  of  extraordinary  high  speeds  (10®cm/sec) 
provided  ^ey  were  sampled  for  sufficiently  short  time  intervals.  For  submicron  and 

devices,  this  translates  Into  time  scales  of  0.1  ps,  and  cutoff  frequencies 
these  predictions  are  clearly  Impressive,  of  serious  concern  is  an 
Inability  of  realizing  the  modest  goal  of  room  temperature  operation  at  lOOGHz,  and  speeds 
of  less  than  5ps,  within  the  framework  of  a  three  terminal  device.  The  difficulty  of 
achieving  anything  close  to  that  offered  by  the  Ruch  calculations  has  been  discussed  by  a 
number  of  workers  over  the  past  several  years.  Perhaps  the  simplest  statement  summarizing 
these  studies  is  that  Ruch's  assumptions  are  not  realized  experimentally.  These 
assumptions  include;  1)  uniform  fields  within  the  device,  2)  zero  rise  time  of  the  electric 
field,  and  3)  neglect  of  displacement  currents.  Because  of  the  importance  of  these 
assumptions,  each  will  be  briefly  reviewed,  along  with  the  Influence  of  nonequilibrium 
holes  and  phonons.  Additionally  the  Influence  of  quantum  mechanical  constraints  on 
transport  will  be  discussed. 

Honuniform  fields 

Electric  fields  in  devices  are  nonuniform,  as  illustrated  for  the  three  terminal 
permeable  base  transistor  (PBT) .  The  PBT  as  currently  configured,  is  a  micron-length 
structure  with  ultra-submicron  features.  A  section  of  the  PBT  is  shown  in  figure  1,  along 
'^ith  its  an  equivalent  planar  FET.  As  shown  in  figure  2^,  for  a  bias,  of  one  volt  between 
the  emitter  and  the  collector  the  potential  distribution  indicates  the  expected  field 
non-uniformity.  For  example,  near  the  cathode,  the  field  is  low,  increasing  to  a  value 
near  20  )cv/cm  and  finally  approaching  a  low  value  near  the  collector  contact.  In  uniform 
field  studies  the  predicted  peak  velocity  approaches  10®ci^sec,  followed  by  a  decay  to 
steady  state.  The  peak  velocity  is  achieved  with  a  carrier  concentration  that  is  close  to 
its  equilibrium  value.  Decay  is  dominated  by  electron  transfer,  with  a  weaker  co^>onent 
due  to  momentum  relaxation.  In  steady  state  the  gamma  valley  carrier  velocity  is  high, 
while  the  concentration  of  gamma  valley  carriers  is  low.  In  the  PBT  the  initial  high  peak 
in  current  associated  with  velocity  overshoot  is  not  achieved.  Rather  the  high  velocity  is 
often  achieved  at  the  expense  of  carrier  concentration.  This  la  displayed  in  figures  3  and 
4  which  show  the  distributon  of  carrier  density  and  velocity  for  the  PBT.  It  is  worthwhile 
noting  that  while  the  distribution  of  charge  and  velocity  are  not  as  encouraging  as  that 
from  the  uniform  field  calculations,  the  current  voltage  levels  are  between  a  factor  of  3 
and  4  above  that  obtained  from  the  drift  and  diffusion  equation  calculations. 
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£initg  fisc  tjpe  contributions 

A  second  issue  is  the  zero  rise  tine.  In  virtually  all  transient  transport 
calculations,  overshoot  phenomena  is  obtained  by  examining  the  resjxjnse  of  carriers  to  an 
electric  field  that  is  suddenly  turned  on.  Under  conditions  of  finite  rise  time, 
relaxation  effects  are  coupled  to  the  increasing  field  and  the  peak  velocity  is  diminished, 
as  shown  in  figure  He  note  that  there  is  a  qualitative  similarity  between  the  finite 

rise  time  anu  the  response  of  carriers  to  a  gradual  spatially-varying  electric  field. 

Qlgplacement  current  effects 

To  carry  the  discussion  of  non-uniform  fields  even  further,  calculations  have  been 
performed  for  the  paradigm  two-terminal  ir'Tfir*'  submicron  structure  shown  in 
figure  6^.  The  structure  results  in  non-uniform  fields  and  when  transients  are  computed 
it  is  a  current  transient  that  is  involved  rather  than  a  velocity  transient.  Indeed,  in 
virtually  all  structures  it  is  a  current  transient  rather  than  a  velocity  transient  that  is 
measured.  As  shown  in  figure  7,  there  is  a  peak  value  to  tho  current  that  occurs  within 
the  time  scale  over  which  velocity  overshoot  may  be  expected.  However,  the  details  of  the 
transient  show  that  the  initial  time  dependent  behavior  is  associated  entirely  with  field 
rearrangement,  and  that  the  transient  is  entirely  a  displacement  current  transient^. 

Hanegulllbrlua  electron  hole  phonon  interaction 

In  overshoot  studies  most  inelastic  collisions  are  treated  assuming  the  emission  of  a 
single  phonon  and  that  the  emitted  phonon  is  in  equilibrium  with  the  lattice.  For 
energetic  electrons,  a  sequence  of  as  many  as  ten  to  fifteen  phonons  maybe  emitted  as  the 
electron  is  relaxing,  and  as  a  result  an  excess  number  of  phonons  may  result  in  a  lack  of 
phonon  equilibrium  with  the  lattice.  This  phenomena  is  illustrated  for  an  unbiased  slab  of 
gallium  arsenide  subject  to  a  uniform  photoexcited  laser  pulse  of  energy  1.8  ev.  Under 
'^is  situation  holes  as  well  as  electrons  are  excited.  The  energy  of  the  excited  electrons 
is  250  Mev,  while  that  of  the  holes  is  40  Mev.  The  time  evolution  of  the  resulting 
electron,  hole  and  phonon  distribution  is  calculated  using  Monte  Carlo  techniques.  Figure 
8  shows  the  relaxation  of  the  electrons  to  equilibrium  with  and  without  hot  phonons.  As 
seen,  the  inclusion  of  the  hot  phonons  results  in  a  delay  of  the  relaxation  to  steady  state. 

The  situation  in  the  presence  of  finite  fields  and  equilbrium  phonons  is  discussed  next 
and  shown  in  figures  9  and  10®.  These  calculations  were  performed  assuming 
nonequilibrium  interactions  between  the  gamma  valley  electrons  and  holes.  L-valley  hole 
interaction  was  not  examined.  The  first  point  of  note  is  that  the  electron  hole 
interaction  is  an  inelastic  collision.  In  addition,  since  both  electrons  and  holes  are 
responding  to  the  presence  of  the  applied  fields,  momentum  scattering  of  the  electrons  and 
holes  can  result  in  a  reduction  of  the  electron  velocity,  and  in  some  cases  turn  the 
electron  in  the  "wrong"  direction.  The  net  result  is  a  decrease  in  the  peak  electron 
velocity.  It  is  impo.'tant  to  note  that  while  the  peak  velocity  is  lower  for  the 
calculation  with  the  electron-hole  interaction,  there  is  a  cross-over  where,  for  a  certain 
interval  of  time,  the  carrier  velocity  in  the  presence  of  holes  exceeds  that  obtained  in 
holes.  This  arises  because  the  electron  hole  interaction  is  an  efficient 
netnod  or  energy  exchange  and  ae  a  result,  the  electrons  are  not  energetic  enough  to 
transfer  into  the  L-valley,  In  other  words,  there  is  a  greater  fraction  of  electrons  that 
are  retained  in  the  gamma  valley  when  holes  are  included  than  when  they  are  not  Included. 

The  possibility  exists  that  the  velocity  of  injected  electrons  in  p-typo  material  may,  for 
a  range  of  values  of  electric  field,  be  higher  than  that  of  thermally  generated  carriers  in 
«-type  material. 

Quantum  contributions 

It  will  be  noticed  that  all  of  our  discussion  has  been  confined  for  electric  fields  and 
It  been  presumed  that  we  are  dealing  with  a  gauge  Invariant  system, 
and  that  the  choice  of  gauge  Is  not  relevant.  As  discussed  by  Aharanov  and  Bohm  and®, 

Me  vector  potential  can  lead  to  observables  such  as  interference  patterns  for  electrons 

traverse  ostensible  field-free  regions,  in  other  words.  Interference  patterns  are 
accounted  for  in  terms  of  differences  in  the  vector  potentials  encountered  by  electrons  on 
two  opposite  sides  of  a  flux  line.  To  see  how  these  effects  enter  into  the  problem  it  is 
useful  to  re-examine  an  approach  to  quantum  transport  introduced  earlier^. 

In  reference  7,  which  Introduced  the  Hiqner  function  as  an  approach  to  quantum  transport 
in  electron  devices,  a  set  of  hydrodynamic  balance  equations  was  introduced.  While  the 
nyarodynamic  balance  equations  were  developed  for  a  system  of  particles,  the  equations  may 
be  reduced  to  examine  transport  for  carriers  not  subject  to  ecattering.  Thus,  the  single 
wave  functions  are  of  interest. 
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Deteralnlng  tha  slngla  pairticla  wave  functlona  ragulras  so]<'tlon8  to  the  single  particle 
Schrodlnger's  equation.  Rewriting  Schrodlnger'e  equation  in  terns  of  a  hydrodynanic  aodel 
requires  solutions  to  the  following  two  equations: 


ae 

-  4  divj  -  0 

dl 


(1) 


du 

m -  -  -9(Q*V) 

dl 


(2) 


Where  V(x)  Is  a  classical  scaler  potential  and  Q  is  the  quantum  potential 


^2 

2m 


(3) 


The  reference  7  study  suggests  that  both  €  »  v'Vand  J-6u  be  .-egarded  as  the  basic 
physical  quantities  representing  properties  of  the  continuous  fluid.  It  is  noted  that  the 
vector  potential  has  dropped  out  of  the  basic  equations,  which  have  a  completely  classical 
form. 


Phlllppldls,  Boho  and  Kaye  (PBK)^  argued  that  the  above  equations  were  incomplete. 
The  most  important  objection  arose  because  certain  subsidiary  conditions  were  not 
satisfied.  To  see  this,  the  wave  function  for  the  system  is  written  as  y -RexpliS/h]  , 
with  the  carrier  velocity  expressed  as 


u 


VS  e 
m  me 


(4) 


If  it  is  now  required  that  the  wave  function  be  single  valued,  then  we  are  forced  to 
admit  the  condition 

^  Vs- dx  -  ^  (mu  ♦ a)  dx  -  nh  (5) 


This  means  that  the  observable  movements  of  electron  transport  in  a  device  must 
necessarily  be  restricted  in  such  a  way  that  the  integral  of  its  velocity  around  a  circuit 
containing  a  flux  line  depends  upon  the  flux  within  the  circuit.  AB  argue  that  there  is  no 
intuitively  clear  justification  for  the  equation  (8)  constraint.  Rather,  they  propose  to 
assume  that  the  vector  potential  has  a  physical  meaning. 

Figure  11  is  a  plot  of  the  quantum  potential  for  an  electron  passing  through  a  dual  slit 
in  the  absence  of  a  magnetic  field.  The  smooth  parabolic  hills  are  in  the  Immediate 
vicinity  of  the  slit.  Particles  subject  to  such  a  potential  accelerate  smoothly  in  the^ 
forward  and  sideward  direction.  This  can  be  seen  from  the  family  of  velocity  trajectories 
displayed  in  figure  12.  Mote,  near  the  slits  where  the  wave  functions  are  localized 
packets,  the  quemtum  potential  exhibits  two  broad  peaks.  As  the  wave  packets  recede  and 
spread,  they  begin  to  overlap  significantly  and  their  interference  properties  come  into 
play.  The  shape  of  the  qu^ultum  potential,  a  long  way  from  the  slits,  is  determined  by  the 
boundary  conditions,  the  shape  and  size  of  the  slits.  It  is  through  the  force  exerted  by 
the  quantum  potential  on  the  particle  that  the  boundary  conditions  are  made  physically 
significant  tar  from  the  slits.  For  the  case  when  a  flux  line  is  place  between  the  slits, 
the  quantum  potential  is  altered,  as  shown  in  figure  13.  The  velocity  trajectories  are 
altered  as  shown  in  figure  14.  The  important  point  to  notice  is  that  the  trajectories  are 
displaced. 

The  above  discussion  reached  Initial  device  Implementation  by  Datta,  et  al.^,  who 
conceptually  altered  the  path  length  of  a  beam  of  particles  as  shown  in  figure  15,  where 
the  slit  in  the  CaAs  between  the  AlGaAs  is  exaggerated.  Additionally  they  altered  the 
phase  through  application  of  a  magnetic  field.  These  alterations  lead  to  interference  in 
the  downstream  carrier  wave  functions.  This  result,  as  discussed  in  reference  9,  is  the 
appearance  of  oscillations  in  the  conductance  of  the  system. 
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An  iaiportant  point  not  addressad  in  refer«nc«  9  is  th«  tamporal  duration  of  th« 
conductance  oscillations.  It  nust  be  considered  that  the  sporadic  appearance  of  phonons 
will  result  in  a  loss  of  coherence  and  limit  constructive  interference  to  time  scales 
shorter  than  100  femtoseconds.  For  time  scales  longer  than  100  femtoseconds,  the 
constructive  interference  effects  will  certainly  be  diminished.  Whether  they  will 
disappear  is  a  matter  to  be  determined.  Here  we  note  that  coherence  is  essential  for 
quantum  well  devices.  There  are  other  issues  that  may  tend  to  mitigate  the  appearance  of 
the  AB  effect  such  as  elastic  scattering  arising  from  surface  roughness. 

Conclusions 

The  driving  force  toward  high  speed  devices  is  limited  by  the  presence  of  carrier- 
carrier  scattering,  phonons  and  submicron  feature  sizes  that  tend  to  reduce  the  coherence 
of  scattered  waves.  While  these  effects  do  not  lessen  the  importance  of  the  noneguillbrium 
transients,  unless  these  effects  are  incorporated  into  the  design  of  the  device, 
applications  may  be  limited. 
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Figure  8.  Energy  relaxation  with  equilibrium 
phonons  (1)  and  nonequilibrium 
phonons  (2) . 
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Figure  9.  Overshoot  witli  and  without 
e-h  scattering.  (From  (51). 


Figure  10.  Carrier  distribution  with  and  without 
e-h  scattering.  (From  (51). 
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Fjfjuro  15.  Proposed  Structure  for  quantum 

interference  effects.  (From  (9)). 


Figure  14.  Trajectories  for  the  AB  effect 
(note  the  shift  of  the  overall 
pattern) .  (From  (0)) . 


